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Dr. George Ellery Hale, distinguished and much beloved astrono- 
mer, died in Pasadena, California, on February 21 after an illness 
extending over several years. A scientist who made many remark- 
able contributions to his chosen field of solar research, the founder 
and almost throughout his life an editor of the Astrophysical Journal, 
founder and first director of the Yerkes and the Mount Wilson ob- 
servatories, organizer of the International Union for Co-operation 
in Solar Research and of its successor, the International Astronomi- 
cal Union, for many years foreign secretary of the National Academy 
of Sciences, and the founder of the National Research Council, 
Dr. Hale has left an impress upon the life of science which will 
endure for generations to come. In later years, when illness pre- 
vented him from carrying on actively the investigations he had most 
at heart, he conceived the project of the 200-inch telescope and se- 
cured the necessary financial support, leaving this great instrument 
as his final gift to the science he had loved so well. 

The name of Hale is associated with so many important organiza- 
tions in the creation and development of which he took the leading 
part that it is not always easy to realize that the instinct for re- 
search was the dominant feature of his life. He was a born investi- 
gator and was never so thoroughly content as when he could work 
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uninterruptedly at the telescope and spectrograph or could analyze 
his results in the quiet of his office or laboratory. For this reason 
there can be little doubt that the happiest years of his life were 
those at Kenwood and the earlier years at the Yerkes and Mount 
Wilson observatories. They were periods of construction and 
development and gave full scope to his creative and experimental 
ability before the heavy responsibilities of the larger organizations 
had begun to be felt to any marked degree. But Hale also realized, 
to an extent much greater than most investigators, that the major 
problems of physical science are usually so extensive and complex 
that progress in their solution requires the co-operative effort of 
many individuals and many institutions. So, although by nature 
a strong individualist in his scientific work and in his outlook on 
life, he was constantly engaged in organizing plans for co-operative 
observations, in encouraging young investigators, and in securing 
the establishment of new institutions and wider and better facilities 
for research. 

Of the conduct of his scientific work it would hardly be possible 
for any two, even of Hale’s closest, associates to write in the same 
way, such were the characteristics of this gifted and many-sided 
man. His active and brilliant mind was constantly at work analyz- 
ing his problems, weighing new hypotheses, and designing new meth- 
ods of attack, always with an enthusiasm resembling the creative 
joy of the artist. He had to a remarkable degree the gift of being 
able to distinguish the essential and salient features of an investiga- 
tion; and it was rarely, indeed, that even a minor research which he 
attempted failed to bring productive results. One very striking char- 
acteristic was his interest in instruments and his remarkable skill 
in the designing of new apparatus. The invention of the spectro- 
hetiograph and the spectrohelioscope and the adaptation of the 
tower telescope and the vertical spectrograph for solar investigation 
are excellent examples of this quality of Hale’s mind. In part it 
was probably due to his engineering training, but much more to 
his appreciation of the enormous importance of the value of im- 
proved instruments to astronomical progress, and to his creative 
ability and love of fine workmanship. He took immense pleasure 
in seeing a new instrument assume form; and, as he often said, 
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“had I not been an astronomer, I should have liked to have been 
an instrument-maker.”’ 

The early years of Hale’s scientific life and development form an 
exceedingly interesting story. We may almost say that, as in the 
case of Michelson, the phenomena of light fascinated him from the 
beginning. From childhood the microscope and telescope, the colors 
of the spectrum, and the infinite variety of forms in which light 
manifests itself in everyday life were an unfailing source of joy and 
satisfaction to him. Probably no boy ever used a diffraction grating 
at so early an age, and the story of Hale’s first meeting with Row- 
land and the latter’s remark that “any grating of mine should be 
good enough for an infant” was always quoted with much amuse- 
ment by Hale in after-life. His father gave him constant encourage- 
ment in these interests, and gradually the equipment was built up 
which came to form the Kenwood Observatory. 

Fortunately, we have in Hale’s own words a most interesting 
statement regarding this formative period of his life and some of 
the influences which led him to take up astrophysical research as 
a career in after years. The statement is contained in some bio- 
graphical notes, written by him in 1933, and represents his mature 
judgment on the factors which affected most greatly the direction 
of his intellectual interests. The notes also throw an amusing side- 
light on one or two astronomical worthies of this period. 

I cannot fix the date of my first interest in astronomy, but it must have been 
when I was thirteen or fourteen years old. I built a telescope; but as I used a 
large single lens, the images were not good. About this time I became acquainted 
with S. W. Burnham, then a stenographer in the Chicago law courts by day and 
an ardent observer of double stars by night. Through him I learned of a second- 
hand Clark refractor, of 4 inches aperture. This was purchased by my father, 
and I mounted it on the roof of the house. The astonishing views it afforded of 
Saturn, the moon, Jupiter, and other objects excited an intense desire to carry 
on actual research. So I attached a plateholder to the telescope and photo- 
graphed a partial eclipse of the sun. I also began to observe sunspots, and made 
drawings of them. Thus I became an amateur astronomer. 

At this period the 18}-inch Dearborn refractor stood at the summit of a 
lofty tower forming part of the castellated building of the old University of 
Chicago (long since demolished, after the removal of the telescope to the 
Northwestern University at Evanston). Professor G. W. Hough, known for his 
observations of Jupiter and double stars, who was in charge of the instrument, 
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was kind enough to let me look through it frequently. I was also permitted to 
aid the feeble gas-engine in turning the dome, or rather drum, of the observa- 
tory. Naturally, my own ambitions were thus stimulated, but neither Hough 
nor Burnham had the slightest interest in astrophysical research, and I could 
not have devoted my life to such work as they were doing, valuable as it was 
to science. The reason lay in the fact that I was born an experimentalist, and I 
was bound to find the way of combining physics and chemistry with astronomy. 
Fortunately, it was not far to seek. 

In Cassell’s book is a description of a spectroscope, with an account of the 
construction of a simple instrument, using either a luster prism or a hollow one 
filled with carbon disulphide. I built both forms, and the odor of the disulphide 
abides with me after a lapse of fifty years. Then my father, always ready to 
encourage serious efforts, enabled me to buy a small spectrometer. This had a 
single prism, and I lost no time in fitting it with a small plane grating. For by 
this time I had learned of the work of Rutherfurd and Rowland and had ac- 
quired some slight conception of the possibilities of high dispersion. Nothing 
could exceed my enthusiasm in observing the solar spectrum and in measuring 
the principal lines. I bought Lockyer’s Studies in Spectrum Analysis and began 
the observation of flame and spark spectra and their comparison with the spec- 
trum of the sun. At last I had found my true course, and I have held to it ever 
since. 

In passing, however, I must emphasize the advantages I gained from my 
earlier work in other fields. It gave me a broad interest in many branches of 
science and taught me something of their mutual dependence. It led me to 
read The Origin of Species and initiated a lasting desire to study evolution. It 
taught me to regard the sun as a typical star, a link in the long evolutional chain, 
and thus helped me to avoid becoming exclusively a specialist in solar research. 
And it was the source of all the work I have done in helping to develop such 
institutions as the National Academy of Sciences, the California Institute of 
Technology, the National Research Council, and the International Research 
Council (now the International Council of Scientific Unions). I still look with 
delight at rotifers and other microscopic objects, and realize that the marvels 
I saw through the little microscope of my boyhood days are as astonishing as 
any revealed in the heavens by the largest of telescopes. 


The early education which Hale received in Chicago was at the 
Oakland Public School and later at Allen Academy, where he was 
allowed to care for the “philosophical instruments” and to assist 
in demonstrations before the student body. At the suggestion of 
Daniel H. Burnham, the well-known architect, who was a close 
friend of Hale’s father, the Massachusetts Institute of Technology 
was chosen for the college career of the son. Hale entered the Insti- 
tute in the autumn of 1886 and was graduated in 1890. His courses 
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were mainly in mathematics, physics, and chemistry, and for astron- 
omy he was obliged to go elsewhere. Through special arrangement 
with E. C. Pickering he was able to serve as volunteer assistant at 
the Harvard College Observatory, where, among other experimental 
instruments, he tried his first form of spectroheliograph. 

Hale’s marriage to Miss Evelina Conklin occurred soon after he 
left the Massachusetts Institute, and on their western trip together 
they visited the Lick Observatory, where Hale was greatly im- 
pressed by the 36-inch telescope and the admirable spectroscopic 
investigations of planetary nebulae by Keeler. On his return to 
Chicago he passed three active years of research at the Kenwood 
Observatory, interrupted only by a trip in 1891 to Europe, where 
he met many astronomers and spectroscopists and studied their 
methods and equipment. The winter of 1893-94 he spent in study 
at the University of Berlin, where the lectures by Planck on thermo- 
dynamics and theories of radiation were of great value to him in 
his later work. He made frequent visits to Potsdam and became 
thoroughly familiar with the researches of Vogel, Scheiner, and 
others. It was during this European journey that he completed his 
plans for the establishment of the Astrophysical Journal. 

After his stay in Berlin and an expedition to Mount Aetna, 
where he tried out his apparatus for photographing the corona 
without an eclipse, he returned to Chicago, where he continued his 
research at the Kenwood Observatory but was also actively engaged 
in plans for the larger institution which had taken form as the 
Yerkes Observatory. 

The University of Chicago was at this time in its great period 
of development under President Harper. A remarkable judge of 
men, he at once appreciated the abilities of this gifted student and 
investigator of twenty-four, and Hale became a member of the staff 
of the University as associate professor of astrophysics in 1892. To- 
gether they formed plans for a great observatory; and, both being 
men of initiative and constructive imagination and enthusiasm, they 
succeeded in interesting Mr. Yerkes of Chicago, who made the gift 
for the project in the same year. Hale carried on the detailed plans 
with characteristic energy, outlining the investigations to be under- 
taken, designing instruments, and assembling a staff. Barnard, with 
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his wide experience and extraordinary skill in visual observations, 
came from Mount Hamilton; Burnham, also with Lick Observatory 
training, from Chicago; Wadsworth, from Michelson’s laboratory; 
Frost, from Dartmouth College; and Ellerman, with his background 
of years of experience with Hale, from Kenwood. The Yerkes Ob- 
servatory was completed and dedicated in 1897, with Hale as its 
first director and professor of astrophysics at the University of 
which it formed a part. He retained his relationship to the Univer- 
sity of Chicago even after he came permanently to California, being 
nonresident professor from 1905 onward. 

The winter of 1903-4, owing to illness in his family, Hale spent 
in Southern California. His attention had previously been attracted 
to Mount Wilson by Professor Hussey, who was investigating pos- 
sible sites for a high-altitude solar observatory in behalf of a com- 
mittee appointed by the Carnegie Institution of Washington, of 
which Hale was a member. In June, 1903, accompanied by Dr. W. 
W. Campbell, Hale made his first visit to Mount Wilson, where 
Hussey was still continuing his observations. Tests of the solar see- 
ing were made by Hale at this time and during several of the suc- 
ceeding winter months. They proved so satisfactory that applica- 
tion was made to the Carnegie Institution for a grant to finance an 
expedition from the Yerkes Observatory and the installation of the 
Snow horizontal telescope. A group consisting of Ellerman, Ritchey, 
and Adams came in the spring of 1904, followed somewhat later by 
F. G. Pease, and formed the nucleus of the staff of the future Mount 
Wilson Observatory. Throughout 1904, systematic records of both 
solar and stellar observing conditions were maintained; and, after 
much careful consideration, Hale recommended to the Carnegie 
Institution the establishment of an independent observatory. A 
grant for this purpose was made in December, 1904; the Snow tele- 
scope was purchased from the Yerkes Observatory; and the Mount 
Wilson Solar Observatory, later changed to the Mount Wilson Ob- 
servatory, as its stellar activities increased, was organized. 

This account of the origin of the Mount Wilson Observatory is 
given somewhat in detail to illustrate the caution and painstaking 
care which Hale combined with his enthusiasm in the development 
of his major projects. Probably no observatory was ever established 
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after so thorough an examination of observing conditions and of 
all the physical considerations entering into its construction, main- 
tenance, and growth. Questions of transportation, water supply, 
electric power, building materials, and living conditions were all 
investigated by him personally; and no subject was too small to be 
given thought and consideration. It was a period in which he en- 
joyed life enormously. The isolation of the mountain top in these 
early years made a strong appeal to the pioneering spirit latent in 
him, and his joy in discoveries of the everchanging beauties of na- 
ture about him was second only to that of his discoveries in science. 

In 110, after six years of active and most productive work, the 
first shadows of physical illness began to fall across Hale’s life. After 
this time he was never completely free from more or less severe 
attacks of a type of brain congestion which occasionally produced 
severe pain and almost complete mental exhaustion. He soon 
learned that such attacks were most apt to occur when his mind 
was occupied with his scientific work, and with profound courage 
and at great sacrifice he attempted to adjust himself to the situation. 
Without completely giving up his research, he devoted more time 
to writing, to executive work, and to the organization and develop- 
ment of projects he had especially at heart. After the World War 
broke out, he spent much time in Washington, where, at the re- 
quest of President Wilson to the National Academy of Sciences, 
he was mainly responsible for the creation and organization of the 
National Research Council. He was foreign secretary of the Na- 
tional Academy from 1g1o to 1921 and gave much thought to plans 
for making the Academy a more active force in national and inter- 
national science. He also visited Europe several times to attend 
meetings of the International Research Council, afterward the Inter- 
national Council of Scientific Societies, and aided greatly in directing 
its policies. Many of these interests he maintained to the end of 
his life. 

On his return to the Mount Wilson Observatory, Hale continued 
his scientific work as his strength permitted, but with increasing 
difficulty; and in 1923 he resigned the directorship. He was ap- 
pointed honorary director by the Carnegie Institution; and, with full 
realization of his physical limitations and with remarkable courage, 
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he faced the problem of how to make the most effective use of such 
years as remained to him. He solved this problem in two ways. 
First, he constructed, largely at his own expense, the small, well- 
equipped solar laboratory which he presented to the Carnegie Insti- 
tution, and in which he could carry on such research as his health 
allowed. Second, he undertook the writing of a number of articles 
planned to bring some of the results of astronomical discovery before 
the intelligent public. In this purpose he was remarkably success- 
ful, his clarity and simplicity of statement forming an admirable 
setting for the spirit of enthusiasm pervading all his writing. 

Two striking, but widely dissimilar, accomplishments mark these 
later years of Hale’s life. The first was the invention of the spectro- 
helioscope, which supplements in many important respects the 
spectroheliograph devised by him so many years ago at Kenwood. 
It is already yielding results of great value at stations all over the 
world and, through the application of moving-picture attachments, 
gives promise of proving the most efficient instrument yet designed 
for the study of rapidly changing phenomena on the sun’s surface. 
The second was the project of the 200-inch telescope. Attracted by 
his article published in 1928, ‘‘The Possibilities of Large Telescopes,” 
and strongly influenced by their personal knowledge of Hale and 
his contributions to science, the president and trustees of the Gen- 
eral Education Board made a grant for the construction of a 200- 
inch telescope, the instrument being given to the California Institute 
of Technology on the basis of a plan of co-operation with the Car- 
negie Institution and the Mount Wilson Observatory. From 1929 
until the closing days of his life, problems relating to this great in- 
strument occupied most of Hale’s thoughts. He was chairman of 
the group having the project in charge, and his active mind was 
fertile in suggestions for the design both of the telescope itself and of 
the equally important auxiliary instruments. As illness increased 
and his strength gradually failed, his mind turned constantly to the 
progress of the telescope and the observatory; and each clear day, 
as the winter rains during the past winter were delayed, he rejoiced 
in ‘one more day for work on Palomar.” 

The first article published by Hale appeared in 1889 in a small 
Chicago periodical, The Beacon, and was entitled “Stellar Photog- 
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raphy.’ This was written at the age of twenty. Succeeding papers 
during this period appeared in The Beacon, Astronomische Nach- 
richten, American Journal of Science, Sidereal Messenger, Publica- 
tions of the Astronomical Society of the Pacific, and especially 
Astronomy and Astrophysics, the astrophysical portion of which he 
edited during the three years of its existence. He had for several 
years considered the desirability of establishing a journal devoted 
to spectroscopy and astrophysics and had sought carefully the opin- 
ions of scientists from all parts of the world. Assured of their strong 
support, Hale founded in 1895, in collaboration with James E. 
Keeler, the Astrophysical Journal, the University of Chicago assum- 
ing the responsibilities of publication and taking over, through pur- 
chase, Astronomy and Astrophysics. The new journal was successful 
from the beginning and was truly international in character, with 
a board of collaborators whose names were eminent in science 
throughout the United States and Europe. Many of the classical 
papers of astrophysics and spectroscopy have first appeared in the 
Astrophysical Journal, and the influence of its founder and first 
editor has been maintained throughout the period of its existence, 
now nearly half a century in length. 

Nothing illustrates more clearly Hale’s breadth of interests than 
a study of his bibliography. In a list of nearly four hundred and 
fifty titles there is found a range from spectroscopic observations of 
solar prominences to an analysis of the intellectual culture of France. 
Included in it are his two interesting nontechnical books, Stellar 
Evolution and Ten Years’ Work of a Mountain Observatory, together 
with compilations, in book form, of his numerous magazine articles 
of an interpretive character. The longer technical monographs ap- 
peared in the publications of the Yerkes Observatory and the 
Carnegie Institution of Washington, but most of his best-known 
papers are found in the Astrophysical Journal. About go per cent 
of his published articles deal with the results of his scientific re- 
searches. 

The earliest work which Hale carried on at the Kenwood Ob- 
servatory consisted of visual observations of solar prominences. 
From the outset he realized the immense value of photography in 
affording a complete and permanent record of these objects, and 
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his ingenious mind immediately undertook to accomplish this pur- 
pose. The problem was twofold in character. First it was neces- 
sary to find a spectral line, other than the hydrogen Ha used for 
visual observations, through which the prominences could be photo- 
graphed, since the photographic plates of this period had little 
sensitiveness to red light. Second, an instrument must be devised 
by means of which the small portion of the image of the prominence 
or the sun’s disk which appeared in the light of the spectral line 
could be made to build up the complete image. Hale solved the 
first part of the problem through his investigations of the H and K 
lines of calcium, already known through the work of C. A. Young, 
to be always bright in the chromosphere and prominences, as well 
as over sunspots. These lines, which Hale was the first to identify 
with spark lines of calcium, were used for most of his work with the 
spectroheliograph until 1908, when photographic plates more sensi- 
tive to red light became available and the Ha line could also be used 
to advantage. It is probable, however, that he never felt inclined 
to revise his early statement about the H and K lines, ‘‘which to 
my mind form the most interesting group in the solar spectrum.”’ 
Two instrumental devices for obtaining photographic records of 
prominences and other solar phenomena were invented by Hale as 
early as 1889. In the second of these, which he finally adopted at 
the Kenwood Observatory, and to which he gave the name of 
“‘spectroheliograph,” two movable slits were used, driven by a spe- 
cially designed clepsydra, one across the image of the sun and the 
other immediately in front of the photographic plate. The instru- 
ment was extraordinarily successful from the outset, and for years 
Hale’s publications described the new and remarkable results ob- 
tained with it. Although his three attempts to photograph the 
corona with the spectroheliograph in full sunlight—at Kenwood, 
Pike’s Peak, and Mount Aetna—proved unsuccessful, they form an 
admirable illustration of his characteristic determination to carry 
the application of an instrument to its extreme limit, a quality which 
proved brilliantly successful on so many occasions. Nearly all our 
existing knowledge regarding the calcium and hydrogen flocculi, the 
forms and motions of prominences, details in the chromosphere, and 
the structure of the solar surface around sunspots is due to the spec- 
troheliograph or to some adaptation of its principles. 
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During this period Hale also carried on some extremely valuable 
investigations of the spectra of prominences and sunspots, identify- 
ing the ultraviolet lines of the Balmer series of hydrogen and numer- 
ous lines of other elements in active prominences, measuring dis- 
tortions of the Ha line and studying its behavior relative to H and 
K, and observing with great care the D3 line of helium both as a 
bright line in the spectrum of the chromosphere or prominences and 
occasionally as an absorption line near active sunspots. The spec- 
trum of sunspots already interested him greatly, and it is quite prob- 
able that he would have anticipated some of his later discoveries 
could his spectroscopic apparatus, necessarily attached to a small re- 
fracting telescope, have had sufficient resolution. Theories of sunspots 
and of the relationship of faculae and flocculi to each other and 
to spots attracted his attention, and he was greatly influenced by 
the work of Secchi and Tacchini in Italy and of Faye in France. 
An interesting summary of many of the results and conclusions of 
his work at this time was published in October, 1892, in Astronomy 
and Astrophysics. 

As Hale moved onward to the wider field of work afforded by the 
Yerkes Observatory, he saw opportunities for extending his investi- 
gations with more powerful apparatus. A large spectroheliograph 
was one of the first requisites, and the Rumford spectroheliograph, 
when attached to the 4o-inch refractor, has remained one of 
the most efficient instruments for certain types of solar research 
ever designed. With it dark hydrogen flocculi were discovered 
and the dark calcium markings which proved to be prominences 
seen in projection against the disk. The dispersion of the instru- 
ment was sufficient to make it possible to set the second slit upon 
different portions of the K line, and thus to map the calcium flocculi 
at different levels throughout the solar atmosphere. A grating spec- 
trograph supplemented the spectroheliograph, and with it Hale made 
the difficult visual discovery of the presence of the green carbon 
band in the spectrum of the chromosphere. He became fully aware, 
at this time, of the possibility of observing the chromospheric and 
much of the “‘flash”’ spectrum without an eclipse, and in future years 
continued such work at Mount Wilson. 

A strangely prophetic remark was made by Hale in an article 
published at the Yerkes Observatory in 1902: “The remarkable 
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peculiarities of the spectra of sunspots seem to deserve more atten- 
tion than they have hitherto received from  spectroscopists.”’ 
Widened lines in the spectra of spots had been observed visually by 
Young in 1872 and were studied systematically by the English 
astronomers Maunder, Cortie, and Lockyer. Hale made some at- 
tempts to photograph the spectrum of sunspots while at the Ken- 
wood Observatory, but his success was only partial because of the 
small size of the sun’s image. With the larger image given by the 
40-inch telescope at the Yerkes Observatory much better results 
were obtained, and not only were the widened lines measured but 
also many of the lines in a green band discovered by Maunder in 
1880. Hale realized that for an adequate investigation of the 
spot spectrum greater linear scale in the spectrograph and a still 
larger solar image were needed. Accordingly, he designed for the 
purpose a large spectrograph used in conjunction with a coelostat 
and long-focus concave mirror, thus planning to bring into regular 
astronomical use an instrument which had been relegated almost 
wholly to eclipse expeditions. This instrument finally took form in 
the Snow telescope, with which in later years at Mount Wilson the 
first of the extensive sunspot investigations was begun. 

Although Hale’s primary interest in life lay in the study of the 
sun, he always emphasized the fact that the sun is an individual star, 
so close that it can be studied in detail, but only one object in the 
vast physical universe about us. For any adequate study of the 
major problems of astronomy, therefore, he felt it necessary to have 
both solar and stellar equipment on an extensive scale, supplemented 
by a physical laboratory to which recourse could be had for the 
interpretation of astrophysical phenomena. This principle underlay 
his plans for each of the observatories he established, and so we 
find great stellar telescopes at Williams Bay, Mount Wilson, and 
on the summit of Palomar. That he developed this point of view 
at an early age and at the same time fully appreciated the striking 
value of that somewhat neglected instrument, the reflecting tele- 
scope, is shown by the gift, through his father in 1895, of the 60-inch 
mirror to the Yerkes Observatory. This is now incorporated into 
the 60-inch reflector on Mount Wilson. 

The 4o-inch telescope at the Yerkes Observatory gave Hale his 
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first opportunity to study stellar spectra. The objects selected were 
the stars of Secchi’s Fourth Type (Harvard type N), which hereto- 
fore had been observed only visually. The carbon bands and a very 
few prominent lines had been recognized, but otherwise the spectra 
were essentially unknown. Observational difficulties were great, 
since all the stars were faint and red and red-sensitive photo- 
graphic plates were very much in their infancy. With the active 
assistance of Ellerman and J. A. Parkhurst, however, Hale obtained 
a series of photographs remarkably good for this period of instru- 
mental development, and discussed the results in detail. The lines 
of numerous elements were identified in the spectra, bright lines 
were observed in several stars, a very significant resemblance to 
certain features of sunspot spectra was discovered, and the physical 
condition of the stars and their place in the sequence of stellar de- 
velopment were considered from numerous points of view. With but 
one or two exceptions the conclusions drawn in this important in- 
vestigation remain without change to the present time, later studies 
with more powerful apparatus supplementing, rather than modify- 
ing, them in any essential way. 

Hale’s great interest in laboratory investigations and his constant 
desire to apply physical tests to astrophysical phenomena led him, 
while at the Yerkes Observatory, to undertake, in collaboration 
with N. A. Kent, an extensive research on the spectrum of the iron 
spark in liquids and compressed gases. At first he thought the dis- 
covery of lines consisting of both emission and absorption com- 
ponents might have application to the phenomena of novae, but 
later he gave up this view. As a physical investigation, however, 
it was important in its bearing on the separation of lines according 
to degree of excitation, the displacement of lines under pressure, 
and the study of selective absorption as a function of wave length. 

The establishment of the Mount Wilson Observatory gave Hale, 
for the first time, the opportunity of building instruments on a 
large scale to fit the problems to be investigated rather than of 
finding the problems to be undertaken with an existing instrument. 
This is the point of view of the physicist and was a principle main- 
tained by Hale throughout his life. The Snow telescope was de- 
signed for certain specific purposes, primarily to make possible the 
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use of larger spectrographs and the violet region of the spectrum; 
the 60-foot tower telescope, to improve seeing conditions on the 
sun’s image and afford convenience of operation and temperature 
control for very long spectrographs; and the 150-foot tower tele- 
scope, to afford a very large solar image for the study of sunspots 
and details on the sun’s surface. The same attitude was taken with 
reference to the stellar instruments, the 60-inch and 1oo-inch tele- 
scopes being designed for use at three foci, so that the magnification 
and the auxiliary instruments could be adjusted to the problem in 
mind. Hale was always greatly interested in the coudé form of the 
reflecting telescope, appreciating from experience the value of large- 
scale stellar spectra and the immense advantages of fixed instru- 
ments operated under laboratory conditions. 

The earliest of Hale’s publications from Mount Wilson are largely 
descriptive, dealing with observing conditions, a subject in which 
he was intensely interested and to which he devoted many experi- 
ments, and with plans for buildings, design of instruments, and a 
program of solar research. The first major investigation was on the 
spectrum of sunspots, and was carried on in collaboration with 
Adams and later with Adams and H. G. Gale in the physical labora- 
tory. The solar equipment used was the Snow telescope and a plane- 
grating spectrograph with a focal length of 18 feet. The photographs 
were greatly superior to any obtained previously, and showed a 
wealth of detail, with many solar lines strengthened or weakened, 
and a great number of lines either not present or very faint in the 
solar spectrum. These latter lines were afterward identified as due 
to bands, those of magnesium hydride by Fowler, and those of 
titanium oxide and calcium hydride by the Mount Wilson observers. 
One result of especial interest was that in the spectrum of many ele- 
ments (such, for example, as iron) some lines were greatly strength- 
ened while others were but slightly affected, if at all. The inquiring 
mind of Hale determined to learn the cause of this behavior, and 
for this purpose had recourse to the physical laboratory. 

It is difficult at the present time, long after the temperature 
classification of spectral lines, the development of the ionization 
theory, and the quantitative analysis of spectra according to energy- 
levels, to realize the situation that existed when these observations 
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of sunspot spectra were made in 1905. Owing mainly to the re- 
searches of Lockyer, the distinction between ‘‘enhanced”’ lines, those 
which first appear or are greatly strengthened in the spectrum of the 
spark as compared with the arc, and the normal arc lines was known 
and its importance recognized; but very little distinction could be 
made among the arc lines themselves. Experiments had shown that 
\ 4227 of calcium could be obtained in the spectrum of a flame, and 
it was sometimes called a “low-temperature line’; but no appre- 
ciable amount of work had been carried on in this field. 

Hale started with the simple and altogether natural working hy- 
pothesis that the differences of intensity between the same lines in 
spots and on the disk of the sun are due to the lower temperature of 
spots. Accordingly, an investigation was begun in the laboratory 
to study the effects of temperature on the spectra of various ele- 
ments. The arc spectrum of iron was photographed, first with a 
rotating arc, then in an arc through which currents of widely varying 
amount could be passed, and finally in the outer flame and in the 
central core of an arc. The results were decisive in showing a great 
difference in behavior among the lines, some being relatively strong 
in the spectrum of the low-current arc and the cooler outer flame, 
while others were strong only in the hot central core and in the 
high-current arc. Comparison with the spectrum of sunspots at once 
showed that the “low-temperature” lines of the laboratory were 
just those which were greatly strengthened in spots while the “high- 
temperature” lines were but little affected. The final evidence that 
the effect was due to temperature was provided by a simple electric 
furnace, which gave results precisely similar to those found with 
the arc. 

This investigation was of great importance not only because its 
results explained many of the most important phenomena of sun- 
spots but also because it had far-reaching consequences. The classi- 
fication of lines according to temperature behavior made at this 
time, and improved and greatly extended by King with the electric 
furnace, laid the foundation for much of the analysis of spectra in 
later years. The discovery of the weakening of the enhanced lines 
in sunspot spectra and the ensuing laboratory investigations led to 
the suggestion that low density might be favorable to the production 
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of enhanced lines, a conclusion so fundamental in the theory of 
ionization. Finally, these results had extremely wide applications to 
stellar spectra, aiding in determinations of temperature and density, 
and through successive stages leading to the spectral differentiation 
of giant and dwarf stars and the discovery of the spectroscopic 
method of determining parallaxes. 

Observations with the spectroheliograph were continued regularly 
by Hale and Ellerman during these years at Mount Wilson, and 
some experimental work with this instrument led to results which 
culminated in Hale’s most brilliant discovery. Sufficient progress had 
been made in sensitizing photographic plates to red light by means 
of dyes, especially by R. J. Wallace, to permit the use of the spectro- 
heliograph with the a line of hydrogen. Photographs with this line 
at once showed a variety of detail, both in dark and bright flocculi, 
which had not been seen previously in observations with the Hé 
line. Most important of all, Hale’s skilful examination at once de- 
tected in the curved form of the flocculi about sunspots evidence of 
vortical motion. Further observations fully confirmed this opinion, 
in one case a very long dark flocculus, which showed gradual curva- 
ture as it approached a double spot, finally forking and being drawn 
into the two centers of the vortex. That he foresaw at once the 
consequences of this discovery is shown by two quotations from his 
article on solar vortices: 

In view of the fact that thé distribution of the hydrogen flocculi frequently 
resembles that of iron filings in a magnetic field, it is interesting to recall the 
exact correspondence between the analytical relations developed in the theory 
of vortices and in the theory of electromagnetism. 

Double lines, which look like reversals, have recently been photographed in 
spot spectra with the 30-foot spectrograph of the tower telescope, confirming 
the visual observations of Young and Mitchell. It should be determined 
whether the components of these double lines are circularly polarized in opposite 
directions, or, if not, whether other less obvious indications of a magnetic 
field are present. I shall attempt the necessary observations as soon as a suitable 


spot appears on the sun." 
On June 25, 1908, Hale obtained a series of photographs with the 
30-foot spectrograph, using a Fresnel rhomb and Nicol prism. These 


* Contributions from the Mount Wilson Observatory, No. 26; Astrophysical Journal, 
28, 311, 1908. 
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gave unmistakable evidence of the Zeeman effect and of the presence 
of a magnetic field in sunspots. 

Following this remarkable discovery, Hale devoted considerable 
time to measurements of the strength of field in spots, comparisons 
with laboratory results, studies of plane polarization across the lines 
of force, and experimental work on vortex models. Since the resolu- 
tion of the components of many of the spot lines was beyond the 
power of the 30-foot spectrograph, he decided to postpone further 
extensive investigations until the completion of the 150-foot tower 
telescope and the 75-foot spectrograph, which were under construc- 
tion in 1909 and 1g1o. The decreasing sunspot activity, however, and 
the scarcity of spots led him at this time to undertake an investiga- 
tion particularly well suited for sunspot minimum but extraordi- 
narily difficult and exacting because of the smallness of the quanti- 
ties involved. This was the problem of the existence of a general 
magnetic field of the sun as shown by the Zeeman effect. Here there 
could be no question of the separation of lines into components, 
but, as computation showed, only of minute displacements, when 
the Nicol prism was rotated, of three or four thousandths of a milli- 
meter even on the great scale of the spectra with the 75-foot spectro- 
graph. A large amount of observational material was required, and 
the measurements were undertaken by van Maanen and other mem- 
bers of the Observatory staff. Positive results were obtained which 
were discussed by Hale in an extensive article written in March, 
1913; but the problem was rendered even more difficult by the ap- 
parent failure of certain lines to show the expected effect. The in- 
vestigation was continued by Hale up to the end of his life, and he 
developed several ingenious and effective devices for measuring the 
minute quantities involved. New series of photographs were ob- 
tained near the sunspot minima of 1922 and 1932; and the spectra 
were measured, but with somewhat inconsistent results. It is very 
difficult to draw a definite conclusion regarding the outcome of this 
long research undertaken with such great skill and patience. The 
presumption for the existence of the general magnetic field seems to 
be strong, but the definite proof may have to await new methods or 


improved instruments. 
In the autumn of 1918, Hale published, in collaboration with 
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Ellerman, Nicholson, and Joy, one of the most important of his 
papers, entitled ‘“The Magnetic Polarity of Sunspots.’”’ The discov- 
ery of the reversal of polarity of sunspots with the spot cycle, fully 
confirmed by later results, should be ranked with that of the Zeeman 
effect in its importance to solar theory. In this article are also found 
the extremely interesting conclusions that the two members of a 
bipolar spot group have opposite polarities, and that the polarities 
of the preceding and following spots in such groups are opposite in 
the northern and southern solar hemispheres. Seven years later, in 
1925, Hale and Nicholson, with more extensive material, discussed 
the law of sunspot polarity as based on 1735 groups observed be- 
tween 1908 and 1925, a period covering two maxima and minima. 
The reversal of polarity with the spot cycle was fully established, 
with only 41 exceptions in the entire number observed. The full 
sunspot period, corresponding to the interval between the successive 
appearances in high latitudes of spots of the same magnetic polarity, 
is, therefore, twenty-three years. The reversal of polarity was found 
to be a sudden phenomenon unheralded by any increase of spots of 
irregular polarity or other exceptional event. It is fitting that a 
series of communications containing such brilliant discoveries in 
solar physics should have come to a close with one of such funda- 
mental and far-reaching importance. 

The invention of the spectrohelioscope and its numerous applica- 
tions are described by Hale in a series of four publications between 
1929 and 1931. “he fascination which changing solar phenomena 
had constantly held for him from his boyhood at Kenwood was 
never stronger than during this last period of his active observing. 
In these articles he gave many results of his own observations, sum- 
marized and discussed progress in studies of the sun, and outlined 
a program for obtaining continuous records of solar activity through 
a series of stations distributed suitably in longitude throughout the 
world. These plans he lived to see largely accomplished. Two of the 
most promising fields for solar study, the relationship of solar erup- 
tions to terrestrial phenomena and the recording of rapidly changing 
prominences and eruptions through continuous high-speed photog- 
raphy, Hale fully recognized and frequently discussed, and had 
health permitted, he would doubtless have developed them exten- 
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sively with the aid of all the powerful resources of modern instru- 
ments and technique. 

Remarkable as were Hale’s personal discoveries in science, his 
contributions in other ways were no less great. Apart from the ac- 
complishments of the institutions he organized and the co-operative 
investigations he initiated, his influence through the opportunities 
he afforded to other scientists to carry on special investigations was 
of immense value. He brought E. F. Nichols to the Yerkes Ob- 
servatory to make the first measurement of the heat of a star, and 
Michelson to Mount Wilson to measure successfully a star’s diam- 
eter and improve so greatly our knowledge of the velocity of light. 
He was in constant touch with scientists all over the world, corre- 
sponding with them, studying their problems and methods, and often 
offering facilities for the extension of the work upon which they were 
engaged. For this reason Kapteyn came to Mount Wilson soon after 
the completion of the 60-inch telescope, to aid in the formulation of 
a program of stellar research; Julius, to investigate the effects of 
anomalous dispersion where observational material was immediately 
available; and Abbot, to undertake measurements of the solar 
radiation. 

The great contributions made by Hale to science were recognized 
by learned societies all over the world, and the record of the honors 
and distinctions conferred upon him is a long and distinguished 
one. He was a member both of the leading academies and the sci- 
entific societies of England, France, Italy, Holland, Sweden, Nor- 
way, Belgium, and several other countries, in addition to those of 
the United States; was elected president of the International Council 
of Scientific Societies in 1931; and served as a member of the Com- 
mittee on Intellectual Co-operation of the League of Nations in 
1922. His foreign membership in the Royal Society of London dated 
from 1909, and his associateship in the Institut de France (Académie 
des Sciences) from 1919. Practically every distinction in the form 
of medals which scientific societies at home and abroad could award 
for contributions to astrophysics was bestowed upon Hale, his ex- 
tensive list of honors culminating in the award of the Copley Medal 
of the Royal Society in 1932. As a trustee of the California Insti- 
tute of Technology and of the Huntington Library and Art Gallery, 
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he also exercised a powerful influence upon the intellectual and cul- 
tural life of the city in which he lived; and an interesting tribute 
to his accomplishments in this direction was the award by the city 
of Pasadena of the Noble Medal for civic service. 

To his own colleagues Hale was an unfailing source of help and 
encouragement. His insight and resourcefulness, his enthusiasm and 
joy in his own work and that of others, his hopefulness, modesty, 
and great personal charm, made any association with him a constant 
delight. So strong was the influence of his personality that even in 
the later years of his illness, when, of necessity, he had to withdraw 
from frequent contact with his associates, the sense of his presence 
and knowledge of his continuing thought of them and their work 
remained as a permanent influence. Finally, in its breadth of inter- 
ests in the problems not only of physical science but of those of the 
imagination and human affairs, Hale’s life was a most remarkable 
one. It is perhaps symbolic of this man of great gifts and wide 
horizons that he who had devoted his life to the study of the nearest 
star should find his last deep interest in an instrument destined to 
reveal the remotest objects of our physical universe. 


CARNEGIE INSTITUTION OF WASHINGTON 
Mount WILSON OBSERVATORY 
March 1938 
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A COMPARISON OF SPECTROSCOPIC AND 
TRIGONOMETRIC PARALLAXES* 


HENRY NORRIS RUSSELL’ AND CHARLOTTE E. MOORE 


ABSTRACT 


The calibration of spectroscopic absolute magnitudes has been made by selecting 
groups of stars with similar intensities for pairs of ‘‘sensitive’’ lines and adjusting the 
calibration-curves so that the absolute magnitudes which these lines give for each 
group shall agree with those determined for the same stars from other data (trigo- 
nometric parallaxes, proper motions, etc.). Since the spectroscopic data are affected 
by accidental errors, this practically unavoidable process gives a regression-curve 
which inevitably underestimates the deviations of the individual absolute magnitudes 
from the general mean for the spectral class concerned, while giving an accurate value 
for this mean. 

The magnitude of this effect can be determined by the comparison of the spectro- 
scopic and trigonometric parallaxes of the same stars, provided that the accidental 
errors of the latter are known. 

Formulae are developed for the determination of the constants involved in this problem, 
and of the statistical errors of the constants derived from samples of finite size. 

The spectroscopic parallaxes of Mount Wilson Contribution No. 511 are then com- 
pared with the trigonometric parallaxes given in Schlesinger’s General Catalogue (Yale, 
1936). For the 1140 stars of the main sequence which are available, the real dispersion 
in absolute magnitude about the mean for the subclass of the Draper classification 
corresponds to a standard deviation of + 34 per cent in the parallax, while the errors 
of the spectroscopic determination (including the effects of real differences between 
stars whose sensitive lines appear the same) give a standard deviation of + 38 per cent. 
The means for all the stars of a given subclass are accurate within 5 per cent, except fora 
few faint M dwarfs. The differences between the individual spectroscopic absolute 
magnitudes and the mean for the subclass should be multiplied by 1.19 to give an im- 
proved calibration for groups selected according to the spectroscopic criteria, but by 
2.6 in order to give values which will be correct in the mean for groups of stars selected 
by criteria independent of the spectroscopic data. When these stars are divided into 
seven groups according to spectral type, the results for all are similar, although with some- 
what larger fluctuations than would be anticipated from errors of sampling. Comparison 
of the mean for groups of stars selected by spectroscopic absolute magnitude and by 
proper motion fully confirms the anticipated differences. 

For the giant stars, 732 in number, the parallaxes are smaller; but good mean results 
can be obtained. The standard deviation in parallax, corresponding to real dispersion 
in absolute magnitude, is +52 per cent (increased by the inclusion of supergiants and 
subgiants), while that arising from the spectroscopic errors is + 38 per cent. The general 
means of the spectroscopic and trigonometric parallaxes agree perfectly. The differences 
from the mean require no correction for grouping according to spectroscopic criteria, 
but a factor of 1.6 when the criteria are independent. 

A special discussion of 57 stars of great absolute brightness shows that the probable 
errors assigned in “‘Schlesinger’s Catalogue’ represent correctly the true accidental errors 
of the determinations. Formulae and tables are given (sec. 13) for the rapid application 
of the corrections to the spectroscopic parallaxes. These are important for only a small 
fraction of the whole number of stars. Their applicability to supergiants is doubtful. 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- 
ington, No. 589. 
' Research Associate, Mount Wilson Observatory. 
389 


} 

| 

| 

| 


390 HENRY NORRIS RUSSELL AND CHARLOTTE E. MOORE 


I. THE CALIBRATION OF SPECTROSCOPIC ABSOLUTE MAGNITUDES 


The validity of the spectroscopic method of determining absolute 
magnitudes and of the resulting parallaxes is fully established; but 
the calibration of the spectroscopic data presents certain complica- 
tions. Suppose that we have, for a considerable number of stars, 
estimates or measures of the spectroscopic criteria associated with 
absolute magnitude, and also determinations of the ‘‘observed’’ 
absolute magnitudes, based on other evidence. An obvious and di- 
rect method of calibration is to divide the stars into groups, within 
each of which the spectroscopic criteria are the same or very similar; 
take the mean of the observed absolute magnitudes for each group; 
and plot these means against the numbers representing the mean of 
the spectroscopic criteria. When calibration-curves of this sort have 
been drawn separately for stars of different sorts, e.g., giants and 
dwarfs, each such class of stars presents, for our present purpose, an 
independent problem. In discussing this, we may assume that we 
are dealing with stars within a narrow range of spectral type as 
ordinarily defined (e.g., dG4) and that their number is sufficient to 
justify a statistical discussion. For main-sequence stars thus se- 
lected the dispersion of the true absolute magnitudes about their 
mean is usually not great, the standard deviation being less than 
1 mag. For a group containing ordinary giants, supergiants, and 
subgiants it is greater. We will here ignore complications which 
may arise if the criteria defived from different pairs of lines are not 
fully consistent, and any systematic errors which may be present in 
the observed absolute magnitudes, and discuss the influence of the 
accidental errors. Both the spectroscopic and “‘observed”’ data will 
contain such errors; and, since the stars have been grouped accord- 
ing to the values of the former, the calibration will give a regression- 
curve, whose course is influenced by these errors in a well-known 


manner. 
2. THE REGRESSION-CURVES 


Let M be the mean of the true absolute magnitudes, for all the 
stars in the group, and 


M=M+w (1) 
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the true absolute magnitude of an individual star. For the ‘‘ob- 
served” absolute magnitudes O let O, be the mean systematic and 
x the accidental error. Then 


(2) 


If the calibration-curves were ideally drawn, we would have for the 
spectroscopic absolute magnitude S 


S=M+wty, 


where y is the accidental error. The curve actually employed may, 
however, require correction in zero point and slope. The former 
leads to a constant correction S,. The latter may be allowed for by 
multiplying the spectroscopically determined difference w + y by a 
factor /, so that 

S=M+S,+l(wty). (3) 


A complete description of changes in the calibration-curve would in- 
clude changes in its shape also; but these cannot be expressed by 
simple equations, and the summary process of correction, which 
amounts to changing only the zero point and horizontal scale of the 
curve, will be adopted here. 

We must now make assumptions regarding the quantities w, x, 
and y. To assume that their mean values are zero amounts only to 
defining M, O,, and S,, so that in the mean of a very large number of 
cases 
O=M+0,, S=M+S,. 


Then 


O-O=w+x, S-S=lwt+y). (4) 


We will next assume that w?, x?, and y? have definite mean values 
W?, X?, and Y?, the same for all large random samples of the group 
(with the usual statistical tolerance); that is, they are drawn from 
statistically homogeneous ‘‘populations.”’ 

Finally, we will assume that the mean values of the products xw, 
xy, and wy tend toward zero for a large number of entries. For «w 
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and xy this may be accepted, since by hypothesis x, the error of O, de- 
pends on observations of some kind quite independent from the rest. 
It is not necessary to assume that the distribution of x is normal, or 
symmetrical about the mean; only that there is no tendency for its 
average value to be influenced by the observational errors of S, or 
by the fact that the star observed is absolutely brighter or fainter 
than the average. 

Correlation between w and y would mean that for stars of great 
and small absolute brightness the effect of absolute magnitude upon 
the spectroscopic criteria was different. This is entirely possible; 
but the mean effect (with which we are here concerned) will be auto- 
matically allowed for in drawing a good calibration-curve, and errors 
in this process have already been taken into account (approximately) 
in the factor /. The assumption that wy = o is therefore a legitimate 
approximation. 

If then we set 


u=S-S=Iwty, w=O0O-O=w+x, (5) 
we have the mean values 
We PW? + PY? | 


— W? + X?, (6) 
LW? . 


| 


We may now define the regression line representing the mean 
values of O for specified values of S, by solving by least squares the 
equations of condition 


where a is the constant to be determined. The solution (in the Gauss- 
ian notation) is 
a[u?] = [uu’] ; (7) 


whence, by (6), the regression line is 


Ww? 


(8) 


u = au 
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This process ignores the possible curvature of the regression-curve; 
but, since we are supposed to start with an approximately calibrated 
set of spectroscopic absolute magnitudes, this should not be seri- 


ous. 
To secure the best agreement between the mean observed and 


spectroscopic absolute magnitudes, for stars grouped according to the 
values of the latter, we must multiply the provisional values of S — S, 
with which we started by the constant a defined by (7). With these 
newly calibrated values, the regression line will be « = wu’ and the 
new value of / will be 


W2 
L = wey: (9) 
This method of calibration involves the well-known tendency to se- 
lect positive errors when S is large, and negative errors when it is 


small. 
For groups selected by some impartial method, such that there is 


no preference for positive or negative errors in either S or O, the re- 
lation is different. To obtain significant results we must have a cri- 
terion which separates—partially, at least——-stars of different real 
brightness. For such a group the means of x and y will tend toward 
zero, but that of w toward some value w,, so that 


tin = lw,; = Wy. 
The mean points will lie on the line 
= (10) 
Finally, for groups selected according to the value of O, we find, for 
the equation of the regression line, 


“ = u (11) 
lW? 


There is nothing new in an analysis of this type.? It should, how- 
ever, be noted that it has not been assumed at any point that the 


2 Cf. Seares, M.N., 84, 15-20, 1923. 
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distribution of the quantities w, «, and y follows the normal error law. 
The assumptions actually made are much less restrictive. If, how- 
ever, the distributions of w, x, and y are normal, it is easy to show 
that the distribution of « and w’ follows the normal law in two vari- 
ates, so that the regression-curves are actually the straight lines (8) 
and (11). 

The equations (6) do not suffice to determine the unknowns. But, 
if the dispersion X can be calculated from a study of the “‘observed”’ 
data, W, Y, and / may be found. It should be noted, however, that 
the regression line (8) involves only [wu] and [ww’] and is independent 
of the value of X (as it obviously must be, since, when taking means 
of O for a fixed value of S, the errors of O tend to vanish in the mean 
of a large number of cases). 


3. APPARENT AND REAL DISPERSIONS 


If we should adjust our calibration so that the slope of the regres- 
sion line (8) was 1, that of the ‘‘impartial’”’ line (10) would be 


I+ 


This slope will be steeper, the greater is the ratio of the mean acci- 
dental error of a spectroscopic determination of absolute magnitude 
to the true standard deviation of the absolute magnitudes about their 


mean. 
The standard deviation of the spectroscopic absolute magnitudes, 


as thus calibrated, would be 


Ws 


The true standard deviation, including the effects of accidental 
errors, is 


Hence, with this adjustment, 


‘ad oW? 
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If the calibration is not thus adjusted, 
o, = + Y2)! = Ie. 


The method of calibration leads therefore to an underestimate of 
the errors of the spectroscopic parallaxes and to the real dispersion 
of the absolute magnitudes. This underestimate may be serious, but 
it is in no sense an error of calculation; it is the inevitable result of 
the method of calibration, and its amount is not diminished by in- 
creasing the number of stars, so long as the accidental errors and 
true dispersion remain the same. This type of calibration appears, 
nevertheless, to be the best that is generally available; for the cor- 
rection from the regression to the ‘‘impartial”’ line can be made only 
when the accidental errors of the individual “observed” values can 
be accurately estimated, which is not always possible, especially 
when mean parallaxes, etc., have been found from proper motions. 
It should be emphasized that the underestimate affects only the 
deviations from the mean for the group (in both directions). The 
group-means themselves are undisturbed. 4 


4. CORRESPONDING RELATIONS FOR THE PARALLAX 


The practical determination of these quantities is complicated by 
the fact that there is no direct method for determining absolute 
magnitude. Some of the best methods for getting mean values, such 
as those depending on proper motions, provide no way for finding 
the errors w’ for individual stars. The mean value of u? may be well 
determinable in other ways, but those of the products wu’ cannot be 
found at all. 

For an investigation of the correlations and regressions, the best 
material is offered by trigonometric parallaxes. It is useless to com- 
pute absolute magnitudes directly from these, since some of them 
are negative. We must transform the spectroscopic absolute mag- 
nitudes into quantities of the dimensions of a parallax. 

The actual spectroscopic parallaxes, s, are inconvenient for this 
purpose, since they are affected by the very large dispersion in the 
distances of the stars. We may obtain convenient quantities by com- 
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puting the values s’ which the spectroscopic parallaxes would have if 
the stars all had an apparent magnitude equal to S,. Then: 
S=S,+5+ 5slogs’, 
and, by (5), 
u=S,—-S+5+ 5logs’ ; 


whence 


10s’ = exp 0.46 (u + S — S,) 
[exp 0.46(S — S,)] (1 + 0.46u + 0.1060? 


+ + 0.002u4....). 


The mean value of s’ is given by 
10s’ = [exp 0.46 (S — S,)](1 + 0.106?) , (12) 


the higher terms being negligible. 


If we set 
s=asto, 
we have 
v u + 0.2 
= = 0.46 3 ( w) ; (13) 
I + 0. 1064? 


In practice S, has been taken nearly equal to S, so that s’ is close 
to o" ro. 

If the dispersion in w is not too great, u and v will be nearly pro- 
portional over most of their range. If the distribution of w is normal, 
that of v will be asymmetrical; but, in the absence of a priori knowl- 
edge of the distribution of absolute magnitudes, consideration of this 
would be a refinement appropriate only to a second approximation. 
If mis the apparent magnitude and we set 


s' =fs, 
we have (14) 
log f = 0.2(m — S,). 


3 In this paper “‘log”’ implies the base ro; ‘‘/m,”’ the base e. 
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For comparison with the ‘‘reduced”’ spectroscopic parallaxes s’, 
we must multiply the trigonometric parallaxes ¢ by the factors f, ob- 
taining 

= ft. 
If p is the true parallax, we may set 


t= Cp(1 + 2), (15) 


where C is a constant factor representing the mean effect of system- 
atic error and z represents the accidental error of observation. 
If M is the true absolute magnitude, 


m—-M=-—5-— 5logp; 
whence, by (1) and (14), 
log f = —1 — logp+o.2(M—S,+w). 


Hence, 


C (1 + 2) expo.46 (M — S, + w) 
C [exp 0.46 (M — S,)](1 + 2 + 0.46w + 0. 106w? 
+0.462w+....). 


The mean value is 


= C [expo.46 (M — S,)|(1 + 0.106w? + ....). 
Setting 
we have 
2 +0.46w + 0.106 (w? — w').... (16) 
1 + 0. 106w* 


If x is the change in absolute magnitude corresponding to the 
change z in parallax, 


1+ 2 = expo. 46x , 
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or 


z= 0.46x + 0.106%7.... 


Hence, with sufficient approximation for our purpose, we have 


0.46 (w+x)+...., | 
(17) 


| 
) 


In an exhaustive investigation, which sought to determine the 
form of the distribution functions of the errors, and of the regression- 
curves, it might be necessary to take account of the terms neglected 
in (17); but the present study represents a first, and not a final, ap- 
proximation. 

If 


= 0.46w, W, = 0.46W , etc., 


we now have in place of (6) 


a 2 2 2 

= + Fi), 

(18) 
vv 

= =/W?. 


The quantities v and v’ may be found from the observed spectro- 
scopic and trigonometric parallaxes. Moreover, the probable errors 
of the latter, if reliable, determine X?7, so that the other unknowns 
may be evaluated. 
5. WEIGHTS 

In applying these principles we must take account of the great 
differences in the percentage errors of the trigonometric parallaxes. 
If ry is the probable error assigned to #, that of the “‘reduced”’ parallax 
t’ will be r’ = fr. The weight p, of ¢’ varies inversely as r’”?. 


|| 
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In a discussion of the masses of the stars, already in progress, the 
unit of weight has been taken to correspond to a probable error of 
+28 per cent of the mean parallax, and this convention has been 
adopted here. The mean value ?’ is usually close to 0” 10, so that we 


may adopt 
= (19) 


The percentage error of a spectroscopic parallax is approximately 
constant. Its value was provisionally assumed to be + 20 per cent, 


TABLE 1 
Exact Assigned Exact Assigned 
Weight Weight Weight Weight 
2.00 0.23 
1.8 0.20 
1.4 14 
1.0 10 
.07 
0.5 05 
0.35 035 
0485 0.29 030 
0.25 0.025 


which gives p, = 2 in all cases. Only a large error in this assumption 
would change the weights enough to require revision. 

To assign different weights to the same star when discussing s’ 
and ¢’ has obvious disadvantages. Both quantities, for any star, 
have therefore been given a weight /, defined by the usual equa- 
tion 


I z I 
=—+-, (20) 
To save labor in the computation, uniform weights were assigned 


to all stars for which r’ lay between certain limits, as in Table 1. A 
coarser subdivision would have been practically satisfactory. 


| | 
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6. FINAL EQUATIONS 
We now have 


N 


N- I 


= [pwi] + + [pri] . 


The factor V/(N — 1) is introduced to take account of the fact that 
the number JN of entries is moderate, instead of very large, as pre- 
viously supposed. 

The mean value of w7 is always W?, while, for an entry of weight 
pi, that of x7 is X7/p., where X, is the standard deviation for unit 
weight. 

Hence (since p/p, = 1-p/p, and Pp, is constant), 


N [po’v’] = 72 p 72 | 72 


To p, = 1 corresponds a probable error r’ = +0.028, and a mean 
error X, = +0.0416 for v’. Hence 


X, = (22) 


and we have 


If we set 
N I 
Br? 
iar (23) 
we have 
+ (plex: = (24) 


N-1 ff 


| 


| | 
alll 
| 
| 
| 
| 
| 
| 
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Similarly, we find 


[pllW? N I NF (25) 
+ = = NE, (26) 
= 2 


where Y, = Y,/ Vp, and is the mean accidental error of a spectro- 
scopic determination of parallax. 

Since p, = 2, and X, = 0.0416, these equations determine W,, 
i, and. 

When several groups are to be combined with one solution, these 
equations may be added as they stand, introducing the proper nu- 
merical values of B and X, into each. If ¢’ is the same for all the 
groups, equations (22) and (23) may be used with the combined 
values of NV and [/]; otherwise this process is slightly inaccurate. 

The group-means s’ and ¢’ will themselves be affected by acciden- 
tal errors. Better values, sof., may be obtained by plotting s’ and 
i’ against the spectral class and drawing smooth curves. If 


As=s,—s, At 


Il 
| 
~ 


we have, referred to the new standard: 


[pvr]. = [pov] + [p]As?, | 
= + [plar, | 
[pvv'], = [pov] + [p]Asdt . 


To retain the factors V/(N — 1) now slightly exaggerates the correc- 
tion for the finite size of the groups; but this is practically harmless. 


7. STATISTICAL ERRORS OF SAMPLING 
With groups of moderate size the quantities | pv], etc., will be sub- 
ject to considerable fluctuations, owing to random sampling—of the 
order of NV’? times their mean values. 
If we set 


e=l(wty), (28) 


| 
| 
| 
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we have, by (24), on the average, 


N [peal 


[pee] = a NE or pee= E, 


where the bar denotes the mean value. 
For an individual group of N data, this will not be exactly true; 


and if we set 


= NP, - yo, (29) 


we may write 


pee = E+ gor Ng = [pee] — NE = [pee] 
~ + 
and similarly, from (21) and (26), 


= +e’, Neg = [pee] — NE 
= [pe'e’'] — N(PW? + OX?), 


pe’ =F+h, Nh=([pee’] — NF =[pee’] — NPIW?. 


(30) 


For the mean-square values we have 
N2g'? = [pe’e’}? — 2NE[pe’e’] + N?E”? = [pe’e’? — N°E” , andsoon; 
also for the product 


N?gh = [peel[pee’] — N?EF,, etc. 


We have here to take the means of a large number of squared 
sums like [pee] with different values of w, x, and y, but the same 
weights p throughout. We have 


[pe’e’? = {[px?] + 2[pxw] + [pw’]}?. 


The mean values of all sums containing odd powers of x or w vanish, 
so that 


[pe’e’ = + [pw]? + + 


| 
| 
| 
| 
| 
| 
| 
| ee | 
| 
| 
| 
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Now 
= + .)? = pixt t+ + 


For any entry (say x.) there is a trigonometric weight /,,, fixed for 
all the successive trials, such that 


Pr, 
For «4 we may assume 
xi = (1 +c) (31) 


For a normal distribution, c = 2. We will at least assume it constant 
for all entries. We then find 


[pep = Xtc xi]? 
Pi pr 


Now p/p: = 1 — p/ps, and p, is constant. Hence, 


Set 
p = (1 + A)(p)? = (1 + ADP? (33) 
and we have 
(34) 


where 


| 
| 
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We find also 
|= - 40. 
pi 
We now have 


= NcX3(Q? + AQ”) + 
Similarly, but more easily, we find 

= + A)P?V2 + 

(pw? = + + , 

= NP(Q — AQ’)XiY;, 

[pxw? = NP(Q — AQ’)XiW;, 

[pwy)? = NP2(1 + A)W2¥?2, 


We now find, easily (setting c = 2), 


N g? = 2{E? + A(E’ — X?)}, 

Nit = EE’ + + — X2)}, 
N gg’ = 2(1 + ADF’, 

N gh = 2(1 + A)EF, 

N gh = 2E'F’ + 2AF(E’ — X?). 


NOX? : [py] = NPY? : [pw] = NPW?. 


(35) 


(36) 


> (37) 


The statistical mean errors of E, E’, and F are, respectively, g, g’, 
and Under the assumed conditions QN7 (=QN3/03) is constant, 
since Q depends on the weights, which are fixed, as is ¢,, while X, 
has been accurately determined from a great number of observa- 
tions, largely extraneous to the present discussion. The mean error 


of the calculated PW; is therefore g’. 


We are especially interested in the errors of the coefficients, a and 
1/l, which define the calibrating lines (8) and (10). We have 


pe’ F+h 
pe 


| 
| | 
| 
| | 
- 
) 
| 
| 
| 


COMPARISON OF PARALLAXES 405 


If a, is the mean value F/E, we have, to the first order, 


whence, by (37), 
NE*(a — a)? = EE! — F? + A(EE’ — EX? — (38) 
Similarly, since 


we find, after some reckoning, 


+ + 40"). (39) 


Apart from the last term in (39) the ratio of the mean errors is 


E(E'— OX?) _ W?+¥2 


which is the ratio of the slopes of the lines. The additional term 
shows that the “impartial” constant 1// is always more poorly de- 
termined than the regression constant a. 


8. THE OBSERVATIONS 


Significant results in a statistical inquiry of this sort can be ob- 
tained only if extensive and homogeneous observational data are 
available. Schlesinger’s General Catalogue of Stellar Parallaxes* pro- 
vides all that could be desired for the trigonometric parallaxes. The 
systematic errors for the various observatories have been very care- 
fully discussed and applied, and the absolute value of the correc- 
tions checked by all available criteria. Schlesinger remarks (p. iv): 
“T regard it as more probable than not that the correction needed 
to make this system absolutely true as a whole does not exceed 
“oo1.”’ The assigned probable errors have also been very carefully 


4 Yale University Observatory, 1935. 
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adjusted, removing the fluctuations due to random sampling, and 
allowing for the nonsystematic discordances between the results of 
different observatories. New evidence will be presented below that 
they correspond closely to the true accidental errors. The trigono- 
metric parallaxes of this Catalogue have been adopted without mod- 
ification in all cases. 

The spectroscopic parallaxes here discussed are those of the latest 
Mount Wilson list. While all these have been derived by the same 
method, many different calibration-curves have been employed, de- 
pending on various pairs of lines. 

We are indebted to Dr. Adams for full details of these curves. 
One set was used for stars belonging to the main sequence and a 
quite independent set for giants, including supergiants and sub- 
giants. The range of spectral class over which a given pair is sensi- 
tive to luminosity varies widely; and, within this range, several 
roughly parallel curves were usually employed for different spectral 
subclasses. The exact spectral criteria employed therefore vary from 
subclass to subclass (and to a smaller degree from star to star); but 
over a moderate spectral range this change is unimportant. 

Few of the calibration-curves are rectilinear, but their curvature 
is smooth, and there are no examples of sharp bends of such a char- 
acter (e.g., in the interval between giants and supergiants) as to 
make it probable that the main part of a necessary revision could 
not be made by changing’ the zero point and slope, as described in 
section I. 

The conditions laid down in section 1 appear therefore to be sat- 
isfied, so long as giants and main-sequence stars are kept separate 
and all giants are discussed together. A few exceptional objects (such 
as the abnormally faint A stars called ‘“‘white dwarfs” in the notes), 
which were calibrated by special methods, have been excluded from 
the discussion. Every star calibrated by the ordinary methods, how- 
ever discordant, has been included. 

The published values of the spectroscopic absolute magnitudes 
have in all cases been adopted. The apparent visual magnitudes 
given in Mount Wilson Contribution No. 511 have been altered in a 


5 Adams, Joy, Humason, and Brayton, ‘The Spectroscopic Absolute Magnitudes 
and Parallaxes of 4179 Stars,” Mt. W. Contr., No. 511; Ap. J., 81, 187, 1935. 
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few cases—with corresponding changes in the spectroscopic paral- 
laxes. Most of these are for double stars of wide separation for which 
the magnitudes given in the Henry Draper Catalogue refer to the 
combined light, but had been assumed, by the compilers of the list, 
to belong to the brighter component. A number of such cases were 
cleared up by reference to other Harvard publications.® For two or 
three stars photometric measures made at Harvard, and published 
in the senior author’s book on parallaxes,’ have been adopted in place 


TABLE 2 
STANDARD SPECTROSCOPIC MAGNITUDES, So 
MAIN SEQUENCE GIANTS 

Sp So Sp So Sp So Sp So Sp So 

At | +0.9 F2 +2.8 G3 +4.4 K4 + 6.4 | Go-Gs | +0.7 
A2 Ev F3 3.0 G4 4.5 K5 6.7 | G6-Gg | +0.6 
A3 F4 G5 4.6 7.3 Ko +0.5 
A4 Fs 3.2 G6 4.8 | Mo 8.3 | Ki-K2] +0.4 
AS ee F6 3-4 G7 5.0} Mr 8.8 K3 +0.3 
A6 1.9 F7 G8 M2 9.4 K4 +o.2 
A7 2.0 F8 3.6 Go 5-3 | M3 10.0 K5 +o.1 
A8 2.2 Fo 3.8 Ko &.5 M4 10.6 |Mo-M2]| —o.2 
Ag 23 Go 4.0 Ki 5-7 | Ms | +11.9 M3 —0.3 
Fo 2.5 G1 K2 M4-M7| —o0.4 


of the corrected Bonner Durchmusterung magnitudes given in the 
Henry Draper Catalogue. One misprint in Mount Wilson Contribu- 
tion No. 511 was discovered: for the star 4386A, a 17"12™1, 6 — 34° 
53’, s should be 0” 166 instead of 07116. 


Q. SUMMARY OF THE DATA 


In applying the equations (14) it is not necessary to start with an 
accurate value of the assumed absolute magnitude S,, since any 
error in this is automatically eliminated from v/s’ or v’/t'; but it is 
convenient to do so. The mean values of S for each spectral sub- 
division were plotted, and smooth curves drawn to represent them, 
with the results in Table 2. The means for the giants are derived 
from normal giants alone. 

6 Harvard Ann., 56, 231 (No. 7), 1912; 64, 165 (No. 6), 1912. 

7 Pub. Carnegie Institution of Washington, No. 147, 1911. 
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For each separate spectral subdivision, the values of s’, t’, r’, and 
the weight p were computed for each star; the weighted means s’ 
and ?’ were taken, and the residuals v and v’ from these means tabu- 
lated. 

For the main sequence the values of s’ gave a smooth plot against 
spectral class, while those of ?’ were sometimes ragged. Smooth 
curves were drawn to represent each, and improved mean values 
s, and ¢, read from these. The quantities | zz], etc., were corrected 
by (27) to give the sums [vz], corresponding to these new means. 
From these the quantities [p]B, NE, NE’, and NF defined in equa- 
tions (23) to (26) were calculated. The resulting data are given in 
Table 3. The second column gives the number of stars, the third, 
the sum of their weights. 


10. PRELIMINARY DISCUSSION 


The mean weight increases along the sequence because the later- 
type stars have, on the average, larger trigonometric parallaxes and 
the percentage errors of these are smaller. The values of s’ show a 
narrow range (except for Ar and A2, where the course of the curve 
adopted for S, is uncertain). The weighted mean of all the values is 
©.1032, for 1140 stars. 

The difference between this and the expected value 0” 100 arises 
mainly from the fact that the values of S,, used in calculating s’, 
represent the means of the absolute magnitudes. If S is the individ- 
ual absolute magnitude 


10s’ = expo. 46 (S — S,) ; 
whence 
= 1 + 0.106 (S — S,)?. (40) 


From a tabulation furnished by Mr. Joy, it is found that, for 1707 
main-sequence stars of spectra Ao to K3, (S — S,)? = 0.208, so that 
we might expect a mean value s’ = 0.1022. The remaining dis- 
crepancy of 1 per cent (if real) would indicate that the stars observed 
trigonometrically were, on the average, 0.02 mag. fainter than the 
general run of those observed spectroscopically—a wholly insignifi- 
cant difference. 
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The tendency of s’ to run higher in the region F5-Gs5 than in G6- 
K6 may be real, but could be removed by altering S, by 0.1 mag. at 


TABLE 3 
OBSERVATIONAL DATA: MAIN SEQUENCE 


Sp N R So lg NE NE’ NF 
100 

6 5.3 | 07157| 07116) 0.98 0"130| o7116| 2.4 | 0.92 2.66] 0.79 
| .109| 0.98 -122| .111| 4.8 | 0.62 0.81|—0.03 
AZ... BO 9.3 .111| 0.98 .114| .106] 10.1 | 0.53 1.66) 0.33 
Aa...) 22] 12-8 -109] .087| 0.98 .102} 15.5 | 0-57 
5. .096| .083| 0.97 -098| 12.6 | 0.32 2.98} 0.63 
A6...| 15 | 10.8 .086| .100| 0.97 .097| 10.2 | 0.16 1.46} 0.14 
A. .095| 0.96 .092| .095| 8.3 | 0.40 3.38| 0.48 
A8.. II 8.5 .096} .099| 0.90 .095| .094| 7-7 | 0-20 3.40} 0.00 
0.4 .098) 0.95 .098| .093} 10.2 0.38 1.55| 0.18 
19 | 14.5 -100] .087| 0.94 -100] .093| 13.6 | 0.71 6.74, 1.62 
Ft... 34 | 19.3 0.93 -102| .094| 27.6 | 1.28 7.98} 1.49 
F2.. 34 | 30.0 .092}| 0.93 .103| .096| 20.4 1.36 5.79} 1.68 
F3....| 49 | 39-5 | -104} 0.94 -104 098] 30.5 | 2.55 | 11-65] 3-77 
FA: .107| .109| 0.95 .105| .102| 30.8 | 1.64 8.88} 2.10 
BS... 60 | 48.0 .105| .097| 0.98 -106| .105| 32.7 | 2.22 | 12.73 2.62 
F6. 59 | 49-3 .102| .105] 1.00 -106| .109| 28.9 | 1.05 | 10.84) 1-59 
2509 .105| .115] 1.02 -107| .113| 17.4 | 0-88 6.88} 0.91 
40 | 35:3 .100] .105| 1.04 107} .116| 21.0 | 1.78 12.48] 2.94 
Fo.. 27 | 20.2 -¥05| 2-05 -108| .118| 12.2 | 0.91 5.72) 0.72 
Go AY} 37-0 .118] 1.04 .117| 16.2 | 0.92 8.15} 1.68 
Gi 39 | 29.7 .109] .122| 1.04 .116) 18.0 | 0.69 6.04} 1.13 
G2 3229 .109| 1.03 -108| .114| 21.3 | 0.64 7.42| 0.82 
G3 35 | 28.3 093} 1.03 .112| 16.6 | 0.44 5.64] 0.32 
G4 27 | 21.6 .110} 13-4 | 1-39 | 5-57 1.84 
G5 .113| .118) 1.02 .105| .108| 16.8 | 2.50] 9-27 2.89 
GO a9 | .100|  .097| I.O1 .103| .106] 14.1 | 0.41 5.19| 0.80 
G7 26 | 26.3 .097| 1.00 .105) 11.6 | 1.04 6.16} 1.79 
G8 21 | 19.4 .103] 1.00 .103] 10.6 | 1.01 5.30) 
Go Fer 16.9 .096| .072] 0.99 .102| 6.4 | 0.89 4:41, 1.20 
Ko 27.1 40:3 .097| .103} 0.99 .100] 11.8 | 0.74 6.60} 0.96 
Kr 4 -102} .109| 0.98 .0gg} 14.6 | 0.71 4.52) 0.31 
K2...| 23 | 24.8 .094| .092| 0.97 100] .097| 11.3 | 0.33 | 4-50] °-32 
| 2 .105| .101| 0.96 100} .096| 5.9 | 0.28 1.04} 0.07 
Kq4...| 16 | 18.4 .086) 0.95 100} .094| 7-7 | 90-2! 1.33] 0.04 
40 | 53.0 .097} .096| 9.95 .093| 15.6 | 1.07 2.95 
47 | -098| .og1| 0.96 too} .092| 20.9 | 1.82 7.82) 1.65 
Mo 33 | 48.6 .103} .099| 0.97 103} .099| 8.9 | 1.30 | 4.26 1.26 
M1 26 | 40.1 .100] 1.00 tot} .100| 6.0 | 0.88 5.09} 0.85 
M2...| 18 | 30.0 .107| 108] 1.03 107} .108} 2.6 | 2.35 2.15| 1.65 
M3...| 12 | 20.8 .102| .119| 1.07 .119| 1.1 | 0.80] 2.69] 1.07 
M4.. 4 =. 103] 1.15 -099| .103| 0.4 | 0-34 | 0-32) 0-34 
Ms.. 7 | 12.6 100) 120}... 0.106] 0.129} 0.4 | 1.14 1.42) 2.08 


the most. The values of /’ are more ragged, on account of the ena 
of the trigonometric parallaxes. They run definitely below 0” 100 
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near Fo, and high between F5 and Gs. The smoothed values s, and 
f, agree with s’ and ?’ within the obvious scatter of the data. 

For the M dwarfs the individual trigonometric parallaxes are of 
high percentage accuracy and the means have not been smoothed. 

The correction to the assumed spectroscopic absolute magnitude 
which is indicated by the trigonometric parallax is AM = 5 log 
(t,/s 9). A smooth curve, drawn to represent the corrected values, 
minor fluctuations being disregarded, gives results which may be 
called S,. 

The Mount Wilson spectroscopic parallaxes of main-sequence 
stars may then be reduced (on the average) to Schlesinger’s trigo- 
nometric system, by multiplying by a factor R, where log RK = 
5(S,—S,). This is given in the sixth column of Table 3. It differs so 
little from unity that it may be ignored, except in refined studies. 
For the fainter M dwarfs it is greater, but here the trigonometric 
parallaxes are far more accurate than spectroscopic parallaxes can 
be. 

It will be noted that the values of R do not exactly equal ¢,/so. 
The curves defining the latter were drawn quite empirically. When 
combined, they indicate ‘‘stillstands’’ near A8 and Go in the de- 
crease of absolute brightness with advancing spectral type, which 
seems improbable. The values of R have been adjusted to give a 
more nearly uniform slope. Whether this adjustment has been over- 
done can hardly be determined from existing data. 

The values s, and ¢, were used in the discussion described below. 
The small changes which would be required to fit the (later adopted) 
values of R would have no sensible effect on the conclusions. 

The general weighted mean of all the values of ¢’ for the whole 
main sequence is 0. 1034—agreeing perfectly with that of s’. 

The data for the giants are given in Table 4. Since these stars are 
much more distant, the percentage errors of the trigonometric paral- 
laxes are large and the weights small, especially for class M. The 
irregularity in the values of s’ arises mainly from differences in the 
proportions of supergiants and subgiants in the various groups. For 
/’, accidental error explains the large variations, especially for the 
small groups, although the small values for the early G’s may be sig- 
nificant. These fluctuations and the uneven distribution of the 


| 
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supergiants make it impracticable to draw mean curves for s, and éo. 
The weighted means for each of the four large groups into which 
the stars were combined have been adopted. The weighted mean of 
s’ for all the 732 stars is o%0949, and of ¢’ also o%0949. The zero 
points of the trigonometric and spectroscopic series are therefore in 
exact agreement. There is no evidence of difference for the various 
spectral groups, except possibly for the early G stars. 


TABLE 4 
OBSERVATIONAL DATA: GIANTS 


[p|B 
SP N (p] s’ v So ty NE | NE VF 
100 

Go-G1. 13 8.2 | 07064) 07044] 07079] 18.7 | 1.72 | 3.85] 0.79 
G2-G3. 24 Il.4 .064] .057| .079| .069] 38.4 |] 5.40 | 16.48] 4.52 
28 8.4 .104; .079] .069} 50.0 | 2.58 | 15.90] 3.98 
G5 32 .092] .074| .079] .069] 54.5 | 1.42] 9.55] 1.83 
Go 42 .10Q} .103] .107] 29.8 | 2.58 7.58) 1.64 
38 13.8 .103] .103) .£07| 27:3 | 2.98) 10.45) 3.47 
G8... 40 13.5 .097} .117| .103] .107| 29.0 | 1.28 6.96] 1.45 
Go.. 26 9.3 .093| .082] .103] .107] 18.6] 0.29] 4.77] 0.66 
58 20.0 .102] .099] .103] .107] 41.8 | 3.92 | 10.38] 3.40 
| 69 22.0 .107| 50.6°} 3.66 57-708 
52 12.8 .103] .116] .096} .098} 49.5 | 2.08 | 12.46] 1.58 
§2 12.6 .100} .102} .096} .098] 49.6 | 0.82 | 13.07) 1.64 
|. 45 9.9 .090] .100] .096} 43.5 | 1.00 | 10.36] 1.54 
70 18.5 .0go} .o80} .096} .098} 65.9 | 1.23 | 14.00] 1.06 
32 .102] .106} .085] .085] 39.4 | 0.73 | 8.26] 0.70 
21 4.0 .0068] .104] .085} .085} 26.3] 1.16 | 4.72] 0.58 
34 7.2 .076} .051| .085} .085| 42.2 | 1.69 | 6.00} 1.30 
27 3.8 .093| .102| .085} .085] 34.8 |] 0.14 7.42] 0.14 
20 1.4 .098] .112| .085} .085| 26.7 | 0.10 3.12] 0.29 
Ms-Mo6... 9 I.2 | 0.072] 0.016] 0.085} 0.085] 11.6 | 0.17 3.80} 0.48 


I1. ANALYSIS FOR CORRELATION 


We may now apply equations (23) to (26). The small groups of 
Tables 3 and 4 were at first treated separately. The results were 
often ragged, and some values of Y?, and more rarely of W?, came 
out negative, showing that the analysis was vitiated by the acci- 
dental accumulation of accidental errors. The results here presented 
depend on means for groups large enough to be statistically signifi- 
cant. The data for the subgroups were combined as described in 
section 6. Table 5 gives, for these groups, the quantities which ap- 
pear in equations (24) to (26) and the derived values of /, W,, and 
Y,. Following the seven groups of main-sequence stars are means for 
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two larger groups, and for the whole. Similar data for the giants ap- 
pear in the lower part of the table. 


For the main-sequence stars as a whole, the effect on the reduced 
parallaxes of the real dispersion of their absolute magnitudes (W?) 
is somewhat greater than that of the errors of the trigonometric de- 

TABLE 5 
CORRELATION ANALYSIS 
Sp N |W3+BxX2| Bx2 | w? | ww? W, | Yo 
Main Sequence 

A1-Ag....| 124] 84.6] 0.250 0.048 10537 
Fo-F4.....]| 187} 141.8] .289 .149] .140] .075 .053 
398.5) .108} .164] .049 .038 . 30 AO) 
Go-G4....| 186] 149.9] .219 .098] .121] .039 .027 32 35 38 
Gs-K1....| 178) 171.2] .242 .087| .155] .054 043 35 .39 | .44 
K2-K6....] 138) 168.1] .132 .063} .069] .023 .022 33 
Mo-M5...} 100] 159.3] 0.100 |0.021/0.079]0.039] 0.043 | 0.50 | 0.28 jo.31 
A1-Fo..... 538] 404.9] 0.274 |0.135]0.139]0.056] 0.046 0.40 | 0.37 |o.38 
Go-Ms.. 602] 648.5] 0.173 |0.067/0.106}0.039] 0.034 | 0.36 | 0.33 |o.38 

All 1140]1053.4] 0.212 |0.093/0.119|0.046] 0.038 | 0.383] 0.344/0.378 

Giants 

G1-Gs.... 97| 40.1] 1.140} 0.696/0.444/0.277} 0.277 | 0.62 | 0.67 |o.52 
GG-Kr....| 194.3] 2252] .166 .156 | 0.66 50 
K2-K5....] 219} 53.8] 0.927) 0.670] .257} .108 .095 | 0.42 
Mo-M6...} 143] 24.7] 1.349] 0.162 |(1.73)|(0.29)]..... 
Gi-K1....| 370] 134.4] 0.771] 0.462]0.309/0.199] 0.192 | 0.64 | 0.57 |o.39 
K2-Mo....| 362] 78.5] 1.060} 0.858]0.20210.119] 0.116 | 0.59 | 0.45 |0.37 

All....| 732] 212.9] 0.877] 0.608]0.269/0.170] 0.164 | 0.63 | 0.52 |o.38 


terminations (BX2). The latter is smallest (though still sensible) 
for the M dwarfs. For the A stars it is twice as great as the real dis- 
persion. For the giants, the trigonometric errors are relatively much 
greater (though no worse in seconds of arc); and for the giant M’s 
they swamp the real dispersion so heavily that its value can hardly 
be ascertained in this way, and the computed value of Y?2 comes out 


negative. 


} 
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The resulting values of W, for the various main-sequence groups 
are in tolerable agreement. The general mean, 0.344, is the standard 
deviation of exp 0.46 (M — M,) divided by its own mean value. 
If M — ./, has a normal distribution, the corresponding standard 
deviation of M is +0.72 mag.® 

For the giants, the values of W, are greater (except the ill-deter- 
mined one for Mo-M6). The general mean, +0.52, corresponds to a 
standard deviation of 1.07 mag. in M. This larger result arises from 
the inclusion of supergiants and subgiants. Its numerical value de- 
pends upon the proportions of these which happened to be included 
in the present list, and has no physical significance. 

The values of Y, are more ragged than those of W,—as might be 
expected, since its determination involves taking another difference 
of squares, followed by division by /, which is none too precisely de- 
termined. Its mean value for the main sequence corresponds to a 
standard deviation +o0.80 mag. in M. The mean for the giants is 
almost identical. 

This is a much greater error than has usually been attributed to 
spectroscopic determinations of absolute magnitude. The reason has 
been explained in section 3. The values previously adopted were the 
apparent dispersions, /Y,. For the mean of the main-sequence stars 
this is 0.145, and for the giants 0.240, corresponding to standard de- 
viations +0.32 mag. and +0.51 mag. in M. 

The factor / averages considerably smaller for the main-sequence 
stars than for the giants, mainly because the real dispersion W, is 
smaller compared with the observational error Y,. 


12. STATISTICAL ERRORS 


The statistical errors which might be expected, on a basis of ran- 
dom sampling from a homogeneous population, are given in Table 6 
for the principal quantities derived in Table 5. The coefficient 


+ Y?) 
SIf y = exp ax, and x has a normal distribution with standard deviation ¢, 
= exp (y— 9)? = exp 2070? — exp a’o?; 
[o(y)/y)?? = exp ato? — 1. 


} 
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is the slope of the regression line (8) and 1// that of the “impartial 
line (10). 

The anticipated mean errors were computed by equations (37) to 
(39) and converted into probable errors. The errors given for the 
larger groups are derived from the sums of [pvv]/s?, etc., for the 


TABLE 6 


STATISTICAL PROBABLE ERRORS 


Sp Wi a r/l 
Main Sequence 

a7 + .o7 1.86+ .18 
0.28 +0.012 | 0.92+0.08 2.02+0.19 
0.37 | 2.47+0.16 

©.344+0.009 | 1.19+0.04 2.61+0.11 

Giants 

0767 £0.13 1.00+0.12 1.60+0.32 
so + .07 1.062 <07 .16 
Mo-M6............. 0.29 +0.66 0.87+0.13 0.58+0.72 


separate smaller groups and do not include the effects of variation 
of s, and ¢, from group to group (which, however, are small). 

For the main sequence the results for each group are statistically 
significant, and also for the giants, except for the M’s. Although this 
group contains 143 stars, the parallaxes are so small that the errors of 
sampling deprive the deduced values of W, and 1// of meaning. 


| 
| 
| | 
| 
| | 
| | 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 


COMPARISON OF PARALLAXES 415 


The value of a, which is much less influenced by these errors, is sig- 
nificant. 

The range of the results for the main-sequence groups shows that 
influences other than errors of sampling cause them to differ. If the 
seven values for W, are treated as independent determinations, with 
weights as in Table 5, the general mean is 0.340 + 0.016 (P.E.). 
The mean value is substantially the same as before, but the probable 
error is 1.7 times greater, indicating that the real differences (meas- 
ured by their squares) are twice the statistical. There is a fairly 
steady decrease in the dispersion of W, (or of the absolute magni- 
tude) with “later” spectral type, except for the group A1—Ag, which 
may not be typical of all stars in this spectral range, as many of them 
are ill-suited to spectroscopic determination of parallax by the 
methods employed for the redder stars. 

For the regression constant a, the weighted mean is 1.18 + 0.06. 
Sampling error contributes two-fifths of the squares of the discord- 
ances. This may account for a good part of the raggedness of the 
individual values; but the remainder seems pretty large to be ex- 
plained by differences in the calibration-curves, which, for the same 
pair of lines, are usually nearly parallel for different spectral classes. 

For 1/1 the weighted mean is 2.69 + 0.16. Once again, the mean 
value is not sensibly affected; but the probable error is much in- 
creased. The mean for the separate values of / is 0.390 + 0.025. 
The weighted mean for Y, is 0.377 + 0.025. The error is greater 
than for W,, as was to be expected. The corresponding statistical 
expectation was too troublesome to compute. 

For the giants, the weighted means and probable errors are W, = 
0.51 + 0.04,a = 1.04 + 0.03, 1/1 = 1.63 + 0.20. Two of these three 
errors are less than the anticipated values, doubtless because of acci- 
dentally favorable grouping of the errors. 


13. FINAL REDUCTION FORMULAE 


In discussing these results, we must remember the two types of 
calibration, of which type I seeks to obtain agreement between the 
mean spectroscopic and trigonometric parallaxes for stars grouped 
according to the spectroscopic absolute magnitude; while type II 
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seeks this agreement for groups selected by some impartial criterion, 
unaffected by the errors of either sort of observations. 

As regards type I, the equality of the mean values of s’ and /’ 
(Table 4) and also the values of a (Table 6) show that for the giants 
the calibration employed in Mount Wilson Contribution No. 511 re- 
quires, in the mean, no correction. Possible differences in the cor- 
rections for different spectral types are too small to be revealed by 
the trigonometric parallaxes. Since the latter have the lowest per- 
centage accuracy for giants, and the published calibration has been 
largely based on other weighty data, such as proper motions and 
radial velocities, this agreement might fully as well be taken as evi- 
dence of the correctness of the system of Schlesinger’s parallax cata- 
logue as of that of the spectroscopic parallaxes. 

For the main sequence, however, the trigonometric parallaxes 
constitute the strongest existing data. Here there is no evidence of a 
general difference of zero point, but the mean value of a shows con- 
clusively that the calibration would be improved if the differences 
s’ — s’ were increased by 19 per cent. It would be still better to ap- 
ply this correction to the absolute-magnitude differences S — S, 
(sec. 2). 

If the spectroscopic parallaxes should be thus adjusted, the agree- 
ment for groups selected by spectroscopic criteria (method I) would 
be exact in the mean. Should this correction be made, the mean 
value of / would be increased from 0.383 to 0.456. The latter value 
represents the extent to which the calibrated differences in absolute 
magnitude are diminished by the method depending on the regres- 
sion line. 

It is doubtful whether any real improvement would be obtained 
by using the individual values of a given in Table 6 for the different 
groups, and the value a = 1.19 may well be used for the whole main 
sequence. The corrections to the zero point indicated by the factor R 
(Table 3) should, however, be made in precise work. 

In order to obtain a calibration of type II, much larger corrections 
are necessary. The zero point may be treated as before; but the dif- 
ferences from the mean for the spectral subclass should be multi- 
plied by 1.6 for the giants and by 2.6 for the main sequence. It is 
essential here to apply this correction to the absolute magnitudes; 
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to treat the parallaxes directly would result in negative parallaxes 
for supergiants brighter than absolute magnitude —2. In view of 
the importance of other data in the calibration for the giants, it is 
possible that a smaller correction factor might secure a better agree- 
ment, for impartially selected groups, between mean spectroscopic 
parallaxes and those from all geometrical methods (including proper 
motions, etc.). The still larger correction factor for the main se- 
quence is, however, strongly determined, and should undoubtedly 
be applied. The evidence hardly justifies using factors varying with 
the spectral type, except as regards the zero point. 

Following equation (10), we have, for the impartial calibration 
line, in the case of the reduced parallaxes, 


S — So 


l 


Let s’’ be the reduced spectroscopic parallax so modified as to agree, 
in the mean, with the trigonometric magnitudes for large impartially 
selected groups. Then 


l 
where 
— to) 
= So t+ 
With 1/l = 2.6, 
s'=s’ .6 (s' — 5; 
Sr = So + 0.62(55 — 


We have also 

tos’ = expo. 46 (S — Sp), 
where S is the tabular spectroscopic absolute magnitude, and S, 
the standard value (Table 2). If S’’ is the absolute magnitude re- 


calibrated according to method II, and S, that corresponding to the 
reduced parallax s,, we have, to the first order, 


10s’ = 1+ 0.46(S — S,), etc., 


' 
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and equations (41) become 
=S§+1.6(S — §,). (42) 


For the main-sequence groups we find as follows: 


SPECTRUM 
ALL 
A1-Ag | Fo-F4 | Fs-Fo | Go-G4 | Gs-K1 | K2-K6 | Mo-M5 
©.105| 0.103] 0.107] 0.108] O.I0I] 0.100} 0.103] 0.1037 
O.IOI] 0.097} O.III] O.114} 0.103] 0.094} 0.106] 0.1040 
0.108] 0.107] 0.104] 0.104} ©O.100] 0.104] O.I0I] 0.1035 
+o™r7 |+o0™1r5 |+0™o0g |+0™og | oMoo |-+0™og |+0™o2 |+0™08 


When the reduced parallaxes have already been calculated, equa- 
tions (41) are more convenient, as they may be applied directly to 
the mean for a group. When individual values are desired, it is 
quicker to use equation (42) or the equivalent form 


log (5) = 0.32 (5 (43) 


which gives the spectroscopic parallax directly. 
For the giants we have similarly, taking 1// = 1.60: 


SPECTRUM 
ALL 
G1-G5 G6-Kr K2-Ks5 Mo-M6 
0.069 0.107 0.098 0.085 0.0950 
0.096 0.096 0.093 0.085 0.0939 


It is obviously sufficient to use the mean values throughout, so 
that 


sl! 


= 5’ + 0.6 (s’ — 0.094) , 
S” = § + 0.6(S — S,), 


” 


log =o0.12(S — 
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14. NUMERICAL TESTS 


As a practical test of the validity of these formulae, the stars of 
each spectral subclass were arranged (a) in order of the reduced 
spectroscopic parallaxes s’, and (b) in order of the similarly reduced 
proper motions, nu’ = fu; and weighted means of s’ and ¢’ were taken 
for four groups of substantially equal numbers. The first grouping 
should give a correlation of type I. If T is the tangential velocity of 
the star (in A.U./yr), we have 


log (10n’) = log T ++ 0.2(M —S,). 


Grouping by yp’ is independent of the spectroscopic errors, and there- 
fore of type II. Despite the large variations to be expected in T, 
it should give a fairly good separation according to M. 

The results for the main sequence are given in Table 7. 

These results are plotted in Figure 1, with the means for grouping 
by wu’ as dots, and those by s’ as crosses, and also the “impartial” 
line (10) and the regression line (8)—the latter dotted—as deter- 
mined by the constants for each group (Table 6). 

The agreement of the means by s’ with the regression lines is in 
all cases close, as it should be, since these are derived from the same 
data. For the means according to yu’ the agreement is unforced. 
Though the results are ragged for some groups, the points invariably 
indicate slopes much greater than those of the regression lines, and, 
in general, agree with the impartial lines as well as might be ex- 
pected. 

A least-squares solution for the four groups in the last column, 
with equal weights (which suffices), gives 


= 2.29(s' — 5’). 


The residuals indicate a probable error + 0.033 for this coefficient— 
which is doubtless too small. 

From these diagrams, as well as from the analyses previously dis- 
cussed, it appears beyond all doubt that the differences of the spec- 
troscopic absolute magnitudes (or the reduced spectroscopic paral- 
laxes s’) from the normal values for the spectral subclass concerned 
require a small increase to obtain the best calibration by the regres- 
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sion method, and a very large increase, by a factor of about 2.5, to 
represent the deviations in groups selected impartially. 


TABLE 7 
GROUP MEANS: MAIN SEQUENCE 


Group A1-Ag | Fo-F4 F5-Fo | Go-G4 | Gs-K1 | K2-K6 | Mo us| All 
Grouping by s’ 
1........] 29 44 58 43 40 33 23 270 
s 0.077 | 0.075 | 0.078 | 0.089 | 0.080 | 0.082 | 0.079 | 0.080 
od 0.078 | 0.060 | 0.074 | 0.074 | 0.076 | 0.073 | 0.086 | 0.074 
2 NV 32 48 57 48 47 35 27 204 
s” ©.099 | 0.100 | 0.100 | 0.103 | 0.096 | 0.092 | 0.093 | 0.098 
t ©.100 | 0.090 | 0.110 | 0.104 | 0.099 | 0.087 | 0.096 | 0.098 
$i: N 32 49 57 48 40 35 26 293 
| | DIO | 146-106 | | 0.207 | 0.100 
r 0.106 | 0.108 | 0.112 | 0.125 | 0.110 | 0.099 | 0.104 | 0.109 
, eee N 31 46 55 47 45 35 24 283 
Grouping by yp’ 

Besant N 20 44° 56 43 40 33 23 208 
3 ©.095 | 0.094 | 0.090 | 0.104 | 0.093 | 0.092 | 0.094 | 0.004 

v 0.068 | 0.074 | 0.082 | 0.092 | 0.086 | 0.085 | 0.089 | 0.084 

RE 4‘ 32 48 57 47 40 34 26 290 
s’ 0.105 | 0.098 | 0.103 | 0.110 | 0.102 | 0.099 | 0.100 | 0.102 

0.095 | 0.087 | 0.108 | 0.115 | 0.101 | 0.091 | 0.107 | 0.101 

g.. N 32 49 57 47 47 35 27 204 
s’ ©.104 | 0.113 | 0.114 | 0.115 | 0.098 | 0.099 | 0.104 | 0.106 
| 0.114 | 0.128: 122. | 0. | 0-110 

4 N 31 46 55 48 45 34 24 283 
oy 0.119 | 0.118 | 0.109 | 0.116 | 0.116 | 0.104 | 0.114 | 0.113 


It is not worth while to make a similar graphical treatment of the 
giants, since the percentage errors of the trigonometric parallaxes are 
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so much larger. The results of Tables 5 and 6 show, however, that in 
this case also the spectroscopic method is successful in finding not 
only the mean absolute magnitudes for the spectral subclass but the 
deviations from this mean. The correcting factors are, in fact, con- 
siderably smaller than for the main sequence. 


15. ERRORS OF TRIGONOMETRIC PARALLAXES 


One result of the present investigation has been to confirm the 
reliability of both trigonometric and spectroscopic determinations 
of parallax. This is conspicuously the case for the trigonometric 
parallaxes of giants. Examined individually, the determinations of 
different observers for the same star differ by so large a percentage 
as to raise doubt as to their significance; but the detailed statistical 
discussion shows that they are substantially as good as they claim to be 
(that is, as is indicated by the probable errors given in Schlesinger’s 
Catalogue). As a further test of this, an analysis was made of the 
trigonometric determinations for the hardest case, the supergiants. 
Trigonometric parallaxes are available for 41 stars in the Mount Wil- 
son list of classes cGo to cKs5, and for 16 others (all but two of 
class M) whose spectroscopic absolute magnitudes (brighter than 
—1.5) mark them as objects of similar character. 

It was assumed that all these stars were of the same absolute 
magnitude (adopted as —2.1 after a preliminary discussion). If s, 
is the parallax corresponding to this assumption, and ¢ the observed 
trigonometric value, the residuals ¢ — s, should exaggerate the trigo- 
nometric errors. For each star, this ratio (¢ — s.)/r was taken (where 
r is the probable error given in Schlesinger’s Catalogue). The individ- 
ual values of this ratio range from —3.8 to +3.5. Their arithmetic 
mean is — 0.13, and their mean square value, 1.46, as against 
(1/0.6745), or 1.48, to be expected in a normal distribution. The 
observed parallaxes then agree with the assumed value slightly bet- 
ter than would be expected from their catalogued probable errors 
alone, leaving a negative value for the effects of real differences in 
absolute magnitude. 

The statistical probable error of a determination of probable error 
for 57 entries is, however, +9 per cent, so that no weight is to be at- 
tached to this result. The probable error of the mean of 57 deter- 
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minations is +0.137, so that the difference of the arithmetic mean 
from zero is also without significance. 

It is evident, however, that the tabular probable errors of Schles- 
inger’s Catalogue are, on the average, reliable measures of the accu- 
racy of the observations. 


In conclusion, it is a pleasure to acknowledge our indebtedness to 
Messrs. Adams and Joy for unpublished data regarding the calibra- 
tion-curves (sec. 8) and the distribution of the tabular parallaxes 
(sec. 10); and also to Miss Myrtle Richmond for material aid in the 
calculations described in section 14. 

CARNEGIE INSTITUTION OF WASHINGTON 
MovuntT WILSON OBSERVATORY, 


PRINCETON UNIVERSITY OBSERVATORY 
November 1937 
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THE PHOTOGRAPHIC DETERMINATION OF 
STELLAR PARALLAXES WITH THE 60- 
AND 100-INCH REFLECTORS* 


SIXTEENTH SERIES 
ADRIAAN VAN MAANEN 


ABSTRACT 

Trigonometric parallaxes determined with the 60- and 100-inch reflectors.—Results 
are given for 25 fields including 28 stars, 15 of which have photographic absolute magni- 
tudes of +10 or fainter; 3 are even +15 or fainter. One white dwarf is included, 
AC+70°8247, for which the absolute magnitude is +11.9. For U Geminorum and 
SS Cygni the absolute parallaxes are +o0”’o10 and —o"o12, respectively. 

The present paper gives the results for 28 stars in 25 fields,which 
have recently been measured for parallax. The work has been car- 
ried out in the same way as for the previously published series.’ The 
results are given in Table r. 

Several of the stars require comment. The first star in the list was 
found by Kuiper to be a close double, A = 0770, Am = 0.3.? The 
components have absolute magnitudes of +13.1 and +13.4, re- 
spectively. 

BD+23°123 is an R-type star with, according to Sanford, a 
velocity of 232 km/sec and dwarf characteristics in the spectrum. 
Using the mean of the parallax determinations by McCormick and 
Mount Wilson, we find an absolute visual magnitude of +5.2. 

U Geminorum and SS Cygni were put on the program as the two 
examples of their type of variation which have fairly large motions; 
from these Parenago and Kukarkin? derived parallaxes of 07016 and 
0”038, respectively. The absolute parallaxes derived here are 
+o%oro and —o%o12, and it is doubtful if the stars should be con- 
sidered as dwarfs. 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 590. 

t Mt. Wilson Contr., Nos. 111, 1916; 136, 1917; 158, 1918; 182, 1920; 204, 1921; 
237, 1922; 270, 1923; 290, 1925; 321, 1926; 356, 1928; 391, 1929; 435, 1931; 468, 1933; 
506, 19353 553, 19306. 

2 Pub. A.S.P., 47, 282, 1935. 

3 Verdnderliche Sterne, V.F.A., Gorki, 4, 249, 1934. 
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A star close to Ross 110 (A = +5”) was found to have the same 
large proper motion of 0”70. Its magnitude is about 0.4 fainter than 
that of Ross r1o. 

Near OF 285 (A = 50”) a star was found which shares the motion 
of the principal double; apparently this is a triple system. The mo- 


TABLE 1 


| Relative No. In 
Object a 1900 6 1900 Wilson Paral- | uw, Rel. | P.E. of stru- 
ag. Spec- = 
lax Exp. ment 
trum 
Anonymous. OF 7428 |-+-68° 47’ 113.7 |.....- +o0"056+0"7704 7 10 18 60 
BD+23°123.... 0 48 54 +23 32 | 8.8t R3 — .005+0.139 6 1 18) 60 
2 742 \+47 +1 6.31| |+ .o12i—0.078 2 18 60 
Ross 30 10 |...... + .046+0.339 5 2, 20) 100 
U Geminorum.. 7 49 10 +22 16 Var.|...... + .o11—0.042 7 10! 18 60 
vM- W 48*..... IO 24. 42 + .o12+0.018 5 2. 18) 100 
BOBS: 9 49 45 |+22 58 |14.2 |...... + .006+0.003 3. =I 16 | 100 
70.30: Sie + .o12+0.254 8 3); 20) 60 
Ross 110 B..... WE + .035+0.55217. 16! 60 
OF 285 A+B...14 41 42 +42 48 | 7.66 F5 + .029—0.05510 2, 16); 60 
OF 285 14 41 37 |\+42 48 |12.6|...... + .co1—0.050 6 =I 14| 60 
Radcl. 60, 600. .|15 57 46 |+30 27 |14.1 |...... + .o59—-0.311 3 5, 18. 60 
OZ 323. ...(17 2 18 |-+47 9.08 K3 |+ 20} 60 
AC+70° 8247, 70 46: 113.2 + .047+0.064 4 7 220; 60 
vM- W 24*.....\19 59 18 |+29 36 |16.5 |...... + .084+0.682 6 2. 100 
vM- W 71*......20 — .003+0.068 4 I 16 100 
vM- W 72*...../20 12 24 |+15 8 |18.3 |...... — .005+0.064 5 2 16. 100 
Ross . 120 27 42: + .063\+0.207'4 7 20, +60 
ADS 14314A...'20 44 13 55 | 7.43) cA7e — .017+0.005 5 1. 18) 60 
Luyten 983. ....|22 33 20 |— 24 35: |... + .066+1.00810 3) 20) 100 
Ross 201...... 6: + .041+0.369 8 10, 20. 60 
vM- W 112*....|\21 37 10|— O 4 |17.38...... + .033—0.21312 20/ 100 
SS Cygni. ...... 21 38 47 +43 8 — .o17+0.112 8 8 20. 60 
Anonymous....|2% §% 2 S4 .018 —0.332, 5 20! 60 
Lal, 22 19 29 +38 4 F8t 025+0.25410 2 16: 60 
*From Mt. W. Contr., No. 412, 1930. t Visual magnitude. 


{Spectrum from Henry Draper Catalogue. 


tions in declination were also measured; and the results are ws = 
+ o"051 and +0%052, respectively. 

Ross 508, vM- W 24, which is a distant companion of Lalande 
38380 (A = 3’), and vM- W 112 have absolute photographic magni- 
tudes of +15.1, +15.4,and +15.0, respectively. Mount Wilson Con- 
tribution No. 507% gives a list of 12 stars whose photographic abso- 


4 AD. J., 8%5.1§2; 1035: 
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lute magnitudes are +15 or fainter. The companion of Procyon 
should be removed from this list, since Kuiper has recently found it 
to be 2.7 mag. brighter than as given there. The companion of Wolf 
860 also should be omitted, as the absolute magnitude is +14.9. The 
stars given in the accompanying table should, however, be added to 


Object a@ 1900 5 1900 B Pg.m Pg. M 
| +22°29’ | +48°] 15.7 | 0.72 | +0.089 | +15.4 
ROSS COS. o.555 15 19 18 +17 51 +52 | 15.4 |] 1.23 | + .088 | +15.1 
Companion of La- 

lande 38380. ..} 19 59 18 +29 36 — 2] 16.5 | 0.92] + .o61 | +15.4 
vM- Wilz2...... 21:37 10 —-0O4 —38 | 17.4 | 0.45 | +0.034 | -15.0 


the list. We therefore know at present 14 stars with photographic 
absolute magnitudes +15.0 or fainter; three of these stars are com- 
panions. The star mentioned by Luyten in Astronomical Journal, 
No. 1070 (46, 133, 1937), may fall in this class too, but the prelimi- 
nary parallax derived by Luyten is too uncertain to include the star 
at present. 

OZ 323 is apparently an optical double; the change in distance 
and position angle shows a relative motion in a straight line, without 
any indication of orbital motion. 

AC+70°8247 is the white dwarf, found by Kuiper and discussed 
by him in Publications of the Astronomical Society of the Pacific, 47, 
307, 1935. The mean of the absolute parallaxes derived at Green- 
wich and here is +0%057, giving an absolute photographic magni- 
tude of +11.9. 

Although the proper motions in a for vM- W 71 and 72 as pub- 
lished in Mount Wilson Contribution No. 412 were +0%079 and 
+o"118, respectively, it should be noted that the stars were 21’ from 
the center of the plates then measured, which were taken at the 
Newtonian focus of the 60-inch reflector. The present measures give 
+0068 and +07064, respectively; and the pair is probably a dis- 
tant system of the 61 Cygni type. 

BD +48°4235 =ADS 14314, a star of the a Cygni type. The 13.7 
companion at A= 273 was not measurable, as it was covered by the 
sector. 
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The motion and parallax of Ross 201 are very much the same 
as for Ross 200; this pair, also, may therefore be a system of the 
61 Cygni type: 

Ross 200. ..... 21536™33 +53°41' w=o"60in p=64° 
Ross 701. ..... 2137 6 +53 33 0.55 75 +0.045 


In recent years I have tried to determine the parallaxes of those 
stars which presented a fair chance of having faint absolute magni- 
tudes; and of the stars on my program thus far finished, 114 fall 
under this heading. Faintness has been no objection, since with the 
1oo-inch reflector we can, with somewhat increased exposures, reach 
the eighteenth photographic magnitude. Altogether we have found 
but five negative parallaxes among the 114 stars; they include, more- 
over, 15 companions. Omitting both these groups, we find the abso- 
lute magnitudes to be distributed as shown in the accompanying 
table. 


M n M n 

20 


The maximum frequency falls at about +12.7 photographic or, 
as practically all these stars are probably of type M, at +10.9 vis- 
ual. If now we look at Figure 1 of Mount Wilson Contribution No. 
562,5 it is evident that, of the stars with large proper motions as a 
whole, we know only the brighter ones. The expected maximum is 
about 2 mag. fainter than the observed maximum. We may there- 
fore conclude that, if we knew also the absolute magnitudes of the 
fainter stars of large proper motion, the maximum frequency would 
be shifted about 2 mag. and therefore be + 12.9, a good confirmation 
of the result derived in Mount Wilson Contribution No. 562.5 

CARNEGIE INSTITUTION OF WASHINGTON 


Mount WILSON OBSERVATORY 
January 1938 


5 Ibid., 85, 26, 1937. 
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ON THE POSSIBILITY OF DETECTING MOLECULAR 
HYDROGEN AND NITROGEN IN PLANETARY AND 
STELLAR ATMOSPHERES BY THEIR ROTATION- 
VIBRATION SPECTRA 


GERHARD HERZBERG 


ABSTRACT 


The detection of molecular hydrogen and nitrogen in planetary or stellar spectra, 
hitherto deemed impossible, can be carried out by means of the rotation-vibration 
spectrum of these molecules. Though H, and N,, as is well known, have no ordinary 
rotation-vibration spectra (since their dipole moment is zero), they do have rotation- 
vibration spectra, owing to their guadrupole moment. 

In the case of H, the 1—o band of this quadrupole rotation-vibration spectrum, 
according to calculations of James and Coolidge, is 8.1 X 10-9 times as intense as the 
1—o band of the ordinary rotation-vibration spectrum of HCl. The minimum absorb- 
ing layer necessary to detect the I—0, 2—0, and 3—0 bands is found to be 2.5, 2.7, and 
13.0 km atm., respectively. This is of the order of magnitude probably available in the 
atmospheres of the major planets. A table of the positions of the lines of the 1—o. 
2—0, 3—0, and 4—o bands as predicted from the ultraviolet 7, spectrum is given. 
The band most favorable for detection is the 3—o0 band at 8500 A. Failure to observe 
this band would at least give an upper limit for the amount of //, present in the atmos- 
pheres of the major planets or of low-temperature stars. 

For V, the predicted positions of the Q branches of the bands are given. Their de- 
tection will probably be more difficult than the detection of the H, bands. 

A further possibility of detecting molecular hydrogen and nitrogen is by the ordinary 
rotation-vibration spectrum of the isotopic molecules 1D and N*4N'5, which are always 
present in natural hydrogen and nitrogen, respectively. 


I. INTRODUCTION 


It has, up to the present, always been considered impossible to 
detect molecular hydrogen or nitrogen in planetary or stellar atmos- 
pheres." The band systems of H, and N, in the visible and the near 
ultra-violet regions have highly excited electronic states (>6 volts 
above the ground state) as their lower states and consequently can- 
not, in general, appear in absorption. If in a high-temperature star 
the thermal energy would be sufficient to excite these levels, at the 
same time it would be sufficient to dissociate the molecules [D(H,) = 
4.45, D(N.) = 7.35 volts], and again no molecular absorption would 
occur. 

On the other hand, according to Wildt,? Russell, and others, it 
seems necessary to assume the existence of large amounts of molec- 


™H. N. Russell, Nature, 135, 219, 1935. 2 Tbid., 134, 418, 1934. 
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ular hydrogen in the atmospheres of the major planets and also a 
certain amount of molecular nitrogen, as indicated by the presence of 
CH, and NH, in these atmospheres. Also, the atmospheres of the 
cooler stars, according to Russell,’ contain considerable amounts of 
H, and N,. It would consequently be of great interest if it were 
possible to detect H, and N, spectroscopically in planetary and 
stellar atmospheres. 

It is the object of this paper to point out a possibility of detecting 
the presence of sufficiently large amounts of molecular hydrogen 
and nitrogen in planetary and stellar atmospheres by their rotation- 
vibration spectra. 


II. GENERAL DISCUSSION OF QUADRUPOLE ROTATION- 
VIBRATION SPECTRA 

It is generally considered that a rotation-vibration spectrum of a 
homonuclear diatomic molecule, such as H, or N., does not exist. 
This consideration is based on the fact that a homonuclear diatomic 
molecule has no dipole moment. But, apart from an alternating 
dipole moment, an alternating quadrupole moment may also lead to 
emission or absorption of radiation (quadrupole radiation); and 
certainly a homonuclear molecule has a quadrupole moment which 
varies during the vibration of the molecule. It is true that quadru- 
pole radiation is very much less intense than ordinary dipole radia- 
tion, and therefore rotation-vibration spectra due to quadrupole 
radiation have thus far never been observed and have not even been 
discussed in any detail. On the other hand, as is well known, quadru- 
pole radiation is responsible for certain electron transitions in atoms 
like those giving rise to the nebular lines. 

Ordinary infrared rotation-vibration spectra, such as those of 
HCl, can be obtained with absorbing layers ranging from a fraction 
of a millimeter (fundamentals) to several meters at atmospheric 
pressure (higher harmonics). It seems at least possible that with 
layers of many kilometers (reduced to atmospheric pressure), as are 
probably available in the major planets, quadrupole rotation-vibra- 
tion spectra may be observed and, conversely, may be used to prove 
the presence of molecular hydrogen and nitrogen. 


3 Ap. J., 79, 317, 1934. 
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Therefore, it appears to be worth while to derive the positions of 
the absorption lines to be expected and also the absorbing layers re- 
quired. From the data on the electronic band spectra of H, and N, 
the positions of the lines can be derived if the selection rules for 
quadrupole radiation are known. 

Figure 1 gives schematically the rotational levels of two different 
vibrational levels of a homonuclear molecule like H, or N, in the 
ground electronic state, which for both 


: a molecules is a 'D} state. The even rota- 
tional levels are symmetric; the odd are 

, * antisymmetric in the nuclei. The former 
3 a have the larger statistical weight in N,; 
2 s the latter, in H,. There is always the strict 
selection rule that antisymmetric levels do 
not combine with symmetric levels, and 

vice versa. Therefore, a dipole transition 

> vg Oo between the two levels is not possible, since 

for dipole radiation we also have the selec- 

tion rule AJ = + 1. But for quadrupole 

5 @ radiation* AJ =o, +1, +2. Of these, 
AJ =o and AJ = + 2, as may be seen 

' * from the diagram, do not contradict the 
3 a symmetry selection rule. These transi- 
7 ; tions, therefore, may occur as quadrupole 
I e radiation of homonuclear molecules. The 


transitions with AJ = o will form a Q 
branch (but J = o —J = o is forbidden), 
the transitions with AJ = + 2 an S branch, and those with 
AJ = — 2an0O branch. 

Just as the intensity of dipole radiation for heteronuclear mole- 
cules like HCI depends upon the rate of change of the dipole moment 
during the vibration, the intensity of the quadrupole radiation of 
homonuclear molecules depends on the rate of change of the quadru- 


pole moment. 


FIG. 1 


4 Cf., e.g., A. Rubinowicz and J. Blaton, Ergebn. Ex. Naturwiss. 11, 176, 1932. 
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III. THE QUADRUPOLE ROTATION-VIBRATION SPECTRUM OF H, 


The most complete data on the rotational and vibrational levels 
of the ground state of H, are given by Jeppesen.’ From his data the 
wave numbers and wave lengths of the lines of the 1 — 0, 2 —o 
3 — o,and 4 — obands can be accurately calculated, taking account 
of the foregoing selection rules. The predicted values are given in 
Table 1. 

TABLE 1 
PREDICTED QUADRUPOLE ROTATION-VIBRATION BANDS OF HYDROGEN 


1—o BAND 2—o BAND 
| | 
| § Branch | Q Branch O Branch | | § Branch | QBranch | O Branch 
| | | | | 
ACA) | ACA) | ACA) A(A) | A(A) | 
4713 | 21233] | 24060).......]...... 1...| | 12624) | 12386)... 
ya 4917 | 20333) 4143 | 24129 3807 | 26263, 2...; 8782 | 11384 8048 | 12422) 7729 | 12935 
or 5108 19571 4126 | 24229. 3568 28019 3.. 8945 11176 8014 12475 7485 | 13357 
5286 18914 4103 | 24360 3339 620041 4..-, 9090 | 10999 7068 12547 7234 13820 
Cee 5448 | 18351, 4074 24541 3001 32342 5... 9213 | 10851 7OII 12637 6079 14325 
3—o0 BAND 4—0 BAND 
| | | 
| S Branch | Q Branch | O Branch S Branch Q Branch | O Branch 
2....| 12420 8049.4 11725 8526.3 11423 8752.2 2... 15832 6314.8 15176 6587.5 14890 6713.9 
3....| 12556 7962.4 11674 8563.9 11173 8047.9) 3... 15940 6272.0 15108 6617.3 14635 6831.1 
4....| 12667 7892.5 11605 8614.2 10911 9162.7 4... 16018 6241.3 15017 6657.2 14362 6961.1 
5...-| 12753 7839.2 11521 8677.5 10639 9396.8) 5... 16066 6222.5 14905 6707.3 14073 7103.9 


Just as in the case of ordinary rotation-vibration spectra, it is to 
be expected that the intensity of the bands will decrease in the series 
I— 0, 2—0, 3 — 0, 4— 0. As to the absolute intensity of the 
bands, in the case of H, a fairly reliable estimate can be obtained, 
as is shown in a paper by James and Coolidge,° who very kindly 
undertook the calculation of the rate of change of the quadrupole 
moment, using their previous calculations on the ground state of H,. 
From their results it follows that the intensity of the 1 — o H, band 


5 Phys. Rev. 44, 165, 1933. 6 Ap. J., 87, 438, 1938. 
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is about 8.1 X 107° times that of the 1 — o HCl band. The 2 — o 
and 3 — o bands of H, are found by James and Coolidge to have 
o.g1 and o.1g times the intensity of the 1 — o band. Since the 
HCl 1 — o band can be obtained with an absorbing length of 0.01 
cm at atmospheric pressure,? one would conclude that the 1 — o, 
2 — o, and 3 — oH, bands may be obtained with absorbing layers 
of 12, 14, and 65 km atm., respectively. 

But there is one factor that will reduce the necessary layer: in 
H,, even at room temperature, only the first three rotational levels 
occur (at room temperature go per cent of all molecules), whereas 
for HCI about twelve levels are of importance. Thus, for H, the 
intensity is concentrated in fewer lines, and, therefore, the lines may 
be more easily detected than one would expect according to the 
foregoing. This effect is still more pronounced for the planetary 
atmospheres, since the temperature is smaller than — 100° C.! Tak- 
ing T = — 123°C., it is found that the foregoing will give a factor® 
of 0.2, so that 2.5, 2.7, and 13.0 km atm., respectively, would be 
sufficient to produce the bands in absorption. It is difficult to esti- 
mate the accuracy of these figures, partly because of the uncertainty 
of the minimum layer required for the HC/ bands and partly because 
the minimum layer for H, will depend greatly on the spectroscope 
and the dispersion used. The influence of the different line widths 
for H, and HCl has also to be taken into account. One would expect 
the line width for H, to be due to Doppler effect only (the half-width 
at 150° K would be 0.025 A), i.e., to be considerably smaller, than 
for the 1 — o HCl band, for which there is strong collision broaden- 
ing. Considering all this, it seems safe, however, to assume the 
possibility of an error of not more than one power of 10 on either side 
of the preceding estimate. 

In Figure 2 the theoretical intensity distribution in an H, rota- 
tion-vibration band is represented, using the matrix elements given 
by James and Coolidge and assuming T = 150° K. The figure refers 
to the 2 — oband. But, of course, the relative intensity distribution 
for the other bands is the same. 

7 Cf. J. D. Hardy, E. F. Barker, and D. M. Dennison, Phys. Rev., 42, 279, 1932. 

§ This takes into consideration the fact that in HC/ every line has two components, 
due to HCI35 and HCl’, 
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As is seen from Figure 2, in the case of planetary spectra, the 
lines S(o), S(t), and Q(1) are the most likely to be observed; but 
with sufficient absorbing layers, the lines with J = 2, 3, and 4 may 
also be found. The observation of many of these lines in planetary 
spectra is, however, made impossible by other, much stronger ab- 
sorptions of the atmosphere of the earth or of the planet. 

Of the 1 — o band, only the line S(1) might be sufficiently free 
from overlapping by the //,0 telluric band at 1.9 uw and by strong 
methane bands extending from 2.15 to 2.45 u.? But with absorbing 


5(4) S(3) S(2) S(t) 


O(2) O(3)  O(4) 


Q(1) | Q(4) 
Q(3) 
| | | | 


gcoo 8500 8000 7500 


layers of CH, like those on the major planets, it is likely that even 
this line is overlapped by methane bands (for example, by difference 
bands or weak parts of the main bands). Besides, an investigation 
in this region with sufficient dispersion would be extremely difficult. 

The 2 — o band lies in the photographic infrared region, and the 
chances of its being observed are consequently better. The strong 
lines S(o) and Q(1) of this band would lie just between the two 
strong H,O bands at 1.1 and 1.4 uw, and could therefore be observed 
in the spectra of low-temperature stars. But in the spectra of the 
major planets they would be covered by strong methane bands ex- 


9 See J. G. Moorhead, Phys. Rev., 39, 83, 1932. 
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tending from 1.13 » to longer wave lengths and appearing with a few 
meters absorbing length only.” The lines S(3) and S(4) are free 
from H,O overlapping and may possibly not be obscured by CH, 
bands. But, as seen from Figure 2, they would be much less intense 
than S(1) and Q(1) at 150° K. 

As may be seen from a spectrogram of Jupiter reproduced by 
Slipher," there is practically no overlapping for the most intense 
lines S(o), S(1), and Q(1) of the 3 — o H, band. Also, it lies in a 
much more convenient spectral region (8000-8600 A). Therefore, 
except for the fact that this band is about five times less sensitive 
than the 2 — o band, it will be particularly suitable for detecting 
molecular hydrogen in stellar and planetary atmospheres. Since this 
band consists of only a few very widely spaced single lines (with dis- 
tances of more than 100 A), it will not be easy to observe it, especial- 
ly because the line width is extremely small. Unless the amount of 
H, is much larger than the minimum amount required for observa- 
tion (13.0 km atm.), considerably higher dispersion than has been 
heretofore used might be necessary to detect it. Failure to observe 
these lines, even with sufficient dispersion, would at least give an 
upper limit for the amount of molecular hydrogen in the atmos- 
pheres of the major planets or low-temperature stars. 

The 4 — oand 5 — o bands would also be in a favorable position 
to be observed, but they would be much less intense. 


IV. THE QUADRUPOLE ROTATION-VIBRATION SPECTRUM OF J, 


Taking the data for the rotational constants of NV, in the ground 
electronic state from Watson and Koontz,” and for the vibrational 
constants from Birge and Hopfield,’ the wave numbers of the lines 
of the N, rotation-vibration bands have been derived. Table 2 gives 
the data for the zero lines of the bands only, which form the heads of 
the Q branches. In consequence of the smaller B and a values for NV, 
than for H,, the lines of the Q branches lie fairly close together, so 
that, except for very large dispersion, they will form an unresolved 


10 W. V. Norris and H. J. Unger, ibid., 43, 467, 1933, and unpublished results of the 
author. 

™ M.N., 93, 657, 1933. 

12 Phys. Rev., 46, 32, 1934. 13 Ap. J., 68, 257, 1928. 
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broad line or band. This coalescence of most of the lines of the Q 
branch will increase to some extent the sensitivity of the test for V,. 
The S and O branches lie to either side of the Q branch, with a spac- 
ing of roughly 8 cm“. 

A reliable estimate for the intensity of the V, quadrupole rotation- 
vibration spectrum is thus far not possible. One may, however, 
expect that the intensity will be at least of the same order of magni- 
tude as for H, but that it will drop more rapidly in the series of har- 
monics than for H,, since the anharmonicity is smaller. 


TABLE 2 


PREDICTED ZERO LINES OF QUADRUPOLE 
ROTATION-VIBRATION BANDS 
OF NITROGEN 


Band v(cm*) 
2330.7 42894 


If the intensity of absorption is large enough, the NV, rotation- 
vibration bands are, of course, to be expected in the solar spectrum, 
owing to the NV, content of the earth’s atmosphere, which corre- 
sponds to about 5 km atm. Only if the V, content of the planetary 
or stellar atmospheres is comparable to, or larger than, that of the 
earth’s atmosphere will it be possible to detect NV, in them. 

The frequencies of the 1 — o and 3 — o bands of N, coincide very 
closely with the intense fundamental v, of CO, and its second over- 
tone; the 3 — o band also coincides with the intense H,O band at 
1.4 u. Even with very high resolution, attainable only for the solar 
spectrum, it would be very difficult to detect the extremely weak NV, 
lines so intermixed with the intense CO, or H,0 lines. 

Fortunately, the 2 — oand 4 — oCO, bands are strictly forbidden 
for dipole radiation, and therefore the 2 — o and 4 — o N, bands 
will not be disturbed by CO, overlapping. They will also be fairly 
free from HO overlapping. In fact, it seems possible that a maxi- 
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mum at 2.16 win Langley and Abbot’s" curve of the infrared solar 
spectrum represents the Q branch of the 2 — o N, band. If this 
identification should prove to be true, it would represent the first 
observed example of a quadrupole rotation-vibration spectrum; and 
at the same time it would show the feasibility of the detection of 
N, and H, in planetary and stellar atmospheres by their quadrupole 
rotation-vibration spectra. 

The 4 — o band of N, lies in the photographic infrared and would 
therefore be much more easily accessible, especially in planetary and 
stellar spectra. But, in order to be observed, it will probably require 
much thicker absorbing layers than the corresponding H/, band, in 
consequence of the effect of the smaller anharmonicity mentioned 
above. Also, it is not certain whether in planetary spectra it will not 
be covered by CH, absorption, even though in the laboratory no 
CH, absorption has been observed in the region of the S and Q 
branches of the NV, 4 — o band with absorbing layers of 16 meters. 


V. DETECTION BY ORDINARY ROTATION-VIBRATION 
SPECTRA OF THE ISOTOPIC MOLECULES 

In addition to detecting H, and NV, by their quadrupole rotation- 
vibration spectrum, the following possibility exists: For the isotopic 
molecules, HD and N'N', respectively, the property “symmetric 
or antisymmetric in the nuclei” no longer exists. Consequently, an 
ordinary rotation-vibration spectrum (consisting of P and R 
branches only) can occur, which, of course, will be very weak because 
of the smallness of the change of dipole moment involved in the 
oscillations. Yet, the intensity of this weak dipole absorption would 
still be considerably larger than that of the corresponding quadru- 
pole absorption. Therefore, it does not seem impossible that the 
small concentration of HD and N‘*N*> molecules in natural hydro- 
gen or nitrogen could be compensated. 

Wick's has calculated the transition moment for the 1 — o band 
of HD. From his figure and the corresponding figure for HC/ taken 
from J. L. Dunham’s work," an intensity ratio of 1:6.6-10° for the 

4 Ann. Astrophys. Obs. Smiths. Institution, 1, 1900. 

1s Accad. Lincei Altti., 21, 708, 1935. 


© Phys. Rev. 34, 438, 1929. 
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1 — o bands of HD and HCl is obtained. Taking the abundance 
ratio of HD in H, equal to 1: 3000, the intensity ratio of the HD 
band in natural hydrogen to the HC/ band follows as 0.5 + 10~°, and 
is thus of the same order as, though somewhat smaller than, the 
figure for quadrupole radiation. The 3 — o and 4 — o bands of HD 
would lie at 9647 A and 7421 A, respectively. 

For N'4N*S the abundance ratio is much more favorable, but the 
mass ratio is much less favorable. The bands would lie very close to 
the N'*NV"4 bands given in Table 2. 


In conclusion I wish to express my sincere thanks to Dr. H. N. 
Russell for his keen and stimulating interest in this work. 
DEPARTMENT OF PHYSICS 


UNIVERSITY OF SASKATCHEWAN 
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QUADRUPOLE ROTATION-VIBRATION 
SPECTRUM OF H, 


HUBERT M. JAMES AND ALBERT SPRAGUE COOLIDGE 


ABSTRACT 


The theory of the quadrupole radiation of homonuclear diatomic molecules is de- 
veloped for the case of transitions between 'Y states. The results are applied in a 
theoretical treatment of the rotation-vibration spectrum of H.2, which, as Herzberg has 
suggested, may be observable in the absorption spectra of planetary atmospheres. The 
quadrupole moment of the system is computed as a function of the nuclear separation 
by the use of electronic wave functions obtained by a variational process. Matrix ele- 
ments of this quantity are calculated, using the vibrational functions of the Morse and 
Péschl-Teller potential-curves. Absolute values are found for the Einstein absorption 
coefficient B for fundamental, first, and second overtone bands, and are compared with 
the corresponding quantities for the dipole spectrum of HCI. It is found that the first 
and second overtone bands have intensities in the ratios 0.91 and 0.19, respectively, 
to that of the fundamental, and that, while this is an exceptional case, one should, in 
general, expect the ratios to be higher in quadrupole than in dipole bands. 


INTRODUCTION 


It has been suggested by Herzberg’ that it may be possible to de- 
tect the quadrupole rotation-vibration spectrum of H, in the at- 
mosphere of planets or cool stars, where the great depth of the ab- 
sorbing material can compensate for the small absorption coefficient. 
For the detailed discussion of this question, as found in Herzberg’s 
paper, it is necessary to Have a reasonably accurate quantitative idea 
of the magnitude of the absorption coefficients, both for the funda- 
mental band and for the lower overtone bands. Heretofore only 
crude estimates of the intensity of the fundamental band have been 
available,? the possible error being of the order of a factor of 10 
either way. No attempt has been made to discuss the relative 
strength of the overtones. 

The material of the present paper falls into two main parts. In 
the first section we apply the general theory of quadrupole radiation 
to transitions between "> states, obtaining transition probabilities 
in terms of the rotational quantum numbers and integrals over the 


Preceding paper. 
2 E. Teller, Hand- und Jahrbuch der chemischen Physik, 9, Part II, 52. 
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electronic and vibrational wave functions of the molecules. In the 
succeeding sections these integrations are carried out for the special 
case of transitions between the vibration-rotation levels of the 
ground state of H,, and absolute values are obtained for the absorp- 
tion intensity in the fundamental band and the first two overtones. 


I. GENERAL THEORY OF QUADRUPOLE TRANSITIONS 
BETWEEN STATES 


It is well known: that, corresponding to the expression 


A(a, b) = 


for the Einstein coefficient for spontaneous dipole radiative transi- 
tions of a system from state a to state 6, one has for spontaneous 
quadrupole radiative transitions 


A(a, 8) = (1) 


Here the Einstein coefficient A (a, b) is defined through the relation 
I,(a, b) = hv - N(a) A(a, (2) 


T,(a, 6) being the intensity of radiation arising from spontaneous 
transitions from a to b, and N(a) the number of systems in the initial 
state a. Nis a dyadic, the instantaneous quadrupole moment of the 
system 


N= (3) 


a 


the summation being taken over all charges; and (a/%/b), = 
(a/N,/b) is the matrix element between a and 6 of Jt., that compo- 


3 See, e.g., Condon and Shortley, Theory of Atomic Spectra (Cambridge; 1935), 
pp. 93 ff. We follow their notation in the main. 
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nent of 9 which is orthogonal to the unit dyadic 3. If we indicate the 
scalar of the dyadic by { },, then 


(a/N/b)2 = f — (4) 


The square of this dyadic indicates the double-dot product of this 
quantity with its complex conjugate; in the usual notation 


| = {(a/N/b)z (a/N/d)z}. (5) 


In general, in the absence of external fields, the energy-levels of 
the system will have a degeneracy of 2/ + 1, associated with the 
different possible orientations of the total angular momentum, and 
the observed spectral line will arise from the transitions from the 
2J 4 + 1 states of the upper level to the 2/z + 1 states of the lower 
level. For the case of natural excitation of the upper level, with the 
same number of systems, \V(a;), in each of the degenerate states, one 
can define 


1,(A, B) = hv N(a,) > A(ai, , (6) 


where A and B represent the totality of degenerate states a; and 
b;, respectively. Defining the Einstein coefficient for transitions be- 
tween the levels A and B by 


I,(A, B) = hv + N(A) - A(A, B), (7) 


one obtains 


(8) 


A(A, B) = 
2Jat 


Finally, defining S(A, B), the strength of the line, by 


S(A, B) = (9) 


| 
| 
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one has 
A(A, B) = (10) 
For absorption, the coefficient B(B, A) is defined by 
A) = hv - N(B) B(B, A) (11) 


1(B, A) being the rate of absorption of energy, and p(c) the density 
of isotropic, unpolarized radiation per unit difference of wave num- 
ber. Between A(A, B) and B(A, B) there exists, of course, the same 
relation as for dipole transitions, and thus 


4nv? S(A,B 
B(B, A) = (12) 

The problem to be considered is the evaluation of S(A, B) for 
transitions between 'D states of a molecule. For these states the 
coupling of the spin to the other co-ordinates of the system is 
negligible; and the quantum number for total angular momentum, 
J, is identical with K, the quantum number for rotation of the nuclei. 
Separating variables in the usual way, we may then write the wave 
function for the molecule as 


Wa = J, == F.(r; Vr, eee Revs (1) Pru (8, ¢) (13) 


Here F, is the electronic wave function, a function of 7,,72,...., 
the co-ordinates of the electrons in a system of axes fixed in the 
molecule, which involves also as a parameter the separation of the 
nuclei, 7. r'Rav(r) is the vibrational function, depending on the 
vibrational quantum number 2, on J, and on the electronic quantum 
numbers e, which determine the effective potential curve in the 
absence of rotation. Py. describes the rotation of the molecule with 
respect to axes fixed in space; in the present case of © states it is 


| 
| 
a 
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simply the normalized surface harmonic, which also describes the 
rotational motion of a single particle in a central field. We take 


Pan = (-1)™ P\"\(cos ; 


(14) 
Then 
(a/M/6). = f if X 


For simplicity we shall deal with the case of H,; the generaliza- 
tion to any other diatomic molecule is immediate and obvious. 


=> lc 


Then, in terms of J, J, K, a set of unit vectors fixed in the mole- 


cule with K along the internuclear axis, we have, with respect to 
the center of gravity, 


Ne = D ded (aX? — V3 — + (2¥} ZIT 
4 (222 — x? — Y)KK + + JD 
+ 3Vi2(JK + KJ) + 3X.Z(IK + KD]. (16) 


For either nucleus we have, independent of the electronic co-ordi- 
nates, X? = Y? = 0; Z? = r*/4. Integrating over the electronic co- 
ordinates, and denoting the resulting averages by a bar, we have, for 
either electron, X? = Y? and X;:V; = YiZ; = XiZ; = 0, because of 
the axial symmetry of the wave functions. The integral over elec- 
tronic co-ordinates which appears in equation (15) is thus 


= N(r)M(6, ¢) . (17) 


+ 
i 
{ 
i 
| 
i 
| 
} j 
} 
| 
i | 
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The first factor, V(r), upon averaging over the vibrational func- 
tions, gives V = N(r) (r)dr. Postponing the discussion 


of this quantity, we consider the averaging of the dyadic IN(8, ¢) 
over the rotational functions. We transform to a stationary frame of 


reference, with fixed unit vectors 7, 7, k, where 


= —singi+cos¢g, 


— cos 6 cos gi — cos sing] + sindk , (18) 


AY NY 


sin @ cos “i + sin 6 sin a cos Ok . 
We can conveniently express the result in terms of the dyadics 
>> >> 
= — + RR}, 
>_> >>> > 
R(+1) = ik + + (19) 


=> 


which are normalized and orthogonal to each other (and to 3). 
Then we easily find 


¢) = V 2P°(cos O)R(o) + 4Pi(cos R(+ 1) 
+ 4Pi(cos 1) + §P2(cos 2) 


+ §P32(cos 2) 
r=+2 


Hence 
T=+2 


Str = 28(7, | (21) 


T=—2 


| 


444. HUBERT M. JAMES AND ALBERT SPRAGUE COOLIDGE 


To evaluate fe ut P,,, Py,m dQ we express P,,, Py, as a sum 


of normalized surface harmonics. One easily finds the necessary 
formulae by iteration of 


cos 6Py,u = \ (2J + 1)(2J + 3) 


M 


(J — M)\(J + M) 
(2J — 1)(2J + 1) 


(2J + 1)(2J + 3) 


sin 


(J — M)\(J — M —1) 
(2J — 1)(2J + 1) 


Py41,M-1 


(2J — 1)(2J +1) 


Using the normalization and ‘orthogonality properties of the Py, x, 
one then gets immediately relations such as 


(Py, = 


(2J + 1)(2J + 3)?(2J + 5) 


J(J +1) — 3M 


J + 2) M) + 2 


(2J + 1)(2J — 1)?(2J — 3) 


3 


J 2) (23) 


i 
| 
| 
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From these formulae one obtains the desired values for (a/2/b),, as 


follows: 
AJ = J'—J =2, 4M =0: 
\ (2J + 1)(2J + 3)?(2J + 5) 
J(J +1) — 3M? — 
AJ = M= VR(o) ; 
AJ = — 2, AM =0: 
V2 N (oJ — 3)(2J — + 1) 
AJ = 2,AM = +1: 
AJ =0,AM=+4+1: 
+ N&K(+ 1) (24) 
AJ =—2,4M= +1: 
[Y= MJ + MU FM FM 
(oF = — +1) 
AJ = 2, 4M = +t 2: 
2 (2J + 1)(2J + 3)%(2J + 5) | 
AJ = 0, AM = +t 2: 
OF = Nes +3) 
AJ = —2, 4AM = +2: 
1 FM FM —2)I FM 
For all other AJ, AM (a/N/b), is zero. a 
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The strength of the transition between any two states, (a/3/b),?, 
is simply the square of the corresponding coefficient of & above, since 
the &’s are themselves normalized to unity. 

The strength of a line (Eq. [9]) is easily obtained by noting that 
>: (a;/N/b;).? must be independent of j; the probability of transi- 


tion out of a state is independent of the quantum number M, so 
long as there is no external field to define a preferred direction in 


space. Then 


S(A, B) = + 1) | 
= (2F4 1) (25) 


Thus, one can arbitrarily choose a; with M = 0, say, and easily 
carry out the indicated summations. One finds 


(FJ + + 2) 
N? X AJ = 2, 
2J (J + 1)(2J + 1) 
S(A, B) = — +3) (26) 
AJ = —2, 
2J —1 
and finally 
+ + 2) 
(2J + 1)(2J + 3)’ liad 
J(J — 1) 
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It will be noted that equations (24) and (26) agree with Condon 
and Shortley’s equations (6a) and (7), if one takes as the values of 
their D, E, F, 


—N 
E=o0, 
F 


(2J — 1)[(2J 3)(2J + ° 


2. QUADRUPOLE MOMENT OF THE GROUND STATE OF H, 


For the discussion of the vibration-rotation spectrum of H, we 
are concerned simply with the quadrupole moment of H, in the 
ground state; in the notation of the preceding section we have e = e’. 
We shall first consider that part of the quadrupole moment that 
gives rise to radiation as an average over the electronic co-ordinates, 
depending on the nuclear separation r as a parameter: 


A later section will be devoted to finding the averages, NV, of this 
quantity over the various radial functions. 
As a first approximation we have evaluated N(r) on the basis of 


the Heitler-London wave function 


where, in the usual notation, 7:0, 7», etc., are measured in Bohr radii 
(az), and S is the overlap integral 


I 
S= ‘f bdr, . 
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The integrations involved in computing V(r) are elementary, and we 
obtain 


+ + S(e~"/60)(60 + 60r + 277? + 773 + r)} 
Sm +r+%); 

S? 

N(r) = ger? 


Values of N(r) for the range of r which is of interest appear in Table 
3. From r = tay to r = 2a, the function is represented to within 


I per cent by 
N(r) = e{o.285 + 0. 263(r — 1.4) — 0.084(r — 1.4)?}. 


The effective quadrupole moment thus appears as a positive 
definite quantity, which arises as the difference between a positive 
contribution er?/4 from each proton and a negative definite contribu- 
tion —e(Z? — X?) from each electron. The cancellation is fairly 
severe, the whole being of the order of one-fourth of the positive and 
negative components separately. This is but one element of the 
situation that makes it impossible to accept the foregoing results as 
reliable. The Heitler-London function is known to be a rather poor 
approximation to the correct function for this molecule, and in par- 
ticular it gives too great a probability density of electrons far from the 
nuclei and also near the nuclear axis beyond the nuclei. Thus, the 
function will be expected to give too large a value for Z? — X?, and 
accordingly too small a value of V(r). Furthermore, it is the rate of 
change of V(r), and the derivatives of this quantity, which finally 
determine the intensity of radiation; since the Heitler-London func- 
tion becomes worse in a systematic way as the nuclei approach, one 
will expect a systematic change in the error of V(r), and consequent- 
ly a large error in the derivatives. For these reasons we have con- 


4H. M. James and A. S. Coolidge, J. Chem. Phys., 1, 825, 1933. Cf. also H. M. 
James, ibid., 3, 9, 1935. 
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sidered it necessary to compute \V(r) by a more accurate method. 
We present the Heitler-London results merely for comparison. 

For the more accurate computation we have employed the method 
which we have previously? devised for a careful study of the ground 
state of H,. The function is of the form 


F= Crnikplmnjk p] 
mnjkp 
where [mnjkp] stands for (ATAZ + 
with 


_ tra + Via — Vib 27 12 


The manipulation of these functions is fully described in the paper 
referred to, and the coefficients C,,n;,, are there given for some quite 


TABLE 1 
COEFFICIENTS IN WAVE FUNCTIONS 
6=0.75 6=0.875 
TERMS 
r=1.1dy T=1.3 r=1.5 r=1.5 7=1.7 r=1.9 
foeode)......... 2.61623 2.465097 2.03888 | 2.37301 2.10568 1.76079 
[00020]... 0.86493 0.75604 0.64527 I. 23006 1.07702 0.93081 
foorto]...... —0.22738 |—0. 26338 |—o. 26011 |—0.42616 |—o.42I101 |—0. 39320 
[10000]... ... —0.30175 |—0.59461 |—0.59251 |—0.59444 |—0.66606 |—0.61730 
[oooor]...... 0.41034 0.25671 0.16236 ©. 33075 0.22241 0.15108 


accurate functions, for values of r ranging from 1.2 ay to 1.7 dy. 
For the present purpose it was sufficient to use simpler, though less 
accurate, functions, such as the five-term function containing the 
terms which we had already found were sufficient, at r = 1.4, to give 
a very fair approximation to the correct energy value. Functions 
were accordingly worked out for five different values of 7, using two 
different choices of 6. Table 1 shows the coefficients defining the 
normalized functions, and Table 2 gives the corresponding calculated 
binding energy, compared with a Morse curve fitted to the experi- 


| 
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mental data of Beutler and Jiinger’ and of Jeppesen.°® It will appear 
from these results that, for a given r, the best value of 6 is somewhat 
greater for the five-term function than it is when more terms are 
employed. The two independent computations at r = 1.5 make pos- 
sible a check on the effect upon V(r) of errors in the wave function. 


TABLE 2 
COMPUTED BINDING ENERGY, IN ENERGY VALUE 
6=0.75 6=0.875 
r=1.1 r=1.3 r=1.5 r=1.5 r=1.7 
Variational function........ 3.95 4.40 4.37 4.10 3.61 
| 


N(r) may be evaluated in terms of the set of quantities employed 
in computing the energy of the functions. Since Xj; = Yj = X}, we 
may write 


NO) = + ARTF VD + TD Zi 


Now, in elliptical co-ordinates we have | 


Xi + Vi = — 3)(1 — Zi = 


Since all the terms in the wave function involve the electrons 
symmetrically, we have relations such as (Ar + Az)[mnajkp| = 
[(m + 2) njkp| + [m(n + 2)jkp], and therefore | 


N(r) = Connie (S(m + 2, ky 


mnjkp 
+ S(m,n + 2,j,k, p; m’,n', 7’, k’, p’) 
+ S(m, n,j + 2,k, p; m,n’, 7’, k’, p’) 
+ S(m, n, j,k + 2, m’, n’, 7’, k’, p’) 
— 35(m + 2,0, j + 2,k, p; n’, 7’, k’, p’) 
— 35(m, + 2,7, 8 + 2, p; m’, 
+ 2S(m, n, j,k, p; m'’, n’, 7’, k’, 
5Zs.f. Phys., 101, 304, 1930. 6 Phys. Rev., 44, 165, 1933. 
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S is here the integral n,j,k, p\[m’,n’, 7’, k’, p'| which may 
be evaluated by the methods given in reference 4. 


The results of these computations appear in Table 3. They can be 
approximated by 


N(r) = ef0.442 + 0.307(r — 1.4) — 0.205(r — 1.4)?}, 


the values of which are also given in Table 3. It is reassuring that no 
serious discrepancy appears between the two calculations at r = 1.5. 
Since the values of 6 used were, respectively, considerably smaller 


TABLE 3 


EFFECTIVE QUADRUPOLE MOMENT XN(r), IN PROTON 
CHARGE X (BOHR RADIUS)? 


r=I1.1dyq F=1.9 rT=4.5 
Computed, H—L function. .198 310 .358 .395 | 0.147 
Quadratic approximation. . 409 .470 .544 |— .576 
Cubic approximation...... 0.331 0.409 | 0.470] 0.515 | 0.544 | 0.178 


and (probably) somewhat larger than the most advantageous choice, 
the errors in the functions are presumably of opposite character, 
particularly with respect to the scale of the electronic distribution as 
a whole. An intermediate value for V(r), such as that given by the 
quadratic approximation, is therefore probably quite near the truth. 
We estimate the error as of the order of 3 per cent in V(r), as given 
by that expression, 5 per cent in dN /dr, and 20 per cent in d?N/dr’. 
It will be noted that, as expected, the value computed with the 
Heitler-London function is too low, though the rate of change is 
about right. On the other hand, the second derivative, which is 
particularly interesting for our present purpose, comes out too small 
by more than half. 

In view of the importance of the second derivative, it is well to 
have some check on the order of magnitude of the third. This can 
be obtained by considering the situation for large values of r, where 


| 
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the Heitler-London function is less in error. Thus, the cubic approxi- 
mation 


N(r) = e{0.442 + 0.303(r — 1.4) — 0.212(r — 1.4)? 
+ 0.028(r — 1.4)3} 


is indistinguishable from the quadratic expression over the range of 
interest but is much better at large 7, as will be seen from Table 3. 


3. EVALUATION OF THE RADIAL INTEGRALS 


In order to discuss the intensities of the fundamental rotation- 
vibration band, as well as the first two overtone bands, we have now 
to consider the evaluation of the radial integrals 


M(e, », J; v', J’) = Rost 


for v = 0, v’ = 1, 2, and 3. These integrals will depend principally 
on v’, and also to a minor extent on the values of J and J’. Following 
Gibson, Rice, and Bayliss,’ we shall neglect the dependence on rota- 
tional quantum numbers, though it should be noted that, in the 
present problem, the values of the integrals are probably more sensi- 
tive to changes in the effective potential curve than in the case which! 
they discuss. We shall, then, use radial functions R,.j, Ryyy, com- 
puted forJ =J’=0. 

It is instructive to regard N(r) as expanded in powers of the dis- 
placement from the potential minimum r,, and to investigate the 
components of V produced by the separate terms. Thus, if 


we have 
Now = NG — + ANG — 


For any potential curve, the contribution from the constant term 
N. is, of course, zero, because different wave functions are orthog- 
onal. The contributions from the higher terms are quite markedly 


7 [bid., 44, 193, 1933- 


| 
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dependent upon the exact form of the potential curve. The simplest 
possible approximation is a parabola adjusted to give the proper 
values of the spectroscopic constants 7, and w,, and it is a property 
of the familiar wave functions of the harmonic oscillator that 
(r — r.)® » vanishes unless m — v’ is even and positive. Thus, to the 
approximation in which the third and higher derivatives of V can 
be neglected, the strength of the fundamental band would be de- 
termined entirely by dN /dr evaluated at r,, the first overtone band by 
d’?N/dr’, with other overtones absent. Explicitly, 


(7 — = — , 
= aly, 

(= = — , 
(7 — = — , 


with R the Rydberg constant, m the mass of the electron, and u the 
reduced mass of the nuclei. In this and in ensuing expressions, the 
signs correspond to wave functions chosen so as to have the first loop 
for small r positive. 

When a better approximation to the true potential curve is used, 
its departure from the parabolic form results in contributions from 
each derivative for every value of v’. For the Morse curve, general 


expressions for (r — r,)%, may readily be found by extending the 


e/o,v 


treatment of Rosenthal.’ We introduce the functions 


with = 1, Fi(x) = p(x), Fix) = (x) + F(x) = 
(x) + (x) (x) + By repeated differentiation of the 
relation I'(x) = (w — 1)I'(x — 1), we readily obtain 


(v — 1)[Fa(x) — — 1)] = — 1). 


5 Proc. Nat. Acad., 21, 281, 1935. 


| 
| 

| 

| (x) 

= 
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The final formulae are found to be 


(r I = »| 2) k)/a)(r — I 
— — 1)((mk)/a)(r — 


+3.2.(n —1)(n — 2) ((Ink)/as(r 


+ (# — — 3)] + (ek) — 


Hn = (nk) = + 


+ ((nk)/a)(r — — — 1) (nk) — 


Here, a and & are a reciprocal length and a number characteristic of 
the Morse potential function. When the curve is adjusted to give 
the correct value of the spectroscopic constant w,x, expressed in 
wave numbers, we have 

We 


a = ; k= 
WeN'e 


Using the spectrocopic constants given by Jeppesen,° we obtain the 
numerical values given in Table 4, where they may be compared with 
the results given by the parabolic approximation. It will be noticed 
that the terms arising from the nonparabolic form are by no means 
unimportant. 

Table 4 also gives the values of V(r) averaged over both para- 
bolic and Morse curve functions and computed both from the quad- 
ratic and the cubic approximations. For each potential function, 
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the effect of including the cubic term is quite unimportant; and it is 
clear that, so far as the first three transitions are concerned, we have 
calculated V(r) over a sufficient range to determine the results with 
fair accuracy. On the other hand, in going from the parabola to the 
Morse curve, we encounter major changes in the ratios of the N’s, 
and still greater changes in the ratios of the predicted band intensi- 


TABLE 4 


VALUES OF INTEGRALS OVER WAVE FUNCTIONS BELONGING 
TO VARIOUS APPROXIMATE POTENTIAL CURVES 


Parabola Morse 
—0o.1644 —o.1668 
° — .0044 
Quad. Cubic Quad Cubic Quad. 
Appr. Appr. Appr Appr Appr. 
No, —0.0505 | —0.0502 | —0.0462 | —0.0459 | —0.0441 
rr —o0.0078 | — .0o081 | — .0137 | — .0136 | — .0166 
° —0,0003 | —0.0037 | —0.0039 | —0.0044 


ties, which are proportional to the squares of the N’s. Analysis of 
the situation shows that the differences are due mainly to the asym- 
metry of the Morse curve, which tends to increase the relative in- 
tensities of the overtone bands. In view of the sensitiveness of the 
results to the asymmetry of the potential function, and the fact that 
the Morse curve is in this respect by no means a good approximation 
to the true curve, it seemed important to examine the results ob- 
tained with some more accurate potential function. 

The most direct indication of the asymmetry of any potential 


| 
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curve is the calculated value of the spectroscopic constant a,, which 
determines the principal term in the energy due to the interaction 
of rotation and vibration, and is proportional to the rate of change 
of the average moment of inertia of the molecule with the vibra- 
tional quantum number. According to Jeppesen,° a, has for H, the 
unusually large value 3.07; whereas the parabola gives a, = 
— 6B?/w,. = — 5.03, and the Morse curve, as here constructed, gives, 
according to Pekeris,? a, = 6x,B,[(B./x.w.)'? — B./x.w,.| = 2.36. 
Now, the curve given by Péschl and Teller’ contains four adjustable 
parameters, one more than the Morse curve. It can therefore be ad- 
justed to give the correct value of a,, as well as B,, w,, and x,w,, 
provided one does not require it to give the proper dissociation 
energy; this proviso (which also applies to the Morse curve used) 
is legitimate in the present case, since it affects only the more remote 
parts of the curve, which are practically without influence on the first 
few wave functions. We have worked out such a curve, but postpone 
for a later communication a discussion of the method. The equation 
for the curve is 


he sinh?a(r —  cosh?a(r — ro) |’ 

with a = 0.52414@y ', fo = —0.0696ay4, = 24.773, = 57-436. 
(The yw in the first factor is; of course, the reduced mass of the mole- 
cule, as before.) 

It is possible to obtain analytical expressions for (r — r,)". ,,, using 
the wave functions of the Péschl-Teller curve; but the results appear 
as a series converging too slowly for practical computation. We 
have been fortunate in being able to use the differential analyzer of 
the Massachusetts Institute of Technology in evaluating the in- 
tegrals V,, y. The technique of this work need not be described here, 
as it was very similar to that previously discussed in connection 
with another investigation." Unfortunately, the curve employed 


9 Phys. Rev., 45, 98, 1934. 
10 Zs. f. Phys., 83, 143, 1933- 
™ Coolidge, James, and R. D. Present, J. Chem. Phys., 4, 193, 1930. 


| 
| 
| 


THE SPECTRUM OF H; 457 


was not exactly that just given, but one based upon Lotmar’s” 
formulae for constructing a Péschl-Teller curve so as to reproduce 
given spectroscopic constants. These formulae contain an error 
which was not detected until after the work had been completed, so 
that the results obtained correspond to a potential curve giving 
spurious spectroscopic constants, slightly different from the true 
values. We have made an estimate of the corresponding errors in 
the V’s by comparing the values computed with Morse curves so 
constructed as to give, respectively, the spurious and the true 
spectroscopic constants. The appropriate corrections (amounting to 
5 per cent for V,,, and less for V,,, and N,, ;) have been applied to 
the values given in Table 4. We believe these may be relied upon to 
I or 2 per cent, except for V,, ,, which may be in error by as much as 
5 per cent, since the extensive cancellation between contributions 
to the integral from alternating positive and negative parts of the 
wave function forv = 3 makes the final result very sensitive to errors 
in the form of the potential curve used and in the operation of the 
machine. It is upon these values of V that the subsequent discus- 
sion will be based. It is to be regretted that we did not realize the 
desirability of obtaining the separate values of the (r — r,)? , for 
comparison with the parabola and Morse results. Greater accuracy 
(as far as machine operation is concerned) could have been secured 
thus, by eliminating the contribution from ” = o at the outset in- 
stead of relying upon the machine to cancel the corresponding term 
in the integral. 


4. RESULTS AND DISCUSSION 


Evaluating the constants of equation (27), using the results of the 
previous section, we have 


(J + + 2) = 
(2J+ 1)(2J + 3)’ 
J(J + 1) 
BB, 1) Cy, wo x 3(2J — 1)(2J + 3)’ AJ oO, 
J(J — 1) 
(2J — 1)(2J +1)’ 


2 Zs. f. Phys., 93, 528, 1935. 
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with 
Co,2 = 0.925 X 10°", 
Co, 3 


Il 


0.067 X Io". 


Special values of B(B, A) for the first line of the S branch in each 
band are as follows: 


4.35X10-, 
Second overtone.................. 0.65X10~. 


With these we may compare the value B = 1.71 X 10° given by 
Bourgin’ for the first line of the R branch of the fundamental band 
of HCl. (The number given by Bourgin is actually greater by a 
factor of c, since he defines B(B, A) in terms of energy density per 
unit difference of frequency rather than of wave number.) 

As the proper measure of the intensity of an absorption line, we 


take the integrated absorption coefficient } a(a)do. It is well known" 


that if one neglects the width of the line, which is, of course, extreme- 
ly small for the very improbable transitions which we are consider- 
ing, then 


a(o)do = (hv/c)B(B, A)N(B) . 


For a given density of absorbing atoms, the integrated absorption 
coefficients for the H, lines thus stand in the ratios 8.1 X 10 °, 
7.3 X 10 °, and 1.5 X 10 °, respectively, to that for the HC? line. 

The ratios of the intensities of the first and second overtones to 
the fundamental are 0.91 and 0.19, respectively. These large values 
are in striking contrast to the corresponding values, 0.016 and 
~ 0.001, for HCI.'5 This fact can be accounted for as follows: The 


expression for f alo)de is proportional to v’ for quadrupole radiation 
and to v for dipole radiation, which increases the ratio of the first 


13 Phys. Rev., 29, 794, 1927. 
14 Condon and Shortley, op. cil., pp. 109 ff. 
18 J. L. Dunham, Phys. Rev., 34, 438, 1929. 
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overtone to the fundamental by a factor of 4. A further factor of 4 
in this ratio comes from the fact that in our case (to judge by the 
Morse curve results) Ni(r — 7.)o,, and 4N%(r — r,)2,, make com- 
parable contributions to V,,,, whereas Rosenthal! finds that for HC1 
the term analogous to the second is probably negligible because of 
the smallness of P? (the value at the potential minimum of the 
second derivative of the electric moment function P{r]). The re- 
maining factor of 3 is due to the exceptional asymmetry of the poten- 
tial curve for H,, for which a/B, = 0.052, as against 0.029 for HCl. 
Of these three factors, the first is certainly, and the second possibly, 
a general characteristic of quadrupole radiation, so that it is to be 
expected that in quadrupole bands in general the overtones will have 
greater intensities, relative to the fundamentals, than in dipole 
bands, although not to the exceptional degree shown in H,,. 


In conclusion, we wish again to express our indebtedness to the 
Electrical Engineering Department of the Massachusetts Institute 
of Technology for the use of the differential analyzer, and to thank 
Dr. E. L. Vernon for assistance in its operation, as well as in some 
of the numerical work of the paper. In addition, one of us (H. M. J.), 
wishes to acknowledge gratefully a grant from the Society of Sigma 
Xi, which materially aided in the progress of the work. 
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ON THE DETERMINATION OF COLOR INDICES OF 
STARS FROM A CLASSIFICATION 
OF THEIR SPECTRA 


W. W. MORGAN 


ABSTRACT 

A discussion of Bottlinger’s photoelectric color indices of the near-by stars of types 
F5-—Mo is made in combination with a classification of spectral type and luminosity. 
It is found that, within the errors of observation, the intrinsic colors of the stars are 
uniquely defined by their spectral types and absolute magnitudes. 

The color indices of stars unaffected by selective interstellar absorption can be pre- 
dicted with accuracy from a classification depending on line ratios within their spectra. 
The suitability of the F5—-Go main sequence stars as color standards is pointed out. 
For these stars color indices can be predicted which have a mean error of the order of 
+0.02 mag. on Bottlinger’s scale. 

1. In an important paper published some years ago, Opik' in- 
vestigated the problem of the practical determination of luminosities 
of stars from their spectral types and color indices. Elaborating on 
the pioneer work of Adams, Seares, Lindblad, and others, he found 
a good correlation between absolute magnitude and color index for 
stars of the same spectral type, except in the case of the A and early 
F stars. For types Gand K, especially, there is a considerable increase 
in color index on passing from the dwarfs to the normal giants and on 
to the supergiants. 

Investigations of the nature and amount of reddening due to inter- 
stellar causes which have been made since Opik’s paper was pub- 
lished, have shown that in the case of B and O stars a considerable 
part of the reddening which is observed with increasing luminosity 
is due to interstellar material. If stars of any given apparent mag- 
nitude are considered, then the more luminous ones are also the more 
distant and will be more likely to be affected by interstellar redden- 
ing. A marked effect is to be expected particularly in the case of the 
B and O stars which possess the characteristics of high luminosity 
and a high degree of concentration toward the galactic plane. In the 
case of the G and K stars, however, the reddening observed with 
increasing luminosity is largely intrinsic, as it is well shown in ob- 


* Tartu Pub., 27, 1, 1929. 
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jects which are so near to the sun that the reddening cannot be ac- 
counted for by interstellar causes. 

It is the purpose of the present paper to reverse the problem treat- 
ed by Opik and to investigate the accuracy with which intrinsic 
color indices (that is, color indices which are unaffected by selective 
interstellar absorption) are defined by spectral type and luminosity 
for stars of types F5—Ks5. Instead of using explicit values of the 
absolute magnitudes of the stars discussed, however, an abstract 
luminosity classification based on spectroscopic criteria will be em- 
ployed. The advantages in the use of such a classification have been 
discussed in a recent paper.’ 

2. The stars suitable for discussion are limited by considerations 
of selective interstellar absorption. If observed color indices are to 
be considered essentially intrinsic in nature, then it is necessary to 
include only the near-by stars. A correction to the observed color 
index of more distant objects as a simple function of the distance is 
not reliable because of the marked irregularity in the distribution of 
the absorbing material; in addition, the characteristics of the inter- 
stellar matter are not, as yet, completely known. In a few regions of 
the northern sky interstellar absorption becomes appreciable at dis- 
tances which are probably somewhat less than one hundred parsecs. 
If we limit ourselves to stars whose apparent distances are less than 
one hundred parsecs, however, selective effects should, in general, be 
quite small because of the large ratio of selective to total absorption. 
A round figure of HR magnitude 4.00 has been set as the faint limit 
for ordinary giants of types F5-K5. The apparent distance of very 
few objects so chosen will be as great as 100 parsecs. A fainter limit 
has been set for the main-sequence stars; the ones included all lie 
within a distance of 50 parsecs. 

Because of the low space density and the great distances of the 
supergiants, it is necessary to include more distant objects if a repre- 
sentative selection is to be made. The limiting apparent visual mag- 
nitude of 4.00 has been adopted. This is the same limit which was 
set for the normal giants; as a result, the mean distance of the super- 
giants will be much greater, and the observed color indices are more 
likely to be affected by interstellar reddening. As some of the bright 


7Ap. J., 85, 380, 1937- 
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supergiants are located at apparent distances of around 500 parsecs, 
their colors may be appreciably affected by selective absorption. 
The majority of the F, G, and K supergiants brighter than the limit 
set are, however, located in regions which are relatively transparent 
up to distances of 300-400 parsecs. 

3. The observational material is of two kinds: accurate color in- 
dices and accurate classifications both of spectral type and luminos- 
ity. Their low accidental error, together with their freedom from 
several troublesome sources of systematic error, makes photoelectric 
color indices most suitable for discussion. The series of photoelectric 
color indices of Bottlinger’ will be used in the present discussion.‘ 
Bottlinger’s original measures have been used in preference to his 
reduction of the color indices to King’s scale because of the compli- 
cated nature of the relation between the two color systems. The 
scale of the photoelectric color indices is approximately one-half that 
of King’s; the average mean error of a single observation of color of 
highest weight is +0.016 mag. on Bottlinger’s scale. Observations 
of full weight only have been used. All stars known to have com- 
posite spectra have been omitted. 

In spite of the fact that all the stars used in the discussion are 
bright and are included in two or more series of spectral classifica- 
tion, there are reasons for re-examining the question of their classi- 
fication. In problems which make use of spectral types of faint stars 
the assumption is usually made that a star of the tenth magnitude, 
for example, which has been classified as Ko, is on the system of the 
bright Ko stars. This, in turn, is dependent on the use of the same 
spectral criteria for bright stars and for faint ones, and to a certain 
extent on the comparability of the spectral dispersions. In general, 
the bright stars are observed with rather high dispersion to obtain 
high accuracy in classification, while low dispersion is employed of 
necessity for faint objects. In the latter case the accidental errors 


3 Veriff. Berlin-Babelsberg, 3, Part IV, 1923. 

4 The carefully compiled catalogue of Bottlinger’s re-reduced color indices and the 
new determinations of W. Becker (Veréff. Berlin-Babelsberg, 10, Part 3, 1933) indicate 
the same conclusions found in the present paper from Bottlinger’s earlier published 
measures. The excellent infrared color indices of Hall (4p. J., 79, 145, 1934) have not 
been used because of their great difference in effective wave length as compared with 
ordinary photographic and visual systems. 
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of classification tend to be higher, and the connection to the system 
of the brighter stars is not as well defined as it would be if the fainter 
stars could be observed with higher dispersion. 

From these considerations certain general characteristics of an 
“optimum” classification for bright and faint stars may be described: 
(a) The dispersion should be low for photographic efficiency. (b) The 
resolving-power of the spectrograms should be high for accuracy in 
classification and for the determination of spectroscopic luminosities. 
(c) There should be no systematic difference in the quality of the 
spectrograms which depends on the apparent magnitude of the stars. 
The use of short-focus cameras on slit spectrographs satisfies all of 
the conditions. Short-focus lenses of the type designed by G. W. 
Moffitt, when used in connection with the modern high-speed, fine- 
grain emulsions, give excellent low-dispersion plates of high resolv- 
ing-power which are suitable for the determination both of spectral 
types and of spectroscopic absolute magnitudes. If care is taken that 
spectrograms of faint stars are not systematically underexposed in 
comparison with bright standard stars, the possibility of serious sys- 
tematic error as a function of apparent magnitude is eliminated. 
The great disadvantage of the slit spectrograph is, of course, that 
only one star can be observed at a time. For this reason, where very 
large numbers of stars are to be observed, the objective prism meth- 
od is the only one to be considered. 

The spectra of all the stars included in the present discussion 
have been obtained with a short-focus Moffitt camera attached to 
the one-prism Bruce spectrograph on the forty-inch telescope. Cri- 
teria of the greatest sensitivity to changes in spectral type were se- 
lected, and a system was set up which is based on the Draper types 
of stars brighter than visual magnitude 4.0. The system is in close 
agreement with the Draper system; there is a maximum correction 
of about two-tenths of a spectral type at G5. The procedure adopted 
at Mount Wilson of placing type Mo immediately after K5 in the 
case of the giants has been followed; spectroscopic evidence and that 
of the color indices agree in making the interval between K5 and 
Mo no larger than that between K4 and Ks. 

The criteria used for the classification of spectral type are listed 
in Table 1. 
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The spectral type at which the components of line ratios become 
equal probably varies considerably with the dispersion used. For 
this reason the system of classification should be considered as being 
defined by the whole catalogue of stars which have been classified, 
rather than by arbitrary measures of certain spectral criteria. 

In addition to the classification of spectral type, the stars have 
been divided into a number of luminosity classes by means of several 
well-known line ratios which have been employed for the determina- 


TABLE 1 
CRITERIA OF SPECTRAL TYPE 


Type Criteria Remarks 
4045/H6 As in a Persei 
er 4045/H6 4045=H6 at Go on spectrograms cf the dis- 
Hb > 4226 persion used 
4045 
Hé=4226* 
G4-Ko...... 4226 is interpolated be- 


tween 4045 and Hé 


Co ae 4045/4226 4045=4226 at Kr on spectrograms of the 
dispersion used 


TiO 6159 Red 7iO band suspected at Ks, faintly 
present at Mo 


* On spectrograms of the dispersion used. 


tion of absolute magnitudes. These ratios are listed in Table 2. The 
luminosity classes used are defined in Table 3. 

A catalogue of the spectral and the luminosity classification of 
114 Stars is given in Table 4. All stars of types F5—Mo brighter than 
4.00 in the Revised Harvard Photometry for which Bottlinger meas- 
ured color indices are included; in addition, a number of fainter 
main-sequence stars have been added. The columns of the table 
give: (1) the serial number; (2) the name of the star; (3) the posi- 
tion for 1900; (4) the HR visual magnitude and Draper type; (5) the 
Yerkes spectral type, with the number of spectrograms on which it 
was determined in parentheses; (6) the total range in the individual 
determinations of the spectral type; (7) the luminosity class; (8) the 
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photoelectric color index as measured by Bottlinger, with the num- 
ber of first-class observations on which it depends. 

The mean error of a determination of spectral type from one plate 
was computed by Schlesinger’s’ method for all stars for which there 
were between three and six plates. The accuracy varies considerably 
with spectral type; the highest precision is obtained for stars of 
types K1—Mo, where half-tenths were commonly estimated. Next in 
order of accuracy are the F4—Go stars, while the classification of 


TABLE 2 
LINE RATIOS USED FOR LUMINOSITY CLASSIFICATION 
Line Ratio Remarks 
7” Range Used 

Fs—Mo Fundamental criterion of luminosity over great 
range in spectral type 

F5-Go Especially useful in detecting F-type super- 
giants 

Fs5—-Mo Lines weak in F stars for dispersion used; espe- 
cially useful in separating G8—K1 subgiants 
from dwarfs 

Fs5—Mo Rather weak in F stars for dispersion used; use- 
ful in detecting supergiants of types Gs5s—K5 

Gs5~Mo Useful in detecting supergiants of types 
G5-Ks5 


Gi-Gog spectra has the lowest accuracy of all. The values of the 
mean error for individual spectrograms are given in units of tenths 
of a spectral class. 


F4-Go G1-Go Ko-Ks5 


Mean error of one plate........ +0.8 23:3 +0.4 


The mean catalogue values have been rounded off to the nearest 
half-tenth; systematic errors due to changes in the personality of 
the observer and the relatively indefinite nature of the definition 
of the classification system tend to make the external errors some- 


5 A.J., 46, 162, 1937. 
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TABLE 3 
LUMINOSITY CLASSES 
Class Description Type Stars Myis Source of M and Remarks 
Ia. Supergiants of F4 e Aur —7.2 | From Kuiper’s solution of orbit* 
types F4-F8 of | F7 p Cas 
highest luminos- | F7+ 6 CMat 
ity 
Ib.......| Supergiants of Fs a Per —4.4 | In Perseus cluster: cluster modulus =6™3§ 
types F4-F8 of | F7 a UMi —4.2 From c.p.m. companion, which is main-se- 
somewhat lower quence F2 star. Absolute magnitude of 
luminosity companion assumed to be +2.5. Mean 
error is equal to dispersion of F2 dwarfs, 
which is around +1 mag. 
| Soeur Supergiants of Go ¢ Mon 
types later than | G4 9 Peg 
F8 G7+ 
Kr+ 
K2 n Per —4.8 | From c.p.m. companion, which is main-se- 
quence Bo star. Absolute magnitude of 
companion is considered 0.0 with a mean 
error equal to the dispersion of Bg dwarfs, 
which is around +1 mag. 
K3+ &Cyg 
Ks5s+ aSco —4.et] From parallax of diffuse nebulosity in which 
Scorpius stars are involved 
(Mr) a Ori 
Stars intermediate | Gr Leo 
in luminosity be- | Go ¢ Cyg 
tween super- K2+ yAnd(A)| —1.5 | From dynamical parallax of y And(BC). B 
giants and ordi- is main-sequence star of type Ao. If abso- 
nary giants lute magnitude of +0.5 is assumed for B, 
then M of y And (A) is —2.6 
K2+ 
K2+ y¥Aql 
fDi Normal giants Fo e Hya +o.1 | From dynamical parallax 
G2 a Aur 0.0 | From spectroscopic and _ interferometric 
measures 
G4+ Peg 
G6 n Her 
Ko 8B Gem +1.2 | Trigonometric parallax =o%100+ 10 
a Boo —o.1 Trigonometric parallax=o.085+ 4 
K5 aeTau —o.6 | Trigonometric parallax=o0.046+ 4 
pAnd +o.6 | Trigonometric parallax =o.044+ 6 
Gs-K2 subgiants | Go n Ser +2.0 Trigonometric parallax =o0.053+ 7 
Kr y Cep +2.6 | Trigonometric parallax=o.067+ 5 
Main-sequence F4 Peg +3.4 | Trigonometric parallax=o.077+ 6 
stars F6+ = 73 Ori +3.8 | Trigonometric parallax=o.128+ 5 
F8 B Vir +3.8 | Trigonometric parallax=o.101+ 7 
Fo+ 1»Cas (A) | +4.9 | Trigonometric parallax=o.182+ 5 
Gi+ Tri +5.0]| Trigonometric parallax=o.098+ 6 
G8 p Cas +5.8 | Trigonometric parallax=o.130+ 6 
G8+ «Dra +6.1 | Trigonometric parallax=o.181+ 4 
K2 ¢ Eri +6.2 | Trigonometric parallax=o.305+ 7 
K3 HR 8832 +6.5 | Trigonometric parallax=o0.146+ 6 
Ks5+ 61CygA | +8.0] Trigonometric parallax=o.299+ 3 


* Ap. J., 86, 575, 1937. 
t Half-tenths of spectral classes are designated by plus signs (+). 
t Absolute magnitude may be affected by interstellar absorption. 
§ Lick Obs. Bull., 14, 170, 1930. 
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TABLE 4 
SPECTRAL TYPES AND COLOR INDICES OF STANDARD STARS 


Star a1goo 51900 | my Spy | Spy wiR Ly 
B Cas of 4™+58°6| 2.42 Fs | F2 (3) | 2| 
6 And © 34 30.3] 3.49 | K3 (3) }o} ITI 
a Cas © 35 56.0] 2.47 Ko | Kr (6) | « | 
nCasA | 043 57.3} 3-64 F8 | Fot (6) | 3] V 
Cas I 2  §4.4| 5.26 Gs | G8 (6) | 4 V 
And I 35-1| 2.37. Ma| Ks5+ (2)} 0] IT | + 
a UMi 123 88.8} 2.12R F8 | F7 (4) | 2 Ib 
n Psc 126 14.8] 3.72. Gs | G8+ (2)] 1 III 
a Tri 147 29.1) 3.58 Fs | Fs (3) | 1 | 
a Ari 2 2 23.0] 2.23 K2|{Ki+ (5)| 2) Il 
6 Tri 211 33.8} 5.07 Go | Gi+ (3) 1 V 
6 Per 237 48.81 4.22 F8 | F7+ (3) 1 V 
n Per 243 3-03 Ko | (4) | 1 I 
a Cet 2 57 3.7| 2.82 Mo ITT |+ 
«Per 3 2 49.2) 4.17 Go| Fot (4) ]1 V 
a Per 317 49.5, 1.90 F5 | F5 (5) | 1 Ib 
o Tau 3 19 8.7| 3.80 Gs | G8 (2)/2| I {- 
v Per 338 42.3}3-03 Fs | (3)/0] Ib | 
Tau 414 15.4| 3.86 Ko | Got (2)/3] IT | 
6 Tau 417 17.3] 3-93 Ko | Gg (2)|} 2] I 
Tau 423 19.0] 3.63 Ko | Ko (3)|o} I 
a Tau 430 16.3] 1.06 Ks5 | Ks (6)| 1} | + 
Ori 4 44 6.8] 3.31 F8 | Fo+ (3) 1 V 
«Aur 451 33.0] 2.900 | K2+ (3)] 1 II {+ 
e Aur 455 43-73-32" Fs Wir; i= 
a Aur 5 9 45.9] 0.21R Go | G2 (6)|3|] I |- 
Aur 512 40.0] 4.85 Go | Go (3) | 2 V = 
x? Ori 549 20.3] 4.62 F8 | Go (3) | 2 V = 
a Ori 5 50 7.4,0.9R Ma] (M1) (3) 1 I 
6 Aur 551 54.3] 3.88 Ko | G7 (2)}o} | 
Gem 638 25.2) 3.18 Gs | G7+ (4) | 3 I 
Gem 640 13.0| 3.40 Fs | F5+ (3)|1] I | 
« Gem 720 28.0] 3.89 Ko | Got (3) ] 1 | 
aCMi | 734 5.5}0.48 Fs | Fat V 
« Gem 738 24.6] 3.68 Gs | G7 (2) | 2 Ill | —- 
Gem 739 28.3] 1.21 Ko | Ko (6)| 3] II |- 
8 Cnc 811 9.5| 3-76 K2 | K4 (3)|o] WI |+ 
o UMa 822 61.1| 3.47. Go | G7 (6) | 2 HI | —- 
Hya 8 42 6.8} 3.482 F8 | Fo (2) | 1 It | 
¢ Hya 8 50 6.3] 3.30 Ko | Go (4)|2| II |- 
JHR 3579 | 854 42.2] 4.09 Fs | Fst (3)]/1] V 
a Lyn 915 34.8] 3.30 Ks | (4) | 1 Ill 
6UMa 9 260 52.1| 3.26 F8 | Fo+ (2) 1 Ill 
11 LMi 9 30 5. 48 ~-Ko | G8+ (2) 1 V 
Leo gh4o™+24°2| 3.12 Go | Gr (4) | 2 Il 


467 
No. CR w 
0.45 (1) 
.07 (2) 
.00 (1) 
.37 (2) 
31 
.26 (1) 
.33 (8) 
9: (1) 
.00 (1) 
-39 (1) 
13. (1) 
“29 (1) 
-35 (2) 
.40 (2) 
17.. .20 (2) 
41 (2) 
(1) 
(1) 
(1) 
22... “25 (2) 
.44 (1) 
-23 (1) 
.44 (2) 
(1) 
27... .32 (1) 
36 (1) 
29... -34 (5) 
30%... (3) 
.15 (8) 
.41 (1) 
(2) 
-43 (1) 
35° .13 (1) 
36... .09 (3) 
(1) 
(1) 
(1) 
er .07 (1) 
43 (2) 
.40 (2) 
44... (1) 
45... .22 (3) 
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No. Star a1goo 51900] my Spy | Spy wi{R Ly CR w 
46. 40 Leo 10'14™+4+20°o| 4.97. Fs | F7+ (2) | 1 V —o.37 (1) 
47.- Leo 10 15 +20.4| 2.61R Ko | K2 (a) (1 
48.. UMa | 1016 +42.0] 3.21 Ks | (3) ] 1 IIT | + .25 (4) 
49.. 36 UMa | 10 24 +56.5| 4.84 F5 | Fo (2) | o V — .39 (4) 
50... 46 LMi 10 48 +34.8] 3.92 Ko | Kr (2) | o III | — .08 (1) 
a UMa | 1058 +62.3| 1.95 Ko | Ko+ (3)] 1 III | — (6) 
52 y UMa | 11 4 +45.0] 3.15 Ko | Kr (2) | o Ill | — .03 () 
Te ¢ UMa | 11 13: +32.1| 4.418 Go | Fot (4) | 1 V — .40 (1) 
UMa | 11 13 +33-6| 3-71 Ko | K3+ (2) 1 Ill | + .16 (1) 
2... x UMa | 11 41 +48.3| 3.85 Ko | Kot (3) | 1 Ill | + .03 (1) 
a6... 6 Vir 11 46 + 2.3] 3.80 F8 | F8 (6) | 2 V — .37 (2) 
B CVn 12 29 +41.9| 4.32 Go | Go V — .36 (2) 
ES. e Vir 1257 +11.5| 2.95 Ko | G8+ (5) | 3 II | — .14 (1) 
a9. .< BCom | 13 7 +28.4} 4.32 Go | Fo (4) | 2 V — .36 (2) 
60... t Boo 13 43 +18.0] 4.51 F5 | F6 (3) | o V — .41 (1) 
61. n Boo 13 50 +18.9| 2.80 Go | Fo Cynic III | — .33 (2) 
62. a Boo 14 11 +19.7| 0.24 Ko | K1r+ (6)] 1 III | + .05 (7) 
63. 6 Boo 14 22 +52.3| 4.06 F8 | F7+ (4) | 2 Vv — .41 (2) 
64. p Boo 14 28 +30.8| 3.78 Ko | K2 (2)|o] I | + .o9 (1) 
65. | 1447 +19.5| 4.048 G5 | Ko 1} V | — -25 G) 
66.. BUMi | 1451 +74.6| 2.24 Ks | K4 (3)} 0] IIL | + .20 (2) 
67.. B Boo 14 58 +40.8| 3.63 G5 | G7 (5) | 2 Wt | — .12 (2) 
68.. 5 Boo 15 12 +33-7| 3-54 Ko | Got (2) | 1 Ill | — .14 (3) 
69.. Dra 15 23 +59-3| 3-47 Ko | Ki+ (3) | 1 III | + .02 (4) 
70. a Ser 15 39 + 6.7| 2.75 Ko | Ki+ (2)] 1 IIL | + .03 (2) 
Ser 15 52 +16.0] 3.86 Fs | Fo+ (6) | 4 V — .43 (2) 
72.. 6 Oph 16 9 — 3.4| 3.03 Ma] Mo (3)]1 Il | + .27 (1) 
73.. n Dra 16 23 +61.7| 2.89 Gs | Got (2) | 1 Il | — .17 (4) 
74.. 6B Her 16 26 +21.7| 2.81 Ko | G8+ (2) | 3 IIL | — .15 (2) 
75.. ¢ Her 16 38 +31.8| 3.00 Go | Gr (6) | 4 V — .32 (1) 
96... n Her 16 40 +39.1| 3.61 Ko | G6 (4)|6| Il | — .17 Gd) 
x Oph 16 53 + 9.5| 3-42 Ko | Ki+ (1) ]... III | — .o1 (1) 
78... a Her 17 12 +36.9| 3.36 Ks | K3 (3) | 2 | II-III} + .16 (1) 
79.. B Dra 17 28 +52.4| 2.99 Go | Fo Cay sha I — .18 (2) 
8B Oph 17 39 + 2.94 Ko | K2 (4) | + .03 (1) 
81. uw Her 17 43 +27.8| 3.48 Gs | Gst+ (3) | 2 V — .26 (2) 
82. Dra 17 52 +56.9| 3.90 Ko | Ki (2)|o} | — .o4( 1) 
83. 6 Her 17 53 +37.3| 3-99 Ko | Ko+ (2) | 1 II | + .12 (1) 
84. & Her 17 54 +29.3| 3.82 Ko | G8 (4) | 2 | WI-IV| — .15 (1) 
85. y Dra 1754 +51.5| 2.42 Ks | K5 (5) | 1 II | + .22 (2) 
86....| 70 Oph 18 o + 2.5| 4.078 Ko | Ko+ (2) | 1 V — .20 (1) 
87... n Ser 18 16 — 2.9| 3.42 Ko | Go (3) | 3 IV | — .15 (1) 
88... x Dra 18 23 +72.7| 3.09 F8 | Fo (s) | 1 V — .40 (2) 
8g....| 110 Her 18 41 +20.5| 4.20 F5 | F6 — .43 (1) 
go.. 6 Dra 1913™+67°5| 3.24 Ko | Ko (2)| 2] | —0.13 (1) 


| 
| 
| 
| 
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TABLE 4—Continued 


No. Star a1900 61900 | myy Spy | Spy w Ly Cp w 
Ql. x Cyg 19'15™+53°2| 3.98 Ko | G8+ (3) | 2 | III-IV| —o.15 (1) 
92. o Dra 19 33 +69.5| 4.78 Ko | G8+ (4) | 3 — .23 (2) 
93: 9Cyg | 19 34 +50.0| 4.64 | F4 (2)| 2} Vi |— .44 
Aql 19 42 +10.4| 2.80 K2|K2+ (4) | 1 II + .25 (4) 
05. e Dra 19 49 +70.0] 3.99 Ko | G7 (2) | 2 Ill | — .21 (1) 
96.. B Aql 19 50 + 6.2] 3.90 Ko | Ko (5s) | 2| IV-V | —.ar (1) 
Q7-...| ySge 19 54 +19.2| 3.71 Ks | Ks+ (4) | 1 | | + .25 (1) 
98... vy Cyg 20 19 +39.9| 2.32 F8 | F8 (5) | 2 Ib — .30 (6) 
99.. B Del 20 33 +14.3| 3.72 Fs | F4 (2) | 2 III | — .40 (1) 
100... e Cyg 20 42 +33.6| 2.64 Ko | Ko (4) | o | IN-IV| — .o9 (2) 
101. nCep | 20 43 +61.5| 3.59 Ko | Got (5) | 3 | IV-V | — .17 (1) 
102. ECyg | 21 1 +43.5| 3.92, Ks | K3+_ (4) | 1 I | + .30 (1) 
103. 61 Cyg 21 2 +38.3| 5.578 Ks | K5+R (2) V + .04 (2) 
104. ¢ Cyg 21 g +29.8| 3.40 Ko | Go (2) II — .10 (1) 
105. Peg 21 39 + 2.54 Ko | Ki+ (6) | 1 I + .25 (1) 
106. a Aqr 22 1 — 0.8] 3.19 Go | Fo (4) | 2 I — .16 (1) 
107. . Peg 22 2 +24.9) 3.96 F5 | F4 (4) | 1 Vv — .45 (1) 
108. ¢ Cep 22 7 +57.7| 3.62 Ko | Ki (6) | 1 I + .24 (1) 
109... n Peg 22 38 +29.7| 3.10 Go | (4)/ 4] HI | — (2) 
110.. Peg 22 45 +24.1| 3.67 Ko | G7+ (2) 1 III | — .18 (1) 
0°: re « Cep 22 46 +65.7| 3.68 Ko | Ko (2) | o III | — .10 (2) 
112....|HR 8832 23 9 +56.6| 5.65 K2 | K3 (2) | o — .08 (1) 
Psc 23 12 + 2.7| 3.85 Ko | G7 | 3 III | — .20 (1) 
Cep 23°35™+77°1| 3.42 Ko | Ki (3) | | —o.10 (1) 

NOTES TO TABLE 4 

7. a@UMi_ The light and color index are both variable to a slight degree. 

25. ¢Aur. Visual magnitude at maximum. 

26. aAur. The apparent magnitude of each of the two components is considered 
to be 1.0. 

29. aOri. Semiregular variable. 

39. ¢€Hya. Apparent magnitude of primary; secondary is probably too faint to 
affect the color or spectral type. 

47. yLeo. Apparent magnitude of primary; secondary may affect spectral type 
and color to a slight degree. 

53. &€UMa. The two brightest components of this multiple system contribute to 
both the observed spectral type and the color index. The apparent 
magnitude is that of the brighter star. 

65. €Boo. The fainter component probably has little effect on the observed spec- 
tral type and color index. 

86. 70Oph. ‘The fainter component may affect the observed spectral type and color 
index slightly. 

103. 61 Cyg. The apparent magnitude of the brighter component is given. As the 


photoelectric color index applies to both stars considered as one, the 
catalogue type is a mean value weighted according to the apparent 
magnitudes of the components. 


Stars having composite spectra have been omitted. 
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what greater than those computed from the internal agreement. 
The values of the mean error have correspondingly been rounded off 
arbitrarily in the following manner: 


Adjusted mean error of one plate...... +10 $2.9 to. 


It is considerably more difficult to give a measure of the accuracy 
of the luminosity classification because of the fact that the mean 
error varies with spectral type and also with luminosity. Because of 
the many marked spectral characteristics of the supergiants, the 
probability is exceedingly small that a good spectrogram of a super- 
giant will be classified in a group of lower luminosity. At the other 
extreme is the case of the separation of the subgiants and dwarfs of 
types G8—K1; the mean error of one plate here is of the order of one- 
half of a luminosity class. A rough set of estimates of the mean error 
of the luminosity classification between successive groups for one 
good spectrogram, in units of a luminosity class, is shown in the ac- 


companying table. 


Class F4-F7 F8-G1 G2-G6 G7-K1r K2-K5 


The mean error of a single observation of photoelectric color index 
for stars in Table 4 is +0.016 mag.; this value was determined by 
Schlesinger’s method from stars which were observed between three 
and six times. The average mean error for a tabular color index when 
the number of observations is taken into account is +0.013 mag. 

4. The photoelectric color indices of the three most important 
luminosity groups are plotted against Yerkes spectral types in Fig- 
ure 1. In addition to the stars in Table 4 the colors of the brightest 
B and A stars are plotted against their Draper spectral types; all 
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B stars brighter than apparent magnitude 2.0 and all A stars brighter 
than 3.0 for which Bottlinger measured color indices have been in- 
cluded. The supergiants (which are denoted by open circles) show 
little or no systematic difference in color from main-sequence stars 
(closed circles) for types earlier than F5; after this they diverge rap- 
idly in color and reach a maximum difference of more than four- 
tenths of a magnitude on Bottlinger’s scale at G8—Ko. The diver- 
gence in color of the supergiants from the normal giants (half- 
closed circles) also passes through a maximum in the nighborhood 
of types G8-Ko and decreases rapidly after that point; by the time 
Ks is reached the difference is relatively small. 

A general inspection of Figure 1 shows the smooth progression of 
change in color index with spectral type for each luminosity group; 
the dispersion in color within a luminosity class at any given spectral 
type is low, and it is noteworthy that there is not a single case, for 
stars of types later than F8, of a dwarf which is as red as a giant of 
the same type or of a giant which is as red as a supergiant. Of par- 
ticular interest is the low dispersion in color for dwarfs of types 
F4—Gr and for normal giants of types K4—Mo. It should not be as- 
sumed, however, that the actual dispersion in color is necessarily 
higher for stars having the same absolute magnitude and spectral 
type in the region G6—K1; the scatter of the plotted points is con- 
siderably greater, but this is due, at least to a considerable extent, to 
the large variation of color with absolute magnitude for stars of 
these spectral types, to the somewhat lower accuracy of the classi- 
fication at Gs5s—Go, and to the more rapid change of color with spec- 
tral type. This tends to superpose an appreciable additional acci- 
dental error on the average errors of classification of spectral type 
and of the color indices. For the F stars the total range in color with 
luminosity is much less; in the case of the K4—-Mo stars the problem 
is simplified by the absence of subgiants and by the small difference 
in color between normal giants and supergiants. 

For these reasons a test of the accuracy with which the intrinsic 
color index is defined by spectral type and luminosity is most readily 
made for the F4—G1 stars and for the giants of types K4—Mo. Be- 
cause of the considerably greater number of F4—G1 main-sequence 
stars, they form the most suitable group for investigation. The color 
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indices and spectral types of these objects are plotted on a larger 
scale in Figure 2. The radius of the plotted circles corresponds to 
the average mean error of the color indices. 

The color indices are well represented by a straight line whose 
equation is 


Cp = —o™429 + 07015(S — 


where Cp is the color index of a main-sequence star on Bottlinger’s 
scale, S is the Yerkes spectral type, and (S — F5) is expressed in 


—0.45 


.40 


—0.30 


F4 F5 Fo F7 F8 Fo Go G1 G2 
Fic. 2.—Photoelectric color indices and Yerkes spectral types of main-sequence stars 
of types F4-Gr. The straight line has the equation: 
Cp = —o™429 + oMors (S—Fs). 


units of tenths of a spectral class. The mean error of a plotted point 
from this straight line is +0.018 mag. Of the twenty-three points 
plotted, two have the maximum residual of 0.04 mag. The repre- 
sentation of the observed colors of the main-sequence stars by their 
spectral types is, therefore, very close. 

The average mean error of a value of Yerkes spectral type for the 
stars in Figure 2is +0.6S. As the change of color index is 0.014 mag. 
per tenth of a spectral class, the mean error of classification is equiv- 
alent to a mean error of +0.008 mag. in color index. The mean error 
of a catalogue value of the color index is +0.013 mag. This corre- 
sponds to a combined mean error of +0.015 mag., as compared with 


Me 

OO 
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an observed mean error from Figure 2 of +o.o19 mag. The differ- 
ence in the two values is slight; when it is considered that corrections 
for components of spectroscopic binaries and close visual binaries 
have not been made and that two systems of internal mean errors 
are being combined, the agreement can be considered to be satis- 
factory. The conclusion seems justified that within the errors of 
observation there is no real dispersion in the intrinsic colors of the 
F4-Gr1 main sequence stars. The color indices of these objects seem 
to be uniquely defined by their spectral types. 


TABLE 5 
COLORS OF Gg-Ko NORMAL GIANTS 
Star Verkes Color Index A 
Type 
Go (4) —o.07 (1) 0.05 
Got (3) .12 (2) 00 
Ko (2) (2) 02 
Ko (7) —o.og (3) 0.03 


The dispersion is correspondingly low in the case of the K4~Mo 
normal giants. The subgtants, which introduce additional uncer- 
tainty in the neighborhood of Ko, are absent, and the supergiants 
are only slightly redder than the normal giants. 

We shall finally consider the normal giants of spectral types Gg- 
Ko to obtain an approximate value of the dispersion in color at an 
unfavorable position, that is, where the change of color with lumi- 
nosity is rapid. The stars are listed in Table 5, together with the de- 
viation of their color indices from the arithmetical mean. The 
mean error of an individual color index in Table 5 is +0.03 mag.— 
about 50 per cent higher than in the case of the F4—Gr main-se- 
quence stars. The increase in the mean error can be accounted for 
by two factors: (1) the change of color index with spectral type is 
more rapid for K giants than for the F5—Go dwarfs; (2) the rapid 
change of color index with luminosity introduces an additional de- 
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gree of uncertainty. The first factor is the more important, as the 
value of dC/dS for normal giants at Ko is about 0.04 mag. 

5. The evidence as a whole, therefore, seems to point to the con- 
clusion reached for the F4—Gr dwarfs: if the accuracy of the obser- 
vations is taken into account, the color indices of the stars of spectral 
types F4-Mo are uniquely determined by their spectral type and 
luminosity. This conclusion makes way for several interesting prac- 
tical applications of calibration diagrams similar to Figures 1 and 2. 
If, for near-by stars, or stars in regions known to be free from ab- 
sorbing clouds, spectral types and luminosities are available of the 
accuracy of those in the present paper, then the color indices of these 
stars may be predicted with considerable accuracy. In the case of 
the Gs5—Kr stars these color indices should be of the same order of 
accuracy as those obtained from a comparison of photographic and 
photovisual catalogues of high weight; for the F4-G1 dwarfs and 
for the K4-Mo giants the accuracy of the predicted colors should be 
somewhat higher. The conclusion is especially important in the case 
of the F5-Go dwarfs; because of their low absolute magnitude, diffi- 
culties with selective absorption should not be encountered until 
stars fainter than the eighth apparent magnitude are considered. 
They are thus capable of serving as “‘calibrated lamps’’ whose spec- 
tral energy distribution may be predicted with high accuracy and 
which can be used with considerable confidence as standards of color. 


Some of the spectrograms were obtained by Dr. P. C. Keenan, 
Mr. E. Ebbighausen, and Mr. G. W. Wares. I am indebted to Miss 
Edith Kellman for computational assistance and for making the 
drawings. 
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IONIZATION AND RECOMBINATION IN THE THEORY 
OF STELLAR ABSORPTION LINES AND 
NEBULAR LUMINOSITY 


S. CHANDRASEKHAR 


ABSTRACT 

In this paper the effect of collisions and photo-ionizations on monochromatic radia- 
tive equilibrium is considered. A formalism is developed which includes the theory of 
stellar absorption lines (Strémgren) and the theory of nebular luminosity (Zanstra) as 
two limiting cases. A second approximation for the nebular problem is also sketched. 
The mathematical analysis is devoted to the treatment of the three-state problem in- 
cluding collisions. When collisions can be neglected, the approximate equations of trans- 
fer used admit of two flux integrals and two K-integrals. 

1. Introduction.—The theory of the formation of stellar absorp- 
tion lines hinges on the problem of writing down the appropriate 
equation of radiative transfer. If we can write down the equation of 
transfer, then we have, in principle, solved the problem. But in the 
considerations leading to the equation of transfer, we can follow one 
of two methods: either argue in terms of macroscopic notions, e.g., 
the mass absorption and emission coefficients, etc., or in terms of mi- 
croscopic notions, e.g., the Einstein transition probabilities, effec- 
tive cross-section for captures, etc. The two methods, if properly 
handled, should of course lead to identical results; but a comparison 
of the two methods is often useful. Thus, in the so-called ‘‘combined 
Schuster-Schwarzschild problem’’ it is assumed that we can charac- 
terize the material by a mass absorption coefficient, x, and a mass 
scattering coefficient, /. The equation of transfer is then written 


down in the form 


dl, 


pds (1) 


where /, is the specific intensity at frequency v, B, the Planck func- 


tion, and 
476 
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In writing down (1), it is assumed that the emission consists of two 
parts: (a) the Kirchhoff emission, «,B,, and (6) the emission due to 
scattering, of amount /,/,. On the other hand, as E. A. Milne" has 
shown, we can derive an equation similar to (1) by an entirely dif- 
ferent type of discussion, which does not involve an appeal to Kirch- 
hoff’s law. The assumptions, (a@) a Maxwellian distribution of veloc- 
ities and (b) the effect of collisions in exciting and de-exciting the 
atoms, lead to the equation of transfer 


dI,,, J, + 
where dt,,, is the optical depth measured in the direction of s, 
I 
Bi, B,,, ’ (4) 


where 6,, and B,, are the probabilities of an atom in state 1 being 
excited to state 2 by collision or by the absorption of radiation of ap- 
propriate frequency. If we write 


Ky 
iz = L (5) 
and 
dt, = (ky + l)pds , (5’) 


we see that (1) and (3) are formally equivalent. We shall return to 
the discussion of the physical meaning of (5) and (5’); but it is al- 
ready clear that Milne’s method, if properly interpreted, will remove 
some of the obscurity involved in the appeal to Kirchhoff’s law in 
the usual treatment. 

However, the theory of the formation of stellar absorption lines 
in the frame either of (1) or of (3) encounters difficulties—among 
others, the matter of the central intensities. 

The physical foundations have been analyzed by B. Strémgren,’ 


™M.N., 88, 493, 1928. 
2Zs. f. Ap., 10, 237, 1935; Handbuch der Ap. 7, 221-235, 1936. 
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who incorporates into the theory at an early stage the influence of 
electron captures and photo-ionizations. Strémgren’s theory is en- 
tirely satisfactory, but at the same time it would be of some interest 
to examine whether Milne’s method can be suitably generalized to 
lead to results equivalent to those obtained by Strémgren’s method. 
We shall see that this is possible. 

There is another aspect to this problem. As is well known, Zan- 
stra’s theory of nebular luminosity utilizes the processes of electron 
captures and ionization as a means for the conversion of the ultra- 
violet energy beyond the head of the Lyman series in the incident 
starlight into energy mainly in Lyman a. The problem of radiative 
transfer associated with the Zanstra theory has been studied by 
V. A. Ambarzumian and the writer.’ The question now arises wheth- 
er we cannot construct a scheme which will include, as limiting cases, 
Strémgren’s equation of transfer at one end and the equations of ra- 
diative transfer for the Zanstra theory at the other end. The main 
object of this paper is to supply such a scheme. 

2. Definitions and fundamental equations.—The method consists 
in extending Milne’s procedure to include cyclical transitions be- 
tween three states, one of which we eventually identify with the con- 
tinuum. 

First let us consider the case of three discrete levels of statistical 
weights, g., g2, g;. We introduce the Einstein coefficients Ann, Bum, 
Bmn, defined with respect fo the intensity of radiation. We have the 
usual relations 


Anm Sm 

(6) 
Bam _ &m 
Bus Bn (7) 


We shall use the abbreviation 


2hv3.n 
(8) 


omn = 


3§. Chandrasekhar, Zs. f. Ap., 9, 266, 1935. 
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If amn(v) is the atomic absorption coefficient for the atom in the 
state m for radiation of frequency v, then 


the integral on the left-hand side being taken through the narrow 
range of frequencies about v,,, within which the atom can absorb. 
We write (9) in the form 


= 40 ’ (10) 


where a,,, is the mean atomic absorption coefficient in the interval 
A’mn. We shall next introduce the probability coefficients 
to take into account the effect of collisions in exciting and de-exciting 
the atoms in the states m and n, (m < n)—more precisely, the prob- 
ability that in unit time an atom in state m, (m), is excited (de-excit- 
ed) by collisions of the first (second) kind is Bmn, (nm). If we as- 
sume that a Maxwellian distribution of the atomic velocities corre- 
sponding to a certain temperature 7 is maintained, then we should 
have 

bmn — /kT 


Now, in a steady state the number of atoms, 1,, ”., 23, in each of 
the states must remain constant. Thus, the constancy of the num- 


ber of atoms in the ground state 1 yields 


(BrJ 12 + by + Bi 3J + b,;) n(A 21 + + 
= n;(Aj + ByJi3 + = 0, (12) 


where the J’s are used to denote 


dw 1 
Inn = Ef sin 6d6 . (13) 


q 
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Using (6), (7), (8), and (11), we can re-write (12) in the form 


§2 12 


My (Br + biz + + — + Ju + | 
5 (14) 


Similarly, the constancy of the number of atoms in state 3 gives 


3 13 
(15) 
+ 4 Ba (vs J x + B., =o. | 


As may easily be verified, (14) and (15) automatically provide for 
the constancy of the number of atoms in the second state. We intro- 


duce the abbreviations 


= ehvmn/kT (m n) ’ (17) 
= Omn + + Fun (18) 


Equations (14) and (15) now take the forms 
ty { 13 + @12) + + @1;)} 


— n> — n;~ Bit = 0, I 
3g, 3613 (19) 


13 + + 23 + @23) 


Br313 + Busts =o, (20) 
§3 §3 


which are our equations of radiative equilibrium. 

So far, we have restricted ourselves to three discrete states. We 
now wish to make state 3 the “‘continuum,”’ i.e., the transitions 
I > 3, 23, now correspond to photo-ionization. As Ambarzu- 
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mian‘ has shown, we can treat the continuum on the same footing as 
a discrete state, provided we define for the weight g, the expression 


)3/? 


where g, is the statistical weight of the normal state of the ionized 
atom and N, is the number of free electrons per cubic centimeter. 
Furthermore, B,,and B,, are now no longer atomic constants but are 
the probability coefficients of the photoelectric transitions defined in 
an appropriate way. 

From (1g) and (20), we derive 


Ny 


£1 82 |, 
B,, 12 B, I + B, 2 + B, I B, Fs + 2 
36 3 3823 f 3813 3( 3 3) 


nN; 


B.,(J 2; + @23){ + + + @1;)} + Bik + @1;) 


3. The equations of transfer.—The equations of transfer are easily 
seen to be (cf. Milne, loc. cit.) 


= (x N2 + N20 12 (23) 
dy 
= (x. nN; I + N39 23 (25) 


4M.N., 95, 469, 1935. 


(22) 


15) 

= 
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where the a»,»’s are given by (10). The optical depths dt,., dt,;, and 
dt,, can be defined by 


Bm — Nn, 
dt,, — ny, = ds , (26) 


83 


dt 
47 


nN, 


ds . (28) 


82 
= — n ds = 
It should be noticed that our optical depths, dt,,,, are measured in 
the direction in which the variation of J is considered. The equations 
(23), (24), and (25) now take the forms: 


al me 
I 2 £2 
n 
dl, 
= Iu oon, (30) 
§3 
82 
dl, 
+ 923+ (31) 
N2—- nN 8 
* 


It will be noticed that (29), (30), and (31) involve the n’s only 
through the ratios (m,/n.), (2/3), respectively. The solu- 
tion (22) of the equations of radiative equilibrium provides the 
values of the required ratios in terms of the J’s; and equations (29), 
(30), and (31) form a system of simultaneous differential equations. 
To evaluate the terms involving the ’s we shall adopt a method 
of approximation introduced by Ambarzumian. We regard J,./¢:. 
and J,;/¢.; aS quantities of the first order of smallness, while we 


| = 
| | 


IONIZATION AND RECOMBINATION 483 


regard J,,/0,; and J,,J,;/0:.0.; aS quantities of the second order of 
smallness. The justification for this is found in the circumstance 
that, near the star, J,,/0,, and J,,/0,; are of the order of magnitude 
and respectively, while J,;/¢,;, for instance, is 
of the order of magnitude et)" | On the other hand, if we 
are far away from the star, the effect of the dilution factor will be 
to make our approximation sufficient for practical purposes. Let 


= NnnBy,,, 5 (32) 
where B, is the Planck function: 
Brinn = — , (33) 
By (18) and (32) we have 


Hence, we have 


= o1,(1 + + first-order quantities , (35) 
£.; = 0.3;(1 + m3) + first-order quantities , (36) 
£1; = 013(1 + m3) + second-order quantities . (37) 


We evaluate the expressions occurring on the right-hand sides of the 
equations of transfer, using (22) consistent with our method of ap- 
proximation. After some rather lengthy calculations, it is found that 


£2 
n, 
Bas + bay Miz + 
nN, — Ny; 
§3 
£2 
83 M13 M13 Ji; 
Mr + + M23 Jz (4 ) 
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where 


Let 


Then 
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= + M12) (41) 
Mr = + m5) , (42) 
a3 = 82B.302;(1 + 723) - (43) 

J + m2B,,, 
Ju = (44) 

Ji3 + 

3 113) (45) 

+ 123B,, 
J. = 6 
+ 123) “4 ) 

+ ba; Miz M23 (47) 
M23) (1 p) Mi2 (47 ) 


By (46) and (47) we can re-write (38), (39), and (40) in the forms 


= + (1 = P) 


= (1 p) + (1 p) , (48) 


(49) 


(50) 


| 
|| 
p 
| 
£2 
Ny — 
§3 
= 
N2 N3 
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The equations of transfer now reduce to 


: 

= — + + — (52) 
ihe ( ad ) J + 13 ( 
dt,, Pors da 53) 


We shall now proceed to consider two limiting cases of the foregoing 


system of equations. 
4. The formation of absorption lines: comparison with Strimgren’s 


results.—According to equations (44) and (51), 


where we have used e€,, to denote 
We re-write (54) in the form 
(x + m2) = — (1+ t+ + 
+ €:2012(1 + 2) Ju (56) 


Jas 


or, again, after some further rearrangement of the terms, 


(1 =- (1 + M2) 112 + (1 = 


+ + €20nB,,, (57) 


where we have used Q,, to denote 


12 12/313 2 


od 
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The equation of transfer (57), which we have derived, includes 
the effect of collisions and of the cyclical transitions (resulting from 
ionization and captures) on monochromatic radiative equilibrium. 
Consequently, it can form the basis of an analysis of the formation of 
absorption lines, and as such we can regard 7 as small, compared 
with unity. This simplifies our expression for Q,,. We can write, ac- 
cording to (45), (46), and (58), 


O12 J 13 O23 


(sg) 
013 Fa 5 
Further, we have 
013 023 


where we interpret 7. as the color temperature of the radiation. On 
the other hand, the Planck function B,,, is defined with respect to 
the temperature of the material. Hence, we have 


We can also simplify our expression for e,,. According to equations 
(41), (42), (43), (47), and (55), we have 


8232302; * Ji; (62) 
~ B + £2B.;0;) O23 


Now, we can further assume that B,, « B,,, since the total probabil- 
ity of photo-ionizations from the excited state can be neglected, com- 
pared with the probability of photo-ionization from the ground state. 
We can thus write 


or, using (6), we have 
BJ, 
(64) 
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The term B,,/,,; has a simple physical meaning: it represents the 
probability of photo-ionization of atoms in state 2 owing to the in- 
cident radiation field. If we denote this probability by C.;, we have 


C, 


Now, Strémgren writes the equation of transfer in the form (Joc. cit., 
Eq. [37]) 


(ky + + (1 + kyBy + (66) 


We verify that Strémgren’s e and Q are identical with our e,, and 
Q,. (cf. our equations [65] and [61] with Strémgren’s [18] and the 
equation at the top of p. 253 of his paper). We thus see that our 
equation of transfer (57) will become identical with Strémgren’s (66) 
if we write in (57) 


dt,, = + l,,.) pds (67) 
= (68) 


We thus see that the correspondence between our treatment and 
Strémgren’s problem is of the same nature as the correspondence 
between Milne’s treatment and the combined Schuster-Schwarzschild 
problem. 

In comparing our method with Strémgren’s, the following points 
should be noted. 

1. In our method we consider (strictly speaking) only a one-con- 
stituent atmosphere, the relevant physical processes being collisions, 
photo-ionizations, recombinations, and radiative resonance transi- 
tions. Consequently, the only processes which give rise to general 
opacity in our model are collisions of the first and second kinds. It is 
in this sense that equation (67) has to be interpreted. 

2. Our method does not take into account the contributions to 
the general opacity arising from factors other than collisions and, 
further, entirely ignores the presence of other constituents in the at- 
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mosphere. In Strémgren’s method, however, it is assumed that if 
the general opacity, x,, owing to all factors and all the constituents 
of the atmosphere, is known, then the appropriate emission is 
x,B,(T), where T is the temperature of the material. 

3. The origin of the identity of our e,, and Q,, with Strémgren’s 
e and (Q is to be found in the circumstance that in both the treat- 
ments the effect of photo-ionizations on monochromatic equilibrium 
is considered in physically equivalent terms. 

It is thus clear that the fact that in the model case we have con- 
sidered—namely, the three-state problem including collisions—our 
method leads to an equation of transfer identical with that derived 
by Strémgren, can be interpreted as implying the validity of the 
form in which use is made of Kirchhoff’s law in Strémgren’s method. 

5. The equations of the planetary nebulae.—Returning to equations 
(51) and (52), we see that in the case of high dilution of radiation 
and in the absence of collisions we can write them in the forms 


where 
Equation (69) can be further simplified by using (10): 


These are precisely the equations of transfer in the form used by 
Ambarzumian and the writer. The probability » introduced in that 
theory is identical with our present /, defined according to (71). We 
thus see that our formalism developed in sections 2 and 3 of this 
paper includes the theory of stellar absorption lines and the theory 
of nebular luminosity as two limiting cases. 


i 
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6. The equations of the planetary nebulae: a second approximation. 
—In the absence of collisions, our equations of transfer are (cf. Eqs. 


[51], [52], and [53]): 


dl; 013 

012023 J; 

dts = + (1 p)J + p (75) 


where p is defined as in equation (71) and 


Mrs = = (76) 
Hence, 
Let 
y= (78) 
Then 
= 


It should be remembered that in the foregoing equations the optical 
depths are measured in the direction in which the variation of J is 
considered. 
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From our definitions of d/,., dt,,, dt.; we verify that 


<= 
dt; B,31;Av2, §3 (82) 
§3 


or, using (49) and (50) and remembering that for the case under 
consideration 7 = 0, we have 


dt;; By3v13 A023 O12 13 O12 (83) 


Hence, equation (81) can be re-written as 


dl, 20, 
O12 dt §2023 La3J 12 + (1 = 13} (84) 
13 


We now see that equations (79), (80), and (84) form a consistent 
system of equations. The right-hand side of equation (79) contains 
terms of the first and the second order of smallness, while the right 
sides of equations (80) and (84) contain terms only of the second 
order of smallness. Consequently, it would be sufficient to find the 
dominant term in dt,,/dt,,. We verify that, to the order of accuracy 
we are working with, we have 


We introduce the abbreviations 


| 
| 
| 
| 
| 
| 
| 
] 
| 
| 
| 
| 2 | I | 12 | 2 | 23 


IONIZATION AND RECOMBINATION 491 


Considering the case of material stratified in parallel planes and 
denoting dr to measure the optical depth in (1, 3) normal to the 
plane of stratification, the equations of transfer reduce to 


dI,, 
+ digi(t — p)Ji;, (88) 
cos qs; + + P)T ’ (89) 
COS ae + (1 a3 J} 12 + PY "J; (go) 


These equations admit of three first integrals. Multiply the equa- 
tions successively by 3 sin 6d6 and 3 sin 6 cos 6d@, and integrate from 
o to 7. Let 


Fan = af sin cos 6d6, (gt) 
Kan = sin? 6 cos 6d6. (92) 
Then 
Also 
(96) 
(97) 
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From (93), (94), and (95) we derive 


dF 12 


dr ’ (99) 
_ p dF,, 
— p) dr (100) 


We first notice that, according to (47), (76), (78), (83), and (86), 
we have 


p 
(101) 


Hence, (gg) and (100) yield the two flux integrals 
F,, + 6,F,; = Constant = S,. + 6,S;;, (102) 


F,; + 6.F,; = Constant = S,, + 6,S,; , (103) 


where we have used S,,, to denote the “incident” flux. From (102) 
and (103) we derive 


2 Sx) = Sas) . (104) 


We now derive a K-integral: from (g6) and (97), we have 


12 
+ 6; + 3) , (105) 
or, by (102), 
d 
Hence, 
+ = — + + constant . (107) 


We easily verify that (107) is equivalent to an expression which the 
author had derived in his paper on the radiative equilibrium of a 
planetary nebula (Joc. cit., Eq. [87]). However, we have now derived 
it under more general conditions. 


| | 
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It is possible to derive another approximate K-integral if 


= . (108) 


Then, according to (104), 6./,, = 6,/.;, and we can re-write (98) as 


dK, 
or, by (96), 
61920 12 6241 Jia (1 10) 


If we make the approximation 


Ju = 3Ku, (111) 
we have 
or 
CK, = : [K?, + constant] , (113) 


which can be written alternatively as 


Kas _ 3 + constant (114) 


023 6192 023 O12 


Since (K,,/¢,2)? is a quantity of the second order of smallness, while 
(K.;,/0.2;) is one of the first order of smallness, it is apparent from 
(114) that, if we neglect quantities of the second order, then K,, is 
a constant. 

The integrals (102), (103), (107), and (113) have been derived 
from the equations of transfer (88), (89), and (go), which are ap- 
proximations to the exact equations of the three-state problem. 
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However, Henyey has since shown that integrals corresponding to 
those we have derived can be isolated for the exact problem. 

We shall next consider a little more closely the case S,, = S.; = 0, 
so that the integral (113) is valid. Let S,, be the incident flux in the 
(1, 3) radiation. We can then write 


Fas = + (115) 


K:, = + Ji3 = + (116) 


where §,,; and &,, are defined for the scattered radiation. Our flux 
integrals can be written as 


= 6,6.[S1;(1 — — = . (117) 
Furthermore, we have 
+ 6:K,; = — 46,S:;(7 + C,) , (118) 
where C, is a constant. From (93) and (96) we derive 


qi dr? = PY [VI 12 (119) 


We shall use the approximation 


= 3K = 3K2; = J (120) 
Then 


Using the K-integral in the form (118), the foregoing equation re- 
duces to 


44: 


X (r+ 3e7+C,). (122) 


dr? 


| 
| 
| 
| 
| 
| 
| 
| 
| 
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Using our second K-integral and denoting by C, the constant of in- 
tegration occurring in (113), we finally have 


2 923 O12 


44: 

which is the fundamental differential equation of the problem. After 

solving (123) for K,., we determine F,, by using (96). Once K,, and 

F,, are known, all the other quantities can be derived successively, 

using the different integrals. Concrete applications of (123) will be 

found in a forthcoming paper. 
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REVIEWS 


Introduction to Physiological Optics. By J. P. C. SourHaty. New York: 

Oxford University Press, 1937. Pp. 426. $5.50. 

This very readable introduction to physiological optics is intended as 
much for laymen who wish to know something about vision as for stu- 
dents who need a broad treatment of a subject finding application in 
physics, medicine, optometry, ophthalmology, cinema technique, and il- 
lumination. Scientific workers who obtain data with visual instruments 
may find the book advantageous for recalling some of the vagaries of the 
eye. What little mathematics is used involves scarcely anything more 
difficult than the algebra of first-order optics and of colorimetry. All con- 
clusions are stated clearly in words instead of symbolic equations. Con- 
siderable interest is provided by the frequent description of little experi- 
ments designed to show such points as the adaptation and accommoda- 
tion of the eye, blur circles on the retina, the Purkinje effect, the blind 
spot, and various aspects of binocular vision, of rod versus cone vision, 
of contrast effects, and of color perception. A number of the classical 
experiments of Scheiner, Wollaston, Helmholtz, and others are included 
among these. The text is often varied by the introduction of appropriate 
passages from earlier books and papers by many authors. There is a 
profusion of drawings. 

Chapter i describes the structure of the eye and some of the subjective 
phenomena of vision, including Fechner’s law, which is here stated to 
be “more honored in the breach than in the observance.’’ Chapter ii is 
a discussion of the imagery of the Gullstrand schematic eye. The next 
three chapters deal with the problem of correcting eye-glasses for the 
mobile eye afflicted with ametropia, presbyopia, or astigmatism. This 
section includes a discussion of the design of the Punktal lens, and has 
several pages on contact eyeglasses. Chapter vi treats in detail of the 
movements of the eyes in binocular vision, of depth perception, and of 
stereoscopy. Rod and cone vision are compared in the next chapter, 
which deals with dark adaptation, color sensitivity, and zonal sensitivity 
of the retina. Chapters viii and ix are devoted to color vision and color 
mixing. The Young-Helmholtz theory of color vision is described and 
criticized, and the theories of Hering and Ostwald are presented. The last 
chapter discusses persistence of vision, flicker phenomena, and after- 
images in monotone and in color. 

Although the topics discussed are essentially the standard ones of 
physiological optics, the author’s treatment combines an exceedingly in- 
teresting approach with judgment in selecting unhackneyed experiments 


and illustrations. Ci. Mel 
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